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PREFACE 


This  report  describes  the  work  performed  under  the  DOT/FAA  High-Velocity 
Jet  Noise  Source  Location  and  Reduction  Program  (Contract  DOT-OS-30034). 

•  Investigation,  including  scaling  effects,  of  the  aerodynamic 
and  acoustic  mechanisms  of  various  jet  noise  suppressors. 

•  Analytical  and  experimental  studies  of  the  acoustic  source 
distribution  in  such  suppressors,  including  identification 
of  source  location,  nature,  and  strength  and  noise  reduction 
potential . 

•  Investigation  of  in-flight  effects  on  the  aerodynamic  and 
acoustic  performance  of  these  suppressors. 

The  results  of  these  investigations  led  to  the  preparation  of  a  design  guide 
report  for  predicting  the  overall  characteristics  of  suppressor  concepts, 
from  models  to  full  scale,  from  static  to  in-flight  conditions,  as  well  as 
a  quantitative  and  qualitative  prediction  of  the  phenomena  involved. 

The  work  effort  in  this  program  was  organized  under  the  following  major 
Tasks,  each  of  which  is  reported  in  a  separate  Final  Report: 


Task  1  -  Activation  of  Facilities  and  Validation  of  Source  Location 
Techniques. 

Task  2  -  Theoretical  Developments  and  Basic  Experiments. 

Task  3  -  Experimental  Investigation  of  Suppression  Principles. 

Task  4  -  Development  and  Evaluation  of  Techniques  for  In-Flight 
Investigation. 

Task  5  -  Investigation  of  In-Flight  Aeroacoustic  Effects  on 
Suppressed  Exhausts. 

Task  6  -  Preparation  of  Noise  Abatement  Nozzle  Design  Guide  Report. 

Task  1  was  an  investigative  and  survey  effort  designed  to  identify  acous¬ 
tic  facilities  and  test  methods  best  suited  to  jet  noise  studies. 

Task  2  was  a  theoretical  effort  complemented  by  theory  verification  ex¬ 
periments  which  extended  across  the  entire  contract  period  of  performance. 


Task  3  represented  a  substantial  contract  effort  to  gather  various  test 
data  on  a  wide  range  of  high-velocity  jet  noise  suppressors.  These  data,  in¬ 
tended  to  help  identify  five  optimum  nozzles  for  in-flight  testing  in  Task  5, 
provided  an  extensive  high  quality  data  bank  useful  to  the  preparation  of  the 
Task  6  design  guide  as  well  as  for  future  studies. 


Task  4  was  similar  to  Task  1,  except  that  it  dealt  with  the  specific 
test  facility  requirements,  measurement  techniques,  and  analytical  methods 
necessary  to  evaluate  the  in-flight  noise  characteristics  of  simple  and 
complex  suppressor  nozzles.  This  effort  provided  the  capability  to  conduct 
the  flight  effects  test  program  of  Task  5. 

Task  6  embodies  the  salient  results  of  Task  2,  3,  4  and  5,  and  combines 
them  with  other  contractor  results  into  a  noise  abatement  nozzle  design 
guide  which  permits  acoustic  and  performance  prediction  of  future  high-speed 
engine-suppressor  installations. 

The  present  volume,  a  supplement  to  the  design  guide,  documents  two  jet 
noise  prediction  methods  developed  under  the  contract:  the  engineering  cor¬ 
relation  of  (M*S)  model  and  the  unified  aeroacoustic  model  (M*G*B)  (each 
capable  of  accounting  for  flight  effects).  The  objective  of  this  report  is 
to  provide  users  with  a  description  of  the  methods  and  associated  computa¬ 
tional  procedures  in  sufficient  detail  that  either  method  can  be  implemented 
and  utilized  as  a  useful  engineering  tool.  The  empirical  M*S  method  is 
capable  of  predicting  static  and  in-flight  acoustic  characteristics  of  multi¬ 
element  suppressors  applicable  to  both  advanced  turbojet  and  variable-cycle 
engines.  The  theoretically  based  M*G*B  method  is  capable  of  predicting  static 
and  in-flight  aerodynamic  and  acoustic  characteristics  of  jets  from  nozzles 
of  arbitrary  shape,  and  as  such  provides  more  insight  into  the  fundamental 
mechanisms  involved  in  a  given  configuration's  noise  signature. 


iv 


TABLE  OF  CONTENTS 


Section 

1.0 

2.0 

3.0 


4.0 


5.0 


Page 


SUMMARY  1 

1.1  The  Engineering  Correlation  (M*S)  Method  1 

1.2  The  Unified  Aeroacoustic  Prediction  (M*G*B)  Method  1 

INTRODUCTION  3 

ENGINEERING  CORRELATION  (M*S)  JET  NOISE  PREDICTION 

COMPUTER  PROGRAM  4 

3.1  Introduction  4 

3.2  Program  Nomenclature  4 

3.3  Description  of  Program  and  Subroutines  4 

3.4  Input  Description  37 

3.5  Output  Description  38 

3.6  Sample  Cases  54 

3.7  Program  Source  Code  Listing 

Section  3  References  7 * 


UNIFIED  AEROACOUSTIC  PREDICTION  MODEL  (M*G*B)  COMPUTER 


PROGRAM  96 

4.1  Introduction  96 

4.2  Program  Nomenclature  and  Symbol  Convention  97 

4.3  Description  of  Program  and  Subroutines  100 

4.3.1  Main  100 

4.3.2  ARCCOS(X)  100 

4.3.3  ERF (X)  100 

4.3.4  LSPFIT  112 

4.3.5  SLICE  112 

4.3.6  CRD  113 

4.3.7  OUTPUT  113 

4.3.8  SHOCK  115 

4.3.9  ATMOS  115 

4.3.10  PNLC  115 

4.4  Program  Usage  and  Logic  116 

4.5  Description  of  Input  118 

4.5.1  Notes  on  Input  126 

4.5.2  Example  Case  Input  Selection  132 

4.6  Output  Description  134 

4.7  Sample  Output  Listing  136 

4.8  Program  Source  Code  Listing 

Section  4  References  170 

CONCLUDING  REMARKS  206 


v 


LIST  OF  ILLUSTRATIONS 


Figure  Page 

3-1.  Computer  Program  Flow  Chart.  12 

3-2.  Nozzle  Types  Included  in  the  Correlation.  50 

3-3.  FORTRAN  Symbol  Convention  for  Acoustic  Arena  Variables.  51 

3- 4.  Definition  of  Cant  Angles  for  Multielement  Nozzles.  52 

4- 1.  FORTRAN  Symbol  Convention  for  Coordinates  and  Geometric 

Variables.  gg 

4-2.  FORTRAN  Symbol  Convention  for  Acoustic  Arena  Variables.  99 

4-3.  Computer  Program  Flow  Chart.  109 

4-4.  Possible  Solution  Types  for  a  Maximum  of  Two  Turning 

Points.  114 

4-5.  Examples  of  How  Boundary  Parameters  are  Specified.  128 

4-6.  Centerbody  Input  Coordinate  Examples.  129 

4-7.  Example  Demonstration  of  Nozzle  Geometry  Specification 

with  a  Generalized  Nozzle  Exit  Configuration.  130 


LIST  OF  TABLES 


Table 

3-1. 

Definition  of  FORTRAN  Symbols. 

3-2. 

Overall  Flow  of  Program  Input  ! 

Format . 

3-3. 

Input  Format. 

3-4. 

Input  Variable  Descriptions. 

3-5. 

Output  Symbol  Descriptions. 

3-6. 

Input  Data  Card  Listing  Sample 

Case. 

4-1. 

List  of  FORTRAN  Symbols. 

4-2. 

Suggested  Input  Format. 

4-3. 

Input  Variable  Definitions. 

4-4. 

Preset  Input  Values. 

4-5. 

Input  Data  Card  Listing  Sample 

Case . 

vii 


1.0 


SUMMARY 


This  supplement  to  the  Task  6,  Noise  Abatement  Nozzle  Design  Guide 
documents  two  computerized  jet  noise  prediction  techniques:  the  engineering 
correlation  method,  and  the  unified  aeroacoustic  prediction  model.  A  com¬ 
plete  description  of  the  computer  programs  are  provided,  including  examples 
of  input  preparation  and  output  cases,  plus  a  listing  of  the  FORTRAN  com¬ 
puter  code. 

1 • 1  THE  ENGINEERING  CORRELATION  (M*S)  METHOD 

A  comprehensive,  empirical,  jet-noise-prediction  method  has  been  de¬ 
veloped  by  correlating  extensive  data  from  this  program  with  available  data 
from  other  published  sources.  This  engineering  correlation  prediction  model 
has  been  designated  as  the  M*S  model  (after  the  authors:  Motsinger  and 
Sieckman)  for  ease  of  reference,  as  well  as  to  distinguish  it  from  the  more 
theoretical  prediction  model  (M*G*B)  developed  by  authors  Mani,  Gliebe  and 
Balsa. 


The  data  were  correlated  by  means  of  basic  engineering  principles  and 
physical  parameters.  The  resulting  M*S  prediction  methods  includes  unsup¬ 
pressed  conical  nozzles;  multitube  and  multichute,  single  and  dual-flow, 
suppressed  nozzles;  and  mul ti tube/chute  nozzles  with  hardwall  and  treated 
ejectors . 


1 . 2  THE  UNIFIED  AEROACOUSTIC  PREDICTION  (M*C*B)  METHOD 

A  unified  aerodynamic/acoustic  prediction  technique  has  been  developed 
for  assessing  the  noise  characteristics  of  suppressor  nozzles.  The  technique 
ucilizes  an  extension  of  Reichardt's  method  so  as  to  provide  predictions  of 
the  jet  plume  flow  field  (velocity,  temperature  and  turbulence  intensity  dis¬ 
tributions).  The  turbulent  fluctuations  produced  in  the  mixing  regions  of 
the  jet  are  assumed  to  be  the  primary  source  of  noise  generation,  as  in  the 
classical  theories  of  jet  noise.  The  altering  of  the  generated  noise  by  the 
jet  plume  itself  as  it  propagates  through  the  jet  to  the  farfield  observer 
(sound/flow  interaction  or  fluid  shielding)  is  modeled  utilizing  the  high- 
frequency  shielding  theory  based  on  Lilley's  equation. 

These  basic  modeling  elements  (flow  field  prediction,  turbulent  mixing 
noise  generation,  and  sound/flow  interaction)  have  been  coupled  together  in 
a  discrete  volume-element  formulation.  The  jet  plume  is  divided  into  elemen¬ 
tal  volumes,  each  roughly  the  size  of  a  representative  turbulence  correla¬ 
tion  volume  appropriate  to  that  particular  location  in  the  plume.  Each 
volume  element  is  assigned  its  men  characteristic  frequency,  spectrum,  and 
acoustic  intensity.  The  sound/flow  interaction  effects  for  each  volume  ele¬ 
ment  are  evaluated  from  the  flow  environment  of  the  element.  The  individual 
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volume  elements  are  assumed  to  be  uneorrelated  with  each  other,  so  that  the 
total  contribution  fa  uie  farfield  is  simply  the  sum  of  the  individual  volume 
element  contributions. 


The  programs  presented  herein  are  primarily  directed  toward  prediction 
of  high-velocity  jet  noise  (1500-2900  feet  per  second)  for  arbitrary  nozzle 
shapes,  including  sound  pressure  level  spectra  at  any  observer  location. 
Static  as  well  as  in-flight  capability  is  included  in  both  models;  however, 
tlie  1  light  data  base'  and  subsequent  verifications  are  somewhat  limited  at 
the  time  of  this  program's  conclusion. 


2.0  INTRODUCTION 


Many  jet  noise  suppressor  nozzles  have  been  designed  utilizing  intuitive 
notions  of  how  to  suppress  jet  noise  which  have  demonstrated  substantial 
noise  reduction,  but  often  at  the  expense  of  considerable  thrust  loss  as  well 
as  increased  engine  weight,  manufacturing  cost,  and  complexity.  Seemingly 
minor  changes  in  suppressor  nozzle  design,  for  the  purpose  of  improving 
chrust  performance,  often  result  in  substantial  loss  of  noise  suppression. 

It  is  therefore  highly  desirable  to  have  available  a  quantitative  prediction 
technique  for  estimating  the  aerodynamic  flow  field  and  acoustic  characteris¬ 
tics  of  suppressor-type  nozzle  configurations,  so  that  design  and  optimiza¬ 
tion  studies  can  be  made  prior  to  construction  and  testing  in  order  to  mini¬ 
mize  the  time  and  cost  of  development.  Ideally,  any  technique  should  be 
sensitive  to  the  controllable  design  variables  and  contain  a  little  empiricism 
as  possible.  When  empiricism  is  necessary,  it  should  be  based  more  or  less 
on  physical  characteristics  (flow,  acoustic  propagation,  etc.)  engineering 
principles  rather  than  on  geometric  parameters. 

The  computer  programs  included  herein  represent  a  conventional  engineer¬ 
ing  correlation  technique  and  a  more  theoretical  approach  derived  from  engi¬ 
neering  principles.  The  design  engineer  can  exercise  either  or  both  models, 
depending  on  the  type  of  results  required.  The  correlation  method  provides  a 
preliminary  design  prediction  of  aerodynamic  and  acoustic  performance;  the 
theoretical  M*G*B  method  provides  a  means  of  assessing  the  relative  impor¬ 
tance  of  various  jet  noise  mechanisms. 

Section  3.0  describes  the  computer  program  for  the  engineering  correla¬ 
tion  jet  noise  prediction  method  (M*S  model);  Section  4.0  presents  the  com¬ 
puter  program  for  the  unified  aeroacoustic  prediction  method  (M*G*B  model). 
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3.0  ENGINEERING  CORRELATION  (M*S)  JET  NOISE 
PREDICTION  COMPUTER  PROGRAM 


3.1  INTRODUCTION 


This  section  documents  the  computer  program  for  the  prediction  of  jel 
noise  by  the  engineering  correlation  method  (M*S) .  The  mathematical  model 
appears  in  detail  in  Reference  1.  A  description  of  the  computer  program  is 
provided  herein  including  examples  of  input  preparation  and  output  cases,  plus 
a  listing  of  the  FORTRAN  computer  code. 

The  computer  program  is  written  in  FORTRAN  Y  language.  It  has  been 
programmed  for  use  on  both  the  GE/Honeywell  6080  and  the  CDC  7600  computers. 

The  range  of  valid  application  of  the  program,  the  limiting  assumptions, 
and  documentation  of  the  data  base  used  for  developing  the  correlation  can  be 
found  in  both  the  Task  3  (Reference  1)  and  Task  6  (Reference  2)  reports. 


3 . 2  PROGRAM  NOMENCLATURE 


Table  3-1  defines  the  FORTRAN  symbols  used  in  the  program.  The  listing 
and  descriptions  of  input  variables  are  given  in  the  Input  Description  section. 


3.3  DESCRIPTION  OF  PROGRAM  AND  SUBROUTINES 


Table  3-2  gives  a  description  of  the  overall  flow  of  the  computer  program 
including  all  routines  used  in  each  step.  Figure  3-1  gives  a  detailed  flow 
chart  of  the  computer  program  logic.  A  description  of  the  main  program  and  each 
of  the  subroutines  is  given  in  the  following  paragraphs. 

M*S  Routine  -  This  routine  reads  the  input  curves  needed  for  the  various 
prediction  routines.  Depending  on  nozzle  type  it  reads  the  nozzle  input, 
initializes  variables,  and  computes  flow  parameters  and  flow  and  physical 
geometries.  The  computation  of  gamma  (ratio  of  specific  heats)  involves  an 
Iteration  using  input  temperature  and  pressure  ratio.  The  output  and  use  of 
prediction  subroutines  are  controlled  by  this  routine. 

Following  the  preliminary  calculations,  control  is  routed  through  the 
multielement,  conical,  or  dual-flow  section  of  the  program.  In  the  multi¬ 
element  part,  calculations  are  first  made  for  the  postmerged  noise.  The 
coefficients  for  the  Potter  and  Crocker  equation  are  set  up,  and,  because  it 
is  a  third-order  equation  (after  simplification),  a  Newtonian  convergence 
routine  is  used  to  determine  the  first  root.  Density  and  diameter  are  then 
calculated  and  a  check  is  made  for  other  possible  roots.  Static  and  total 
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Table  3-1.  Definition  of  FORTRAN  Symbols. 


FORTRAN 

Related 

Symbol 

Meaning 

Subroutines 

A 

Ejector  treatment  parameters 

MS,  EJECTS 

AA8,  A8 

Inner  nozzle  flow  area 

MS,  SHKSUB 

AJ 

Acoustic  angle,  degrees 

MS,  SUB 3,  SUB 5 

EXTP,  SHKSUB,  EJECTS 

AJA 

Jet  plume  spreading  angle,  radians 

AJR 

Acoustic  angle,  radians 

MS,  EXTP,  EJECTS 

ALT 

Input  altitude  or  arc  distance 

MS,  EXTP 

AN 

Noy  Weighting 

PNLPT 

ANl 

Number  of  elements 

MS 

ASK 

Intermediate  variable 

PNTT8 

AO 

Ambient  speed  of  sound 

MS,  SUB1 

SHKSUB,  PNTT8 

A1 

Intermediate  variable 

MS,  EJECTS 

A1 

Ratio  of  merged  to  exit  area 

MS 

A2 

Ratio  of  merged  to  exit  area 

EJECTS 

A3 

Single-flow  nozzle  total  exit  area 

MS 

A3 

Intermediate  variable 

EJECTS 

A4 

Intermediate  variable 

MS 

A4 

Ejector  treatment  PWL  Insertion  loss 

EJECTS 

A5 

Area  of  multielement  merged  stream 

MS 

A5 

Ejector  treatment  SPL  insertion  loss  at 

given  acoustic  angle 

EJECTS 

A6 

Ratio  of  ejector  inlet  area  to  nozzle 

total  area 

MS,  EJECTS 

A7 

Multielement  nozzle  area  ratio 

MS 

A9 

Outer  nozzle  flow  area 

MS 

B 

Shock  strength  parameter,  B 

SHKSUB 

B1 

Intermediate  variable 

EXTP 

B2 

Intermediate  variable 

EXTP 

B3 

Intermediate  variable 

EXTP 

B8 

Tube  or  chute/spoke  cant  angle,  radians 

MS 

B9 

Tube  or  chute/spoke  cant  angle,  degrees 

MS 

C 

Normalized  OASPL  jet  mixing  noise  curve-fit 

coefficients 

MS,  SUB1 

CJ 

Ten  dB  down  value  for  EPNL 

PNTT8 

CMAX 

Intermediate  tone  correction 

TPNLC 

Cl 

Jet  mixing  noise  OASPL  corrections 

MS,  SUB1 

C1J 

Intermediate  variable 

EXTP,  SHKSUB 

C2 

Jet  mixing  noise  relative  velocity 

exponents 

MS,  SUB1 

C3 

Inner  stream  specific  heat 

MS 

C4 

Outer  stream  specific  heat 

MS 

C9 

Local  speed  of  sound 

MS,  SHKSUB 
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Table  3-1.  Definition  of  FORTRAN  Symbols  (Continued). 


FORTRAN 

Related 

Symbol 

Meaning 

Subroutines 

D 

Intermediate  variable 

MS,  PNTT8 

DE 

Hard-wall  ejector  reference  effect  at  0y 

EJECTS 

DF.K 

Flight  Effect  at  90°  on  Shock  Cell  Noise 

SHKSUB 

DEN 

Density  correction  (pj/p0)w 

SUB1 

DIS 

Intermediate  variable 

EXTP 

DJ 

Characteristic  element  dimension 

MS 

DN 

Nozzle  outer  diameter 

MS 

DOP 

Doppler  Factor 

EXTP 

DT 

Tube  diameter 

MS 

DUM 

Intermediate  variable 

SUB1 

DO 

Shock-noise  normalization  parameter 

SHKSUB 

D1 

Reference  far-field  distance 

MS,  EXTP,  SHKSUB 

D2 

Hard-wall  ejector  reference  effect 

EJECTS 

D3 

Ejector  radius  or  diameter 

EJECTS 

DA 

Equivalent  area  diameter 

MS,  EJECTS 

D5 

Merged  flow  diameter 

MS 

D7 

Initial  time  for  EPNL 

PNTT8 

D8 

Nozzle  characteristic  dimension  for  shock 

noise 

MS,  SHKSUB 

D9 

Final  time  for  EPNL 

PNTT8 

E 

Jet  mixing  noise  spectral  distribution  at  6 

SUB1 

E 

Intermediate  Variable 

EXTP 

El 

Ejector  effect 

EJECTS 

E3 

EPNL 

PNTT8 

E9 

EGA  indicator 

MS,  EXTP,  PNTT8 

F 

Center  frequency 

MS,  EXTP,  SHKSUB 
PNTT8,  EJECTS 

F 

Intermediate  variable 

TPNLC 

FP 

Peak  frequency 

EJECTS 

FO 

Critical  frequency  for  effective  number  of 

elements 

MS 

FI 

Intermediate  variable 

MS,  SHKSUB 

F2 

Intermediate  variable 

MS,  SHKSUB 

F3 

Intermediate  variable 

SHKSUB 

G 

Shock-cell  noise  prediction  input  curve 

MS,  SHKSUB 

GJ 

Critical  refraction  angle  indicator 

MS 

G1 

Intermediate  variable 

SHKSUB 

G2 

Outer  stream  ratio  of  specific  heats,  y 

MS 

G3 

EGA  at  output  distance 

EXTP 

G8 

Intermediate  y 

MS 

G9 

Inner  stream  ratio  of  specific  heats,  y 

MS 

H 

Output  sideline  or  arc  distance 

MS,  EXTP,  PNTT8 

HI 

Intermediate  variable 

SHKSUB 

I 

Index 

MS,  SUB1 ,  SUB5, 

SUB 4,  SUB2,  SUB6, 
EXTP,  SHKSUB,  TPNLC 
PNTT8,  EJECTS 
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Table  3-1.  Definition  of  FORTRAN  Symbols  (Continued). 


FORTRAN 

Related 

Symbol 

Meaning 

Subroutines 

IDCASE 

Case  Description 

MS 

I  DENT 

Run  Description 

MS 

1M 

Intermediate  variable 

MS 

IP 

Intermediate  variable 

EJECTS 

II 

Indicator 

TPNLC 

HAS 

Noise  component  identification 

MS,  PNTT8 

I I CASE 

Case  Description 

MS,  PNTT8 

IIP 

Intermediate  variable 

MS 

ISPLF 

Intermediate  variable 

TPNLC 

J 

Index 

All  Subroutines 

JJ 

Index 

PNTT8,  EJECTS 

K 

Index 

SUBl ,  SUB 3 

KK 

Jet  mixing  noise  spectral  distribution 

curve-fit  coefficients 

MS,  SUBl 

KSTART 

Index 

SHKSUB 

KT 

Intermediate  variable 

PNTT8 

KO 

Intermediate  variable 

MS 

K1 

Extrapolation  indicator 

MS,  SUB 3 

K2 

Intermediate  variable 

MS 

K6 

Intermediate  variable 

SUBl,  EJECTS 

K7 

Shock-noise  case  indicator 

MS 

K8 

Index 

SHKSUB,  EJECTS 

K9 

Print  Indicator 

MS 

L 

PNL  calculation  coefficients 

MS,  PNLPT 

LI 

Output  acoustic  range 

EXTP 

L2 

Reflected  axial  source  location 

EJECTS 

L3 

Ejector  length 

EJECTS 

L8 

Ejector  length  effect 

EJECTS 

L9 

Ejector  length  to  suppressor  nozzle 
equivalent  diameter 

M 

Mach  number 

MS,  EJECTS 

MP 

Maximum  PNL 

PNTT8 

MM 

Intermediate  variable 

MS 

N 

Number  of  elements  in  nozzle 

MS 

NFI.T 

Flight  Effects  Exponent  Indicator 

MS,  SUBl 

N1 

Angle  indicator 

MS,  SUBl 

0 

OASPL 

SUBl,  SUB 3,  PNTT8 

OJ 

Critical  refraction  angle 

MS,  EJECTS 

09 

OAPWL 

SUB5,  SUB6,  PNTT8 

P 

PNL 

SUB 3,  PNTT8 

PA 

Air  attenuation 

EXTP 

PJ 

Intermediate  variable 

MS 

PTCOR 

Tone  correction 

TPNLC 

P0 

Ambient  static  pressure 

MS 
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Table  1- 1  .  Definition  of  FORTRAN  Symbols  (Continued). 


FORTRAN 

Related 

Symbo 1 

Meaning 

Subroutines 

PI 

it  (3.14159) 

EXTP,  SHKSUB 

P3 

Frequency 

EXTP,  EJECTS 

PA 

Inner  nozzle  total  to  ambient  pressure 

ratio 

MS 

P5 

Outer  nozzle  total  to  ambient  pressure 

rat  io 

MS 

P9 

Nozzle  total  to  ambient  pressure  ratio 

MS 

Q 

Spherical  spreading  effect 

EXTP 

Q1 

Intermediate  variable 

MS,  PNTT8 

Q2 

Jet  mixing  noise  normalization  parameter 

SUB1 

R 

Intermediate  storage  variable 

SUB4 ,  SUB 6 

RJ 

Ambient  density 

MS,  SUB1 

RJ1 

Intermediate  variable 

SUB1,  PNTT8 

RP 

Centerbody  plug  radius 

MS 

RS,  RR 

Specific  resistance 

EJECTS 

RVK 

Flight  Kf foots 

SUB1 

RX 

Specific  reactance 

EJECTS 

R 1 

Tube  equivalent  radius 

MS 

R2 

Nozzle  outer  diameter 

MS 

R) 

Inner  flow  density 

MS 

R4 

Chute/spoke  outer  flow  width 

MS 

R5 

Outer  flow  density 

MS 

R6 

Chute/spoke  inner  flow  width 

MS 

R7 

Outer  nozzle  duct  height 

MS,  SUB1 

R8 

Outer  nozzle  radius  ratio 

MS 

R9 

Centerbody  plug  radius 

MS 

S 

Predicted  SPL 

MS,  SUB1 ,  SUB 3 
SUB 5,  SUB4,  SU1 
SUB 6,  SHKSUB, 
PNTT8 

SBAR 

Intermediate  variable 

TPNLC 

SC 

Intermediate  variable 

TPNLC 

SJ 

Intermediate  variable 

MS,  PNTT8 

SL 

Input  sideline  distance 

MS,  EXTP 

SP 

Intermediate  variable 

TPNLC 

SPI. 

Intermediate  variable 

TPNLC 

SPLP 

Intermediate  variable 

TPNLC 

SPLPP 

Intermediate  variable 

TPNLC 

SS 

Outer  chute/spoke  width 

MS 

SX 

Source  location 

MS 

SI 

Shock-cell  noise  prediction  input  curves 

MS,  SHKSUB 

SLJ 

Outer  element  spacing  to  characteristic 

diameter  ratio 

MS 

S2J 

Relative  source  strength 

EJECTS 
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Table  3-1.  Definition  of  FORTRAN  Symbols  (Continued) 


FORTRAN 

Symbol  Meaning 


Related 

Subroutines 


S6 

Nozzle  outer  radius 

MS,  EJECTS 

T 

Temperature 

SUB1 

T 

PNL 

SUB  3 

T 

Flyover  time 

PNTT8 

TC 

Cutoff  effect 

MS 

TC2 

Intermediate  variable 

TPNLC 

TC3 

Intermediate  variable 

TPNLC 

TJ 

Intermediate  variable 

PNLPT,  PNTT8 

TT 

Intermediate  variable 

PNTT8 

TT3, 

T3 

Nozzle  total  temperature 

MS 

TT4 , 

T4 

Inner  nozzle  total  temperature 

MS 

TT5 , 

T5 

Outer  nozzle  total  temperature 

MS,  SUB1 

TZ 

Initial  time  for  EPNL 

PNTT8 

TO 

Ambient  temperature 

MS,  SUB1 ,  PNTT8 

T1 

Intermediate  variable 

PNTT8,  EJECTS 

T2 

Intermediate  variable 

MS 

T8 

Total  temperature 

MS,  SUB1 

U 

Arc  or  sideline  indicator 

MS,  EXTP ,  PNTT8 

U3 

Nozzle  fully  expanded  velocity 

MS 

U5 

Outer  nozzle  fully  expanded  velocity 

MS 

V 

Intermediate  variable 

SUB 3,  PNLPT 

VJ 

Suppressor  merged  velocity 

MS 

VO 

Aircraft  velocity 

MS,  SUB1, 

SHKSUB ,  PNTT8 

VI 

Ratio  of  merged  velocity  to  exit  velocity 

MS 

V6 

Intermediate  variable 

MS 

V7 

Intermediate  variable 

MS 

V8 

Fully  expanded  jet  velocity  input  to  jet 

mixing  noise  routine 

MS,  SUB1 

V9 

Fully  expanded  jet  velocity  input  to 

shock-cell  noise  routine 

MS,  SHKSUB 

W 

Density  exponent  curve-fit  coefficients 

MS,  SUB1 

WE 

Density  exponent 

SUB1 

WJ 

Intermediate  variable 

SUB1,  PNTT8 

W4 

Inner  stream  weight  flow 

MS 

W5 

Outer  stream  weight  flow 

MS 

W8 

Weight  flow 

MS,  SUB1 

X 

Source  location 

MS,  EJECTS 

X 

SPL 

SUB 3,  EXTP,  PNLPT 

XJ 

Intermediate  variable 

SUB1 ,  EJECTS 

XM 

Point  of  merging 

MS 

XO 

Potter  and  Crocker  equation  coefficient 

MS 

XI 

Potter  and  Crocker  equation  coefficient 

MS 

X2 

Potter  and  Crocker  equation  coefficient 

MS 
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Table  3-1.  Definition  of  FORTRAN  Symbols  (Concluded) 


FORTRAN 

Related 

Symbol 

Meaning 

Subroutines 

X3 

Potter  and  Crocker  equation  coefficient 

MS 

X4 

Specific  reactance 

EJECTS 

Y 

PWL 

SUB 5,  SUB4,  SUB6 

PNTT8 

YJ 

Intermediate  variable 

SUB 5,  EJECTS 

Y1 

Intermediate  variable 

MS,  SUB 4,  SUB 6 

Y1J 

Intermediate  variable 

MS 

Y2 

Intermediate  variable 

MS 

Y9 

Nozzle  type  indicator 

MS,  SUB1 

ZI 

Intermediate  variable 

SHKSUB 

ZJ 

Intermediate  variable 

EXTP,  EJECTS 

ZK 

Intermediate  variable 

SHKSUB 

ZZ 

Effective  number  of  elements  effect 

MS 

Zl 

Intermediate  variable 

SUB1,  PNTT8 

Z2 

Intermediate  variable 

MS 

Z3 

Intermediate  variable 

MS,  PNTT8 

Z5 

Number  of  rows  of  tubes 

MS 

Z8 

Effective  number  of  elements  adder 

MS 

Z9 

Total  number  of  elements  adder 

MS 

Z9 

Constant 

MS,  SUB 2 
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Table  3-2.  Overall  Flow  of  Program. 

1 .  Read  Input  Curves  (M*Sl . 

2.  Read  Input  and  Calculate  Flow  Parameters  for  each  Stream  (M*S). 

The  Following  through  11  are  used  or  Skipped  as  Necessary. 

3.  Determine  Postmerged  Noise  (M*S,  SUB1 ,  SUBS). 

4.  Determine  Premerged  Noise  (M*S,  SUB1). 

5.  Determine  Premerged  Cutoff  and  Shielding  Effects  (M*S). 

6.  Calculate  Ejector  Effects  and  Correct  the  Premerged  Noise  (M*S,  EJECTS, 
SUB5). 

7.  Sum  the  Premerged  and  Postmerged  Noise  (SUB6). 

8.  Calculate  Shock  Noise  for  Outer  Stream  and  Apply  Cutoff,  Shielding, 
and  Ejector  Effects  (M*S,  SHKSUB,  EJECTS,  SUB5). 

9.  Add  to  the  Sum  of  Premerged  and  Postmerged  (SUB6). 

10.  Calculate  Shock  Noise  for  Inner  Stream  (M*S,  SHKSUB,  SUB5) . 

11.  Add  to  the  Sum  of  Premerged  and  Postmerged  and  Outer  Stream  Shock 
(SU B6). 

12.  Extrapolate  and  Calculate  OASPL,  PNL  and  PNLT  (this  may  be  done  after 
each  Component  is  Calculated  for  Print  Purposes)  (SUB3). 

13.  Print  Output  and  Calculate  EPNL  (PNTT8) . 
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Figure  3-1.  Computer  Program  Flow  Chart  (Continued), 
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Figure  3-1*  Computer  Program  Flow  Chart  (Continued). 


Figure  3-1,  Computer  Program  Flow  Chart  (Continued). 
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Figure  3-1.  Computer  Program  Flow  Chart  (Continued). 


19 


19 


20 


Figure  3-1.  Computer  Program  Flow  Chart  (Continued). 
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Figure  3-1.  Computer  Program  Flow  Chart  (Continued) 
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Figure  3-1.  Computer  Program  Flow  Chart  (Continued). 
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b)  SUBROUTINE  SUBS 


Figure  3-1.  Computer  Program  Flowchart  (Continued) 


i  )  SUBROUTINE  SUIJS 


Figure  3-1.  Computer  Program 

Flowpath  (Continued). 
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d)  SUBROUTINE  SUB4 


e)  SUBROUTINE  SUB2 


/  DO  FOr\ 
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/ RETURN \ 
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/  RETURN  \ 


Figure  3-1.  Computer  Program 

Flowchart  (Continued). 
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i )  SUBROUTINE  SUB6 


Figure  3-1.  Computer  Program 

Flowchart  (Continued) 
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Figure  3-1.  Computer  Program 

Flowchart  (Continued) 
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i)  SUBROUTINE  EJECTS 


Figure  3-1.  Computer  Program  Flowchart  (Continued). 
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Figure  3-1.  Computer  Program  Flowchart  (Continued). 
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NO  FLOW  CHART  IS  SUPPLIED  FDR  THE  FOLLOWING 
ROUTINES  BECAUSE  OF  THE  COMPLETE  NATURE  OF 
THEIR  DOCUMENTATION  IN  PUBLISHED  LITERATURE 

•  SHKSUB 

•  PNLPT 

•  TPNLC 


Figure  3-1 .  Computer  Program  Flow 
Chart  (Concluded). 
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temperature  are  determined,  the  input  variables  to  the  conical  nozzle  noise 
routine  are  set,  the  noise  is  calculated,  and  flight  effects  are  applied  if 
necessary.  This  component  is  then  extrapolated  and  (if  desired)  printed. 

The  pretnerged  noise  is  then  calculated.  The  effective  number  of  tubes 
and  the  critical  angle  are  determined.  Then  the  length  of  the  potential 
core,  Xc,  the  point  of  merging  (used  for  cutoff  only),  and  the  radius  ratio 
are  determined.  The  axial  location  of  the  beginning  of  peak  noise  generation, 
Xp,  and  the  critical  frequency  for  absorption  are  calculated  before  entering 
the  frequency  loop  to  calculate  source  locations,  absorption  effects,  and 
cutoff  effects.  These  are  then  applied  to  all  angles  forward  of  critical 
with  angles  aft  of  critical  set  equal  to  critical  angle  SPL.  Ejector  effects 
are  determined  and  applied  before  extrapolation  and  (if  desired)  printed. 
Shock-cell  noise  (if  applicable)  is  determined  after  summing  the  premerged 
and  postmerged  components.  It  is  then  corrected  for  ejector  effects  and 
flight  effects,  whereupon  multielement  corrections  are  applied,  extrapolated, 
printed,  and  added  to  the  other  components.  The  total  is  then  extrapolated 
(if  required)  and  printed,  and  a  return  is  made  for  the  next  case. 

The  conical  part  of  the  routine  calculates  the  conical  mixing  noise  and 
shock  noise,  applies  flight  effects,  extrapolates  and  prints  them  separately 
if  desired,  sums  them,  and  prints  the  total;  after  which,  a  return  is  made 
for  the  next  case. 

The  coannular  part  uses  the  premerged  and  postmerged  routines  of  the 
multielement  part  if  a  suppressor  is  involved.  Variables  are  set,  and,  if  a 
suppressor  is  involved,  the  postmerged  routine  of  the  multielement  part  is 
entered  to  calculate  merged  flow  conditions.  Mixed  conditions  are  then 
determined  and  the  merged  noise  is  calculated,  extrapolated,  and  printed  (if 
desired) . 

The  premerged  noise  is  now  calculated  in  accordance  with  whether  a  sup¬ 
pressor  is  present  or  not.  This  component  is  extrapolated,  printed  if  desired, 
and  added  to  the  postmerged.  Outer-stream,  shock-cell  noise  is  determined, 
depending  on  whether  a  suppressor  is  present  or  not,  extrapolated,  printed 
(if  desired),  and  added  to  the  other  components.  Finally,  the  inner  stream 
shock  is  computed,  extrapolated,  printed  (if  desired),  and  added  to  the  other 
components.  The  total  is  then  extrapolated  as  required,  and  printed;  and 
control  is  returned  for  the  next  case. 

SUB1  Subroutine  -  This  subroutine  provides  SAE  ART  876  (1975  revision) 
conical  nozzle  noise  predictions  and  determines  and  applies  mixing  noise  flight 
effects.  Use  and  limitations  are  as  described  in  the  aforementioned  documents. 
Output  from  this  routine  is  on  a  one-foot  arc.  Basically,  polynominal  curve 
fits  of  the  data  in  SAE  ARP  876  (1975  revision)  were  used.  A  correction  was 
made  to  the  predicted  OASPL  to  increase  the  accuracy  of  the  routine  based  on 
available  data  on  suppressor  nozzles.  This  correction  amounts  to  +1  dB  at 
all  angles  and  frequencies. 
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SUB3  Subroutine  -  This  routine  resets  the  variables  for  input  into  the 
extrapolation  and  PNL  calculation  subroutines.  It  determines  whether  extra¬ 
polation  is  required  and  calls  LXTP.  PNLPT  is  called  to  determine  PNL  and 
OASPL.  TPNLC  is  called  from  PNLPT  to  determine  PNLT.  The  variables  are  then 
reset  maintaining  the  newly  calculated  values. 

SUB5  Subroutine  -  This  routine  calculates  sound  power  level  from  sound 
pressure  level  for  each  one-third-octave  band,  and  then  antilogarithmically 
sums  them  to  obtain  the  overall  levels. 

SUB4  Subroutine  -  This  routine  places  previously  calculated  sound  pres¬ 
sure  level  and  sound  power  level  in  other  variable  name  storage  for  future 
use  in  the  program. 

SUB2  Subroutine  -  This  routine  adds  a  constant  value  to  the  one-third- 
octave  band  SPL  at  all  angles  and  frequencies. 

SUB6  Subroutine  -  This  routine  antilogarithmically  sums  different  SPL 
and  PWL  spectra  to  obtain  a  total  spectra,  and  then  sums  the  total  PWL 
spectrum  to  obtain  OAPWL. 

EXTP  Subroutine  -  This  routine  extrapolates  an  input  spectrum  to  a 
desired  acoustic  range  using  the  inverse-square  law  (spherical  spreading), 
air  attenuation  per  SAE  ARP  866  (Reference  3),  and,  if  desired,  extra  ground 
attenuation  (EGA)  per  the  routine  presented  in  SAE  AIR  923  (Reference  4).  A 
curve  fit  of  the  59°  F,  70%  relative  humidity,  standard-day-air  attenuation 
is  used,  as  well  as  curve  fits  for  EGA.  The  routine  automatically  accounts 
for  range  changes  from  angle  to  angle  on  a  sideline  and  includes  the  option 
of  a  100-ft  layer  of  EGA,  full  EGA,  or  no  EGA  as  per  SAE  AIR  923. 

SHKSUB  Subroutine  -  This  routine  predicts  shock-cell  noise  by  the  pro¬ 
cedure  defined  in  SAE  ARP  876  (1976  proposed  revision).  Output  from  this 
routine  is  on  a  one-foot  arc.  The  definition  of  D8  was  varied  to  allow 
calculations  for  nonround  nozzles.  Shock-cell  noise  flight  effects  are 
determined  and  applied  in  this  section. 

PNLPT  Subroutine  -  This  routine  sums  the  SPL  in  a  given  sp .ctrum  anti¬ 
logarithmically  to  obtain  OASPL  and  uses  the  procedure  defined  in  SAE  ARP 
865A  (Reference  31)  to  calculate  PNL. 

TPNLC  Subroutine  -  This  routine  calculated  tone-corrected  PNL  via 
Section  B36.3  of  the  FAA  Noise  Certification  Document  (Nov.  17,  1969)  as  a 
function  of  the  uncorrected  one- third-octave  spectrum  SPL. 

PNTT8  Subroutine  -  This  routine  sets  the  format  and  prints  the  noise 
output  from  the  main  program.  It  prints  the  identification  of  the  noise  out¬ 
put  and  one-third-octave  band  SPL  and  PWL  for  24  frequencies  and  15  angles 
(20°  to  160°  to  the  inlet)  as  well  as  OASPL,  PNL,  and  PNLT  for  each  angle. 

The  second  part  of  the  routine  calculates  EPNL  (if  required)  according 
to  the  procedure  described  in  EAR  Part  36,  using  PNL  rather  than  PNLT.  Times 
associated  with  given  icoust ic  angles  for  a  level  flyover  (assuming  the 
engine  centerline  is  parallel  to  the  ground)  are  determined  first.  Peak  PNL, 


the  associated  time,  and  the  10-dB  down  levels  are  determined.  Initial  and 
final  times  are  then  determined  by  linear  interpolation  (using,  when  necessary, 
extrapolation  using  the  first  or  last  two  points).  The  PNL  history  is  then 
integrated  between  the  10-dB  down  points  by  summing  half-second  increments 
(determined  by  linear  interpolation)  to  obtain  the  duration  correction.  This 
is  added  to  the  maximum  PNL  to  obtain  EPNL;  the  EPNL  is  then  printed.  It 
should  also  be  noted  that  the  program  automatically  calculates  an  EPNL  for 
static  sideline  cases  assuming  a  300  ft/sec  flyover  velocity. 

EJECTS  Subroutine  -  This  routine  first  determines  the  effect  of  a  hard- 
wall  ejector  of  given  geometry  in  terms  of  the  reference  SPL.  Directivity 
and  spectral  effects  are  then  determined.  If  no  treatment  is  present  in  the 
ejector,  control  is  returned  to  the  main  program.  If  treatment  is  present, 
an  impedance  prediction  routine  for  SDOF  treatment  (single  degree  of  freedom) 
is  entered.  The  resistance  and  reactance  of  the  treatment  panel  is  determ¬ 
ined;  this  yields  a  coefficient  of  absorption.  The  location  of  a  given 
source  and  the  strength  relative  to  the  peak  are  then  calculated.  The  co¬ 
efficient  of  absorption  multiplied  by  the  number  of  reflections  for  a  given 
acoustic  angle  plus  the  relative  source  strength  when  summed  over  all  sources 
yields  an  SPL  reduction.  This,  when  integrated  over  all  angles,  gives  a  sound 
power  insertion  loss.  This  reduction  is  log-averaged  over  the  lower  limiting, 
center,  and  upper  limiting  frequencies  for  the  given  one-third-octave  band. 

The  sound  power  insertion  loss  is  then  converted  into  a  delta  SPL  for  each 
acoustic  angle  and  added  to  the  hard-wall  effect.  Control  is  then  returned 
to  the  main  program. 

3 . 4  INPUT  DESCRIPTION 

The  input  data  are  supplied  through  NAMELIST  input  format.  Any  number  of 
successive  cases  can  be  run  consecutively,  limited  only  by  the  users  execu¬ 
tion  time  available.  Each  successive  case  requires  only  the  INPUT  NAMELIST. 

The  data  from  preceding  cases  remain  in  storage;  thus,  only  those  variables 
which  are  to  be  changed  from  the  preceding  case  input  value  need  be  included 
in  the  INPUT  file  of  succeeding  cases. 

The  input  format  is  given  in  Table  3-3.  The  definitions  of  each  of  the 
input  variables  given  in  Table  A-3  are  given  in  Table  3-4.  All  variables  are 
preset  to  zero  before  the  first-case  input  is  read.  Only  the  input  variables 
listed  under  a  nozzle  type  in  Table  3-3  need  be  input  for  any  case.  Notes  on 
the  input  follow  the  tables.  Further  descriptions  of  input  variables  are 
given  in  Figures  3-2  and  3-3. 

3.4.1  Notes  on  Input 

1.  The  ALT  variable  is  used  as  the  main  distance  indicator;  therefore,  for 
ground  static  arc  or  sideline  cases  the  distance  of  interest  is  input  through 
this  variable,  and  the  SL  variable  is  set  to  zero.  In  flyover  cases,  ALT  is 
used  as  the  altitude  indicator,  and  SL  is  used  as  the  sideline  distance. 

2.  EGA  is  "Extra  Ground  Attenuation"  as  defined  in  SAE  AIR  923  "Method  for 
Calculating  the  Attenuation  of  Aircraft  Ground  to  Ground  Noise  Propagation 
During  Takeoff  and  Landing."  The  "100-ft  layer"  is  defined  in  Figure  3  of  the 
above-mentioned  document. 


3.  Major  nozzle  dimensions  are  input  in  feet;  element  or  ejector-treatment 
dimensions  are  input  in  inches.  This  alleviates  inputting  very  small  numbers 
(i.e.,  0.1  inch  versus  0.0083  foot). 

4.  Cant  angles  for  multitube  and  multichute/spoke  nozzles  are  defined  in 
Figure  3-4. 

5.  The  "A"  variables  are  input  as  10  if  treatment  other  than  SDOF  is  used. 

In  this  case  RR  and  RX  must  be  input. 

6.  The  specific  resistances  and  reactances  of  the  treatment  used  in  the  ejec¬ 
tor  are  input  through  the  RR  and  RX  variables.  Values  at  the  lower  limiting, 
upper  limiting,  and  midpoint  frequencies  are  used.  For  ease  of  input,  the 
program  assumes  the  value  at  the  upper  limiting  frequency  of  one  one-third- 
octave  band  to  be  equal  to  the  value  at  the  lower  limiting  frequency  of  the 
next  highest  band.  Therefore,  only  49  values  must  be  input. 

3.5  OUTPUT  DESCRIPTION 

The  output  format  is  generally  self-explanatory.  The  input  is  printed 
out  using  the  nomenclature  defined  in  Table  3-5.  Output  flow  conditions 
follow.  Finally,  SPL  and  PWL  spectra,  OASPL,  OAPWL,  PNL,  PNLT,  and  EPNL  are 
printed  as  required. 

A  warning  flag  is  built  into  the  iterations  for  gamma  and  merged  vel¬ 
ocity.  The  flag  message  for  either  iteration  is:  DTD  NOT  CONVERGE;  and 
when  it  appears  the  run  terminates.  Cross  input  errors  have  been  the  only 
cause  of  this  message  encountered  in  the  development  of  the  program. 

At  the  beginning  of  each  run,  an  unlimited  number  of  cards  can  be  input 
for  the  run  identification.  (A  case  identification  card  is  available  before 
each  case  also).  The  format  for  each  card  is: 

60  -  Character  Title  Card,  Columns  1-60 

To  enter  the  case  section  of  the  input  the  following  card  is  required: 
CASES  (Starting  in  Column  2) 

The  run  or  case  identification  cards  may  be  omitted  but  the  "CASES" 
card  must  be  present.  The  case  identification  is  saved  and  will  be  printed 
on  succeeding  cases  unless  another  case  identification  card  is  read. 


v  r 


Table  3-3.  Input  Format. 


(FOR  CONICAL  NOZZLES) 


Column 

2 

♦ 


(60-Character 

Identification 

Card, 

Columns  1-60) 

$  INPUT  Y9  = 

1, 

P9  = 

TT3  = 

9 

A9  * 

9 

K9  = 

ALT  = 

9 

SL  = 

,  * 

u  = 

E9  = 

9 

VO  = 

,  NFLT  = 

$ 
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Tabic  3-3.  Input  Format  (Continued). 
(FOR  SINGLE- FLOW,  MULTITUBE  NOZZLES) 


Column 

2 

♦ 

(60-Character  Identification  Card,  Columns  1-60) 


$  INPUT 

V  9  =  2, 

N  = 

,  RP  = 

,  B9  = 

DT  = 

,  A7  = 

,  25  «  , 

S1J  - 

,  TT3  = 

,  P9  = 

K9  = 

,  ALT  = 

,  SL  = 

U  = 

,  E9  = 

,  VO  = 

A6  = 

,  L9  = 

,  NFLT  = 

Table  3- 


.  Input  Format  (Continued). 


Column 

2 


Table  3-3.  Input  Format  (Continued). 
(FOR  SINGLE-FLOW,  MULTICHUTE/SPOKE  NOZZLES) 

Column 

2 

♦ 

(60-Character  Identification  Card,  Columns  1-60) 


$  INPUT 

Y9  =  3, 

N 

,  RP  = 

,  B9  = 

y 

RA  = 

,  R6  = 

,  SS  = 

> 

II 

y 

TT3  = 

,  P9  = 

,  NFLT  = 

K9  = 

,  ALT  = 

,  SL  = 

.  u  = 

y 

E9  = 

,  VO  = 

,  A6  = 

.  L9  = 

y 

A  =  _ 

y  y 

♦ 

y 

RR  and  RX  as  per  the  multitube  nozzle  case. 


Table  3-3.  Input  Format  (Continued). 


(FOR  DUAL-FLOW  NOZZLES  WITH  MULTITUBE 
SUPPRESSORS  ON  THE  OUTER  STREAM) 


Column 

2 

♦ 

(60-Character  Identification  Card,  Columns  1-60) 


$  INPUT 

RP  = 

Y9  =  5, 

,  DN  = 

,  AA8  * 

> 

\o 

n 

..  -  -  9 

TT4  = 

■O 

II 

,  TT5  = 

,  P5  = 

9 

N 

,  DT  = 

> 

w 

,  B9  = 

9 

Z5  = 

,  S1J  = 

,  NFLT  = 

9 

K9  = 

,  ALT  = 

,  SL  = 

,  u  = 

9 

E9  = 

,  VO  = 

,  A6  = 

,  L9  = 

9 

A 

*  9 

9 

_ _» 

RR  and  RX  as  per  multitube  case. 

$ 
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Table  3-3.  Input  Format  (Concluded) 


(FOR  DUAL-FLOW  NOZZLES  WITH  MULTICHUTE/SPOKE 
SUPPRESSOR  ON  THE  OUTER  STREAM) 


Column 

2 

♦ 

(60-Character  Identification  Card,  Columns  1-60) 


$  INPUT  Y9  = 

:  6, 

RP  = 

,  DN  = 

AA8  =  _ 

,  A9  = 

9 

TT4  = 

,  PA  = 

TT5  =  _ 

II 

9 

N  = 

,  B9  = 

,  NFLT  = 

9 

II 

,  R6  = 

SS  =  _ 

,  A  7  = 

9 

K9  = 

,  ALT  = 

SL  =  _ 

,  u  = 

9 

E9  = 

,  VO  = 

_ ,  A6  =  _ 

,  L9  = 

9 

A 

»  9 

9 

_ * 

RR  and  RX  as 

per  multitube 

case. 

$ 
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Table  3-4.  Input  Variable  Descriptions. 


(FOR  CONICAL  NOZZLES) 


Variable  Note 

P9 
TT3 
A9 
K9 


ALT 


1 


SL 


1 


U 


E9 


2 


VO 

NFLT 


Description 

Nozzle  Total  to  Ambient  Pressure  Ratio 

Nozzle  Exit  Total  Temperature,  °  R 

Nozzle  Exit  Flow  Area,  ft^ 

Print  Indicator: 

0  =  Total  Nozzle  Noise  Only 
1  =  Nozzle  Component  and  Total  Noise 

Altitude,  Ground  Sideline,  or  Arc  Dis¬ 
tance  at  which  Prediction  is  to  be  made, 
ft 

Sideline  Distance  at  Which  Prediction  is 
to  be  made,  ft  (Used  for  Flyover  Cases 
Only) 

Arc  or  Sideline  Indicator 

1  =  Predictions  to  be  made  on  an  Arc 

2  =  Predictions  to  be  made  on  a  Sideline 

(or  Flyover) 

EGA  Indicator 
0  =  No  EGA 

1  =  Full  EGA 

2  =  100-ft  Layer  of  EGA 

Aircraft  Flight  Velocity 

Flight  Effects  Indicator 

1  =  "Free  Jet" 

2  =  "True  Flight" 


-15 


./ 


Table  3-4.  Input  Variable  Descriptions  (Continued). 


(FOR  SINGLE-FLOW,  MULTITUBE  NOZZLES) 


Variable 

Note 

Description 

N 

Number  of  Tubes 

RP 

3 

Centerbody  Plug  Radius,  ft 

B9 

4 

Tube  Centerline  Cant  Angle,  degrees 

DT 

3 

Tube  Diameter,  in. 

A7 

Nozzle  Area  Ratio 

Z5 

Number  of  Rows  of  Tubes  Counting  Center 
Tube  (if  Present)  as  zero 

S1J 

Tube  Centerline  Spacing  to  Tube  Diameter 
Ratio 

TT3 ,  P9, 
SL,  U,  E9 

K9,  ALT, 

,  VO 

Same  as  Conical  Nozzle 

A6 

Ratio  of  Ejector  Inlet  Area  to  Nozzle 
Total  (or  Annulus)  Area  (Input  Zero  for 
no  Ejector) 

L9 

Ratio  of  Ejector  Length  to  Suppressor 
Nozzle  Equivalent  Diameter 

A  ( 1) 

3,5 

Ejector  Treatment  Faceplate  Thickness,  ii 

A(2) 

3,5 

Ejector  Treatment  Hole  Diameter,  in. 

A(3) 

3,5 

Ejector  Treatment  Cavity  Depth,  in. 

A(4) 

3,5 

Ejector  Treatment  Open  Area  Ratio 

RR 

6 

Ejector  Treatment  Specific  Resistance, 
Ray Is  (49  Values  Required) 

RX 

6 

Ejector  Treatment  Specific  Reactance, 
Rayls  (49  Values  Required) 
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Table  3-4.  Input  Variable  Descriptions  (Continued). 
(FOR  SINGLE-FLOW,  MULTICHUTE/SPOKE  NOZZLES) 


Variable  Note 

Description 

N 

Number  of  Elements 

RP 

3 

Centerbody  Plug  Radius,  ft 

B9 

4 

Chute/Spoke  Exit  Cant  Angle,  degrees 

R4 

Outer  Circumferential  Flow  Dimension,  in 

R6 

Inner  Circumferential  Flow  Dimension,  in 

SS 

Outer  Circumferential  Element  Dimension, 
in. 

A7 

Nozzle  Area  Ratio 

TT3 

SL, 

,  P9,  K9 ,  ALT, 

V,  E9 ,  VO, 

Same  as  Conical  Nozzle 

A6, 

L9,  A,  RR,  RX 

Same  as  Multitube  Nozzle 

47 


Table  3-4.  Input  Variable  Descriptions  (Continued). 


(FOR  DUAL-FLOW  NOZZLES  WITH  A  MULTITUBE 
SUPPRESSOR  ON  THE  OUTER  STREAM) 

Variable  Note 

RP 
DN 
AA8 
A9 
TT4 
P4 

TT5 
P5 

N,  DT,  A7,  B9, 

’  Z5 ,  S1J,  A6 ,  L9 , 

A,  RR,  RX 

K9 ,  ALT,  SL,  U,  E9,  VO 

i 

> 
s 

> 
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Description 

Centerbody  Plug  Radius,  ft 

Nozzle  Outer  Diameter,  ft 

Inner  Nozzle  Flow  Area,  ft^ 

o 

Outer  Nozzle  Flow  Area,  ft 

Inner  Nozzle  Exit  Total  Temperature,  °  R 

Inner  Nozzle  Total  to  Ambient  Pressure 
Ratio 

Outer  Nozzle  Exit  Total  Temperature,  °  R 

Outer  Nozzle  Total  to  Ambient  Pressure 
Ratio 

Same  as  Multitube  Nozzle 

Same  as  Conical  Nozzle 


Table  3-4.  Input  Variable  Descriptions  (Concluded). 


(FOR  DUAL-FLOW  NOZZLES  WITH  MULT I CHUTE /SPOKE 
SUPPRESSORS  ON  THE  OUTER  STREAM) 


Variable  Note 

RP,  DN,  AA8,  A9, 

TT4,  P4 ,  TT5,  P5 

N,  B9,  R4,  R6,  SS,  A7 

K9 ,  ALT,  SL,  U,  E9,  VO 

A6 ,  L9 ,  A,  RR,  RX 


Description 

Same  as  Dual-Flow /Multitube 

Same  as  Multichute/Spoke 
Same  as  Conical 
Same  as  Multitube 
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Conical 


i _ I _  jjmmx 


Multitube 


Figure  3-2.  Nozzle  Types  Included  in  the  Correlation. 


5i 


a 


Jet  Axis 


input  as  "ALT." 


3-3.  FORTRAN  Symbol  Convention  for  Acoustic  Arena 
Variables . 


1 


Multitube  Nozzles 


Multichute/Spoke  Nozzles 


Figure  3-4.  Definition  of  Cant  Angles  for  Multielement  Nozzles. 
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Table  3-5.  Output  Symbol  Descriptions. 


Symbol 

Description 

ARD 

Suppressor  Nozzle  Area  Ratio 

AT 

Area  of  an  Individual  Flow  Element 

A5 

Merged  Flow  Area 

A6 

Mixed  Flow  Area 

A8 

Inner  Nozzle  Flow  Area 

A28 

Outer  Nozzle  Flow  Area 

DUCT  H 

Outer  Nozzle  Duct  Height 

D5 

Diameter  of  the  Merged  Flow  Stream 

PO 

Ambient  Pressure 

PT8/P0 

Inner  Nozzle  Pressure  Ratio 

PT28/P0 

Outer  Nozzle  Pressure  Ratio 

RH05 

Density  of  the  Merged  Stream 

RH08 

Density  of  the  Inner  Stream 

RH028 

Density  of  the  Outer  Stream 

TO 

Ambient  Temperature 

TT5 

Total  Temperature  of  the  Merged  Stream 

TT6 

Total  Temperature  of  the  Mixed  Stream 

TT8 

Total  Temperature  of  the  Inner  Stream 

TT28 

Total  Temperature  of  the  Outer  Stream 

U5 

Fully  Expanded  Merged  Velocity 

U6 

Fully  Expanded  Mixed  Velocity 

U8 

Fully  Expanded  Inner  Stream  Velocity 

U28 

Fully  Expanded  Outer  Stream  Velocity 

W6 

Mixed  Stream  Weight  Flow 

PWL 

-13 

Sound  Power  Level,  dB  re:  10  watts 

OASPL 

Overall  Sound  Pressure  Level  re: 

2  dynes /m^ 

OAPWL 

Overall  PWL 

PNL 

Perceived  Noise  Level,  PNdB 

PNLT 

Tone-Corrected  PNL,  PNdB 

EPNL 

Static  Effective  Perceived  Noise  Level 
EPNdB 
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3.6  SAMPLE  CASES 


Example  cases  for  a  conical  nozzle  with  and  without  EGA,  a  dual-flow 
nozzle  with  a  multitube  suppressor  and  a  treated  ejector,  and  a  dual-flow 
nozzle  with  a  multichute  suppressor  are  given.  The  input  data  cards  are 
listed  in  Table  3-6  as  per  the  format  given  in  Table  3-3. 


Table  3-6.  Input  Data  Card  Listing  Sample  Case. 


AR  SItCKMAN  TASK  3  HIGH  VELOCITY  JET  NOISE  PRCGRaM 

GENERAL  ELECTRIC  CO.  RLOG  100  BIN  79  M.D.  H77  X22h  1 
MS  —  ENGINEERING  CORRELATION  MQOEl  .-  cnc  VFRSICN 
CASES 

CONICAL  NOZZLE  CHECK  CASE 
SINPIIT  Y9A1, 

P9A3.247,  TT3A1380,  A9A2.346,  RP*0,  K9#l, 

ALT*2U00,  UA2,  F9A0.  V0A350,  A6*0,  l 9*o,  AA4*0, 

S 

JINPUT  E9 A2$ 

DUAL  ELCW  MULTI-TU8F  CHECK  CASE 
SINPIJT  V9A5, 

RPA1.423,  PNA6.  687  ,  AA8A7.649,  A9A5'.O03,  TT4A1010, 

P4A1.567,  TT5A1632,  PSA3.2 78,  K9#|,  NA*9, 

DTA1.672,  A7A2.7S,  R9A0,  ZSA 1, 

SI J *2. 818#  ALTA320#  UA1.  F9*0,  V0*0. 

A6A0,  L 9 A 0  »  AA4*0, 

A6A1 .103,L9A3.9S2, AAa«10, 

RRAO9*0 ,31 1 , 

RXA-B7. 1  35, -77. 549,. 69.239, -61 . 153, -54'.  949, -«8.463, 
-«3.269,-3B.7fe7,-3fl’f6l  1  ,-31  .  008  , -27).  683 , -?4  .  ?  i  9 , -2 1 . 6?o , 

-19. 367, -J7. 287# -15. 484, -13.8 19, -12. 277 #-10.954 #-9.652, -8. 600, 
-7. 702, -6. 864, -6, 088, -5. 370 #-4.762, -4. ?32,-3.771, -3. 34?, 
-2.968,-2.619,-2.251  ,-1.970,-1  .722,-1.(|R7,-1  .278,-1  .077, -.882, 
-.704, -.515, -.34  7,-1  185, -.010,.  185, '.401  ,. 703,  1  .  1  ,  1.794, 4. 09  7, 

$ 

DUAL  FLOW  MULTI-CHUTE  CHECK  CASE 
fINPUT  Y9 A  6 , 

RPA.624,  DNA2.671,  AA0A.8I1,  A9A1.5B5,  TT4A1470, 

P4A1.490,  TT5A1750,  P5A3.97,  K9A0,  NA20, 

B9A0,  R4A2.874,  R6A2.060,  SSA2.155, 

A 7 A i  . 75 ,  ALTA2400,  UA2,  E9#0,  V0A35O, 

A6 A  0 ,  L9A0,  A  A  4  *  0 , 

S 


NOTE:  The  symbol  A  indicates  an  equal  sign  (=). 
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XX  XX  XXXXXX 

XXX  XXX  XXX.XXX . 

X  X  X  x  XXXX  XX 

XX  XXXX  XX  XX 

XX  XX  xx  XXXXXX 

XX  XX  XXXXXX 

XX  XX  XX 

XX  XX  XX. 

XX  AX  XXXXXX 

XX  XX  XXXXXX 


X  X  X  X  X  x  X  X X  X  X  X  XXXXXX  XXxXXX  A  X 

X  .  X  x  X  X  X  XX  XX 

XXX  XX  XXX  X  X  X  XX  X 

_  XXX  X  XXX  X.  __XXXX  __X  X 

A  A  A  X  X  >  X 

X  X  X  A  X  X  X  X  X  x  a  X  XXXXXX  A  X 


H  j(iH  v  F  IOC  I  TV  JET  NOjSF  PKOGPAN  (CONTRACT  OOT-OS.  7Q034) 

m-.i  --  kngjnfehtho  copwhation 


A^  'ItCKMAK  TA^k  1  H I  AH  VELOCITY  .  IE  7  NOISF  PROGRAM 

OF  A1'  PAL  FLt'CTPfC  CO.  PLOG  ;i  0  L-  PIN  79  m.O.  077  X??6) 

PS  —  ENGINEERING  COPRFlaTIUN  mqqel  —  cue  VERSION 


. 


I 


i 


wri  oCTTf  Jr  t  l<. :  I  St  PROGRAM  -  l  fiC.  iNf-  f  >.  Ii-.r,  CRRRFt  A  T  I  ON 


i  c  n  f-  c  rv| 
la  ^  ip  1^1 


a  ^^-jNrlir^o  -  rv^ri 

•  •  •  •  •  *  •  •  j  • 

in  a.  a  <r  •—  xir-  >o  <r  <r  <*i<v 

U“-  iP  >C  X  vOj  s£  n  £h  <C  X  XI  X 


c  o  a  ^  i/l  ^  X) 

*  u;  i  3  ^ 


POO  uj 

o  o  a. 

-J  U.  Ol 
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3.7  PROGRAM  SOURCE  CODE  LISTING 


This  section  contains  the  FORTRAN  IV  source  code  listing  for  the  engi¬ 
neering  correlation  computer  program,  suitable  for  running  on  the  CDC  7600 
computer.  The  listing  of  subroutines  is  as  follows: 

(1)  Main  Program  (MS) 

(2)  SUB1  (Contains  SUB1  through  SUB6) 

(3)  EXTP 

(4)  SHKSUB 

(5)  PNLPT 

(6)  TPNLC 

(7)  PNTT8 

(8)  A  block  data  listing 

(9)  EJECTS 
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4.0  UNIFIED  AEROACOUSTIC  PREDICTION  MODEL  (M*G*B)  COMPUTER  PROGRAM 


4 . 1  INTRODUCTION 


This  section  describes  the  computational  algorithims  and  associated 
computer  program  that  provide  the  necessary  link  between  the  symbolic 
representation  of  the  M*G*B  model  and  the  actual  numerical  results  of  the 
prediction  method. 


The  computer  program  is  written  in  FORTRAN  IV  language.  It  has  been  run 
on  both  the  GE/Honeywell  6080  and  CDC  7600  computers,  and  can  easily  be  modi¬ 
fied  for  running  on  other  systems.  The  program  subdivides  the  jet  plume 
utilizing  a  built-in  grid  system  which  requires  minimal  input  for  specifi¬ 
cation.  This  grid  system  can  be  superseded  by  the  user  through  more  complex 
input  if  desired.  The  nozzle  geometry  is  input  through  discrete  point  coor¬ 
dinates  for  each  nozzle  element  boundary,  and  up  to  109  elements  can  be  input 
for  a  given  case.  A  maximum  of  24  axial  stations  along  the  jet  plume  is 
permitted,  and  up  to  200  radial  points  per  axial  station  can  be  accommodated. 
These  limits  can  be  changed  if  so  desired  by  modifying  the  appropriate 
DIMENSION  and  COMMON  statements  in  the  program  logic. 

The  limiting  assumptions  made  in  developing  the  method  have  been  dis¬ 
cussed  in  Reference  1,  but  it  is  appropriate  to  summarize  them  here  to  warn 
against  indiscriminate  violation  of  these  limitations.  They  are  as  follows: 

1.  The  exhaust  nozzle  elements  should  be  coplanar;  that  is,  each  tube  or 
chute  of  a  multielement  configuration  should  exhaust  at  the  same  axial 
plane.  However,  nozzle  element  exit  planes  can  be  staggered,  provided 
that  the  mixing  layer  of  a  given  element  does  not  impinge  on  the  wall 
of  another  element. 

2.  The  jet  exhaust  gases  must  all  be  of  the  same  constituent,  for 
the  calculation  cannot  accommodate  gas  mixtures  or  species  concen- 
t  rations . 

3.  Within  any  nozzle  element,  the  flow  is  assumed  to  be  uniform  at  the 
exit  plane. 

4.  The  time-averaged  static  pressure  is  assumed  to  be  constant  and 
uniform  throughout  the  jet  flow  field  and  surrounding  ambient 
field. 

5.  The  exhaust  nozzle  elements  must  discharge  axially,  radial  mean 
flow  and  swirl  are  neglected  in  the  model. 
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6.  The  effects  of  shock  formations  on  mixing  and  turbulence  levels  are 
neglected. 

These  assumptions  and  limitations  are  those  which  pertain  to  the  types 
of  problems  which  can  be  analyzed.  There  are,  of  course,  additional  assump¬ 
tions  that  went  into  the  formulation  of  che  model  itself  which  may  restrict 
the  accuracy  of  the  model,  but  which  do  not  restrict  the  type  of  problem 
which  can  be  analyzed.  The  user  is  advised  to  consult  Reference  2  for  the 
details  of  the  model  formulation. 


4.2  PROGRAM  NOMENCLATURE  AND  SYMBOL  CONVENTION 

Thu  jet  plume  and  nozzle  geometry  coordinates  are  computed  in  the  MAIN 
routine.  The  jet  plume  is  divided  into  KX  axial  slices,  specified  by  KA 
(1  £  KA  £  KX) .  The  FORTRAN  symbol  variables  for  the  various  coordinate 
parameters  and  indices  are  shown  in  Figure  4-1.  Note  that  the  radial  sub¬ 
division,  specified  by  index  M  (1  £  M  £  NR),  proceeds  in  increments  DSIG(KA), 
from  SIC  =  RMIN(KA)  to  the  maximum  value  set  by  NR.  The  value  of  NR  is 
determined  during  the  calculation  from  the  location  where  the  axial  momentum 
flux  is  within  a  certain  tolerance  of  being  equal  to  the  ambient  level,  i.e., 

| RU2  -  RU2E (1) |  £  RU2M 

where  RU2M  is  a  specified  input  tolerance.  The  maximum  allowable  value  of  NR 
can  be  specified  by  the  input  variable  IQUIT.  The  program  dimension  sizes 
Limit  KX  and  IQUIT  to  the  following  maximum  values: 

KX  £  24  IQUIT  £  200 

The  nozzle  geometry  itself  is  input  as  a  number  (NEST)  of  boundary  ele¬ 
ments.  Each  element  is  specified  by  coordinate  pairs  RA(I,J)  and  DALP(I,J), 
where  RA(I,J)  denotes  the  radius  and  DALP(I,J)  denotes  the  angular  increment, 
as  shown  in  Figure  4-1.  The  index  I  denotes  the  boundary  contour  point  number, 
and  the  index  J  denotes  the  boundary  number.  The  reference  angular  location 
for  each  boundary  is  given  by  ALPO(J).  For  each  boundary,  the  exit-plane 
values  of  total  pressure  PT(J)  and  total  temperature  TT(J)  are  also 
specified.  Boundary  Number  One  (J=l)  is  always  considered  to  be  the  ambient 
f ield. 

The  farfield  acoustic  calculations  are  performed  on  either  a  constant- 
radius  arc  or  a  sideline  parallel  to  the  jet  axis,  according  to  whether  the 
input  variable  NUMANG  is  set  equal  to  1  or  2,  respectively.  For  NUMANG  =  1, 
the  input  DIST  is  the  arc  radius;  for  NUMANG  =  2,  DIST  is  the  sideline  dis¬ 
tance.  The  acoustic  arena  geometry  specification  in  terms  of  FORTRAN  vari¬ 
ables  is  shown  in  Figure  4-2.  Note  that  a  distinction  is  made  between  the 
source-to-observer  distance  RSTAR  and  the  nozzle-to-observer  distance  RADIUS. 
The  observer  angle  relative  to  the  jet  axis  THETA  is  always  in  units  of 
radians,  while  the  observer  angle  relative  to  the  inlet  axis  THETD  is  in 
units  of  degrees.  The  farfield  sound  pressure  level  SPL(I,J)  is  computed  at 
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Figure  4-2.  FORTRAN  Symbol  Convention  for  Acoustic  Arena  Variables. 


every  1/3-octave  frequency  from  FM1N  to  FMAX,  at  ten-degree  increments  trout 
THETD  =  20°  to  160°. 

A  list  of  the  important  FORTRAN  symbols  used  in  the  computer  program  is 
given  in  Table  4-1,  along  with  their  algebraic  equivalents  where  possible.  A 
complete  description  of  all  of  the  input  variables  and  examples  of  input 
preparation  are  given  in  Section  4.5. 


4.3  DESCRIPTION  OF  PROGRAM  AND  SUBROUTINES 


A  flow  chart  of  the  computer  program  logic  is  shown  in  Figure  4-3.  It 
indicates  the  sequence  of  operations,  the  interconnections  of  various  portions 
of  the  program,  and  their  functions.  A  description  of  the  main  program  and 
each  of  the  subroutines  is  given  in  the  following  paragraphs. 


4.3.1  MAIN 

The  main  program  initiates  the  computation  and  controls  the  sequence  of 
operations.  It  reads  the  input  data,  computes  the  grid  system  for  the  aero¬ 
dynamic  flow  field,  and  computes  the  various  required  nozzle  exit  plane  flow 
parameters  such  as  velocities,  Mach  numbers,  momentum  and  enthalpy  fluxes, 
etc.  The  main  program  prints  out  all  input  data,  nozzle  exit  conditions, 
nozzle  geometry,  and  coordinate  system  parameters. 

The  main  program  controls  and  executes  the  jet  plume  flow  field  compu¬ 
tation.  After  each  axial  slice  has  been  evaluated,  the  MAIN  program  calls 
subroutine  SLICE  to  perform  the  requested  acoustic  calculations.  Upon  com¬ 
pleting  the  calculations  at  all  axial  slices,  MAIN  then  calls  OUTPUT  to  per¬ 
form  some  final  calculations  and  print  out  the  farfield  noise  levels.  If 
additional  cases  are  requested,  the  entire  procedure  is  repeated,  beginning 
with  reading  of  input  data;  otherwise  the  execution  is  halted. 


4.3.2  ARCCOS (X) 

This  is  a  function  subroutine  which  computes  the  principal  value  of  the 
arc  cosine  of  the  variable  X.  It  is  used  in  MAIN  in  evaluating  certain 
angles  relating  boundary  coordinate  points  and  flow  field  location  points. 


4.3.3  ERF (X) 

Tills  function  subroutine  evaluates  the  error  function  of  argument  X 
using  polynomial  approximations  as  given  in  Reference  3.  It  is  used  in 
MAIN  for  evaluating  flow  field  integrands. 
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Table  4-1.  List  of  FORTRAN  Symbols 


I'UHTHAN 

Symbol 

Meaning 

Related  Subroutines 

AA 

Air  attenuation  factor 

ATMOS 

AAA 

Intermediate  variable 

LSPFIT,  MAIN 

ABDTH 

|  A<J>| 

MAIN 

ABLE 

Intermediate  variable 

MAIN 

ABPA 

1  ‘P~a  | 

MAIN 

ACH 

Mach  number  M 

MAIN 

ACHM 

Average  mach  number 

MAIN 

ACH2 

M2 

MAIN 

AK 

Sound  level  constant 

MAIN,  OUTPUT 

AL 

Lighthill  parameter 

MAIN 

ALFA 

Frequency  constant 

MAIN 

ALP 

Angle 

MAIN 

ALPHT 

Convection  constant  a t 

SLICE 

ALP0 

Reference  boundary  angle 

MAIN 

AMU  IN 

Input  turbulence  constant  Pt 

MAIN 

AMULT 

Intermediate  value  for  pt 

MAIN 

AO 

Speed  of  sound  Ca 

MAIN 

ATOTAL 

Total  flow  area 

MAIN 

B 

Intermediate  variable 

LSPFIT 

BETA 

Shock  strength  parameter  8 

SHOCK 

BETAIN 

Input  turbulence  constant  Bt 

MAIN 

BETAMC 

Convection  constant  8mc 

MAIN,  SLICE 

BK 

Intermediate  variable 

SLICE 

BKR 

Intermediate  variable 

MAIN 

BOT 

Intermediate  variable 

LSPFIT 

BUG 

Intermediate  variable 

MAIN 

C 

Constant 

LSPFIT 

CH 

Spreading  parameter  C^/Cu, 

MAIN 

CHX 

Spreading  parameter  Chx 

MAIN 

C.TOCO 

Ratio  of  Cj/Ca 

SLICE 

CM 

Spreading  parameter  C^, 

MAIN 

UMAX 

Intermediate  variable 

TPNLC 

CMC 

Intermediate  variable 

MAIN 

CMMC 

Spreading  constant  C^ 

MAIN 

CMVR 

Spreading  constant  C2 

MAIN 

CNST 

Constant 

SLICE 

CO 

Ambient  speed  of  sound  Ca 

MAIN,  SLICE,  SHOCK 

COEF 

Conversion  factor 

OUTPUT 

CONV 

Convection  factor 

SHOCK 

CONVF 

Flight  dynamic  factor 

SLICE 

CONVO 

Convection  factor 

SLICE 

CONV  2 

Modified  convection  factor  C 

SLICE 

CON  1 

Constant 

SLICE 
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Table  4-1.  List  of  FORTRAN  Symbols  (Continued). 


FORTRAN 


Symbo 1 

Meaning 

Related  Subroutines 

C0N2 

Constant 

SLICE 

COST 

Cos  4> 

MAIN 

COSTO 

Cos  $ 

MAIN 

CP 

Specific  heat  Cp 

MAIN 

CT 

Cos  9 

SLICE,  CRD 

CTSQ 

Cos^  g 

SLICE 

CTH 

Cos  6 

SHOCK 

CVR 

Intermediate  variable 

MAIN 

DALP 

Boundary  coordinate  Aa 

MAIN 

DDTHE 

Tolerance  on  A0,  radians 

SLICE 

DDTHED 

Tolerance  on  A0 ,  degrees 

SLICE 

DELRA 

Transformed  boundary  radius  Av 

MAIN 

DELSIG 

Transformed  radius  Ar 

MAIN 

DELTA 

Turbulence  constant  <5t 

MAIN 

DELTIN 

Input  array  of  6t 

MAIN 

DEQ 

Equivalent  diameter  Deq 

MAIN,  SLICE,  SHOCK 

DIA 

Reference  Deq 

MAIN 

DIRECT 

Directivity  factor 

SLICE 

DIST 

Sideline  or  arc  distance 

MAIN,  SLICE 

DJET 

Reference  diameter 

MAIN 

DPHI 

At 

MAIN 

DRMIN 

Ar  -  minimum  value 

SLICE 

DS 

Source  strength  amplitude 

SLICE 

DSIG 

Ar 

MAIN,  SLICE 

DSPL 

Mixing  noise  pressure 

SLICE,  OUTPUT 

DSPL1 

Intermediate  variable 

SHOCK 

DSPL2 

Intermediate  variable 

SHOCK 

DTHED 

A0,  degrees 

SLICE 

DTHM 

Maximum  increment  of  t 

MAIN 

DU 

Intermediate  variable 

MAIN 

DUDR 

9U/3r 

MAIN,  SLICE 

DV 

Eddy  volume  dV 

SLICE 

DX 

Axial  step  size  Ax 

MAIN,  SLICE 

EF 

Enthalpy  flux 

MAIN 

EFE 

Enthalpy  flux 

MAIN 

EM 

Mach  number 

SLICE 

EMACH 

Exit  Mach  number 

MAIN,  SLICE,  OUTPUT 

F 

Intermediate  variable 

LSPFIT 

FAC 

Intermediate  variable 

PNLC 

FC 

Center  frequency 

SLICE 

FIRSTU 

Flight  velocity  Ua 

MAIN,  SLICE 

FIS 

Intermediate  variable 

MAIN 

FM 

Mass  flow 

MAIN 
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Table  4-1.  List  of  FORTRAN  Symbols  (Continued). 


FORTRAN 

Symbol 

Meaning 

Related  Subroutines 

FMAX 

Maximum  observed  frequency 

MAIN,  OUTPUT 

FMIN 

Minimum  observed  frequency 

MAIN,  OUTPUT 

FO 

Observed  frequency 

SLICE,  SHOCK,  OUTPUT 

FP 

Peak,  frequency 

SHOCK 

FR 

Frequency  ratio 

SLICE 

FRSQ 

Intermediate  variable 

SLICE 

FS 

Source  frequency 

SLICE 

GAM 

Specific  heat  ratio  y 

MAIN,  SHOCK 

GAMA 

Gas  constant  parameter 

MAIN 

GEXP 

Gas  constant  parameter 

SHOCK 

GM 

Shielding  function 

CRD 

GOSQ 

Shielding  function 

CRD 

G2 

Shielding  function 

SLICE,  CRD 

HF 

Spectrum  function 

SLICE 

HPSI 

Intermediate  variable 

MAIN 

HTR 

Stagnation  enthalpy 

MAIN 

1 

Index 

ALL 

IC 

Index 

LSPFIT 

I  DENT 

Title  (80-characters) 

MAIN 

11 

Index 

TPNLC 

IMH 

Index 

MAIN 

IQU1T 

Maximum  number  of  points 

MAIN 

IS 

Index 

MAIN 

ISSY 

Index 

MAIN 

1  SAVE 

Index 

LSPFIT 

ISYM 

Symmetry  indicator 

MAIN 

IT 

Symmetry  indicator 

MAIN 

J 

Index 

ALL 

JMAX 

Maximum  band  number 

OUTPUT,  SHOCK,  SLICE 

JM1N 

Minimum  band  number 

OUTPUT,  SHOCK,  SLICE 

J1 

Index 

CRD 

Jll 

Index 

CRD 

J2 

Index 

CRD 

J21 

Index 

CRD 

J211 

Index 

CRD 

K 

Index,  also  wave  number 

MAIN,  SLICE,  PNLC 

KN 

Surrounding  boundary  index 

MAIN 

KNCAS 

Case  counter 

MAIN 

KNK 

Surrounding  boundary  index 

MAIN 

KX 

Number  of  axial  slices 

MAIN 

L 

Index 

MAIN 

LAVG 

Shock  spacing 

SHOCK 

LEAF 

Number  of  boundary  leaves 

MAIN 
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Table  4-1.  List  of  FORTRAN  Symbols  (Continued). 


FORTRAN 

Symbol 

Meaning 

Related  Subroutines 

LEAV 

Number  of  boundary  leaves 

MAIN 

LINE 

Printout  counter 

MAIN 

LPHI 

Number  of  flow  field  leaves 

MAIN 

LQ 

Index 

MAIN 

M 

Index 

MAIN,  SLICE 

MACH 

Mach  number 

SLICE 

MAXNOY 

Maximum  noy  value 

PNLC 

MC 

Convection  Mach  number 

SLICE,  SHOCK,  CRD 

MCIN 

Input  array  of  Mc 

SLICE 

MIN 

Input  array  of  Mc 

CRD 

M.I 

.let  exit  Mach  number 

SHOCK 

N 

Index,  also  number  of  shocks 

MAIN,  SHOCK,  LSPFIT 

NBREF 

Reference  boundary  number 

MAIN 

NCASK 

Number  of  cases 

MAIN 

NCBDY 

Number  of  centerbody  points 

MAIN 

NCELL 

Number  of  shock  cells 

MAIN,  SHOCK 

N COUNT 

Counter 

LSPFIT 

NN 

Acoustic  calculation  selector 

MAIN,  SLICE 

NODE 

Intermediate  variable 

MAIN 

NOV 

Minimum  number  of  points 

MAIN 

NOY 

Noy  value 

PNLC 

NPAOK 

Page  counter 

MAIN 

NPR 

Printout  counter 

MAIN 

NPRINT 

Printout  selector 

MAIN,  SLICE 

NPTS 

Number  of  points 

LSPFIT 

NR 

Number  of  points 

SLICE,  CRD 

NRI 

I  ndex 

SLICE 

NTP 

Number  of  turning  points 

SLICE,  CRD 

NUM 

Number  of  boundary  points 

MAIN 

NUMANil 

Arena  selector 

MAIN,  SLICE 

NHMK 

Number  of  boundary  points 

MAIN 

nxc 

Index 

LSPFIT 

OAPWI. 

Overa 1 1  powe r  1  eve  1 

OUTPUT 

OASPL 

Overall  sound  pressure  level 

OUTPUT,  PNLC 

OBSTN 

Observed  Strouhal  number 

OUTPUT 

0MEC1R 

Source  radian  frequency 

SLICE 

PAA 

Ambient  static  pressure 

MAIN 

PC 

Intermediate  variable 

PNLC 

PCC 

(las  constant  parameter 

MAIN 

PHI 

An  g  1  e  {> 

MAIN 

PI 

Constant  n 

MAIN,  SLICE,  OUTPUT 

P 102 

"/ 2 

CRD 

P 12 

2n 

MAIN 
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Table  4-1.  1. 1st  of  FORTRAN  Symbols  (Cent I lined ) . 


FORTRAN 

Symbol 

Meaning 

Related  Subroutines 

PNDB 

PNdB 

PNLC 

PNL 

PNL 

OUTPUT,  PNLC 

PNLT 

PNLt,  tone-corrected  PNL 

OUTPUT 

POWER 

Exponent 

MAIN 

PS 

Ambient  static  pressure 

MAIN,  SHOCK 

PSQ 

Square  of  acoustic  pressure 

OUTPUT 

PSQM 

Mixing  noise  p2 

SHOCK 

PSQS 

Shock  noise  p2 

SHOCK 

PSQT 

Total  noise  p^ 

SHOCK 

PT 

Stagnation  pressure 

MAIN,  SHOCK 

PWL 

Power  level 

OUTPUT 

PWR 

Sound  power,  watts 

OUTPUT 

Q 

Intermediate  variable 

MAIN 

RA 

Boundary  coordinate  radius 

MAIN 

RAD 

Flow  integration  variable  R0 

MAIN 

RADO 

Flow  integration  variable  R0 

MAIN 

RADIUS 

Nozzle-to-observer  radius  R 

SLICE,  OUTPUT,  ATMOS 

RADX 

Argument  R0/Cmx 

MAIN 

RCBDY 

Centerbody  radial  cooordinate 

MAIN 

PRCRIT 

Critical  pressure  ratio 

SHOCK 

RCRC 

Intermediate  variable 

MAIN 

RFO 

Intermediate  variable 

OUTPUT 

RHO 

Density  p 

MAIN 

RHOE 

Ambient  density  pa 

MAIN,  SLICE 

RHOESQ 

pa 

Az imuthally-averaged  p 

SLICE 

RHOR 

MAIN,  SLICE 

RIN 

Input  radius 

SLICE,  CRD 

RJET 

Reference  jet  density  ratio 

MAIN 

RMIN 

Minimum  value  of  r 

MAIN 

RMINEX 

Exit  plane  value  of  RMIN 

MAIN 

RMINSQ 

Square  of  RMIN 

MAIN 

RMNSQE 

Square  of  RMINEX 

MAIN 

RMP 

Dummy  variable 

MAIN 

RND 

ROOT 

Normalized  radius  r/Deq 

Root  (zero)  of  g^ 

MAIN 

SLICE 

R00T2 

Si 

SLICE 

RO 

Source  radius  r0 

CRD 

RSIG 

Turning  point  radius  rCT 

SLICE,  CRD 

RSKJl 

rol 

CRD 

RSIC2 

ro2 

CRD 

RSORSQ 

Source  location  correction  (R*/R)2 

SLICE 

RSTAR 

Source-to-observer  radius  R* 

SLICE 

105 


Table  4-1.  List  of  FORTRAN  Symbols  (Continued) 


FORTRAN 

Symbol 

Meaning 

Related 

Subroutines 

RU 

Mass  flux  pU 

MAIN 

RU2 

Momentum  flux  plj2 

MAIN 

RU2E 

Exit  plane  value  of  plj2 

MAIN 

RU2M 

Minimum  value  of  plj2 

MAIN 

RU2REF 

Reference  value  of  pl|2 

MAIN 

R1 

Intermediate  variable 

CRD 

R2 

Intermediate  variable 

CRD 

S 

Intermediate  variable 

TPNLC 

SA 

Intermediate  variable 

MAIN 

SAC 

Intermediate  value  of  x^ 

MAIN 

SACO 

Intermediate  value  of  x^ 
Intermediate  value  of  xr 

MAIN 

SAR 

MAIN 

SARO 

Intermediate  value  of  xr 

MAIN 

SAX 

Intermediate  value  of  xx 

MAIN 

SAXO 

Intermediate  value  of  xx 

MAIN 

SBAR 

Intermediate  variable 

TPNLC 

SDl) 

Intermediate  value  of  3U/3r 

MAIN 

SEFE 

Integral  of  enthalpy  flux 

MAIN 

SGN 

Sign 

LSPFIT 

SGN1 

Sign 

CRD 

SGN2 

Sign 

CRD 

SGI 

Intermediate  variable 

CRD 

SG2 

Intermediate  variable 

CRD 

SHIELD 

Shielding  integral 

SLICE, 

CRD 

SIC 

Intermediate  value  of  x^ 

Radius  r 

MAIN 

SIG 

MAIN 

SIGN 

Sign 

ERF 

SIGSQ 

r2 

MAIN 

SIGR 

Radius  r 

MAIN,  SLICE 

SINT 

Sin  6 

MAIN 

SINTO 

Sin  0O 

MAIN 

SIR 

Intermediate  value  of  !r 

MAIN 

SIX 

Intermediate  value  of  xx 

MAIN 

SPL 

SPL  array 

ALL 

SPLL 

Intermediate  variable 

TPNLC 

SPLMAX 

Maximum  SPL 

SHOCK 

SPLP 

Intermediate  variable 

TPNLC 

SPLPP 

Intermediate  variable 

TPNLC 

SPLU 

Intermediate  variable 

TPNLC 

SRU 

Mass  flux  integral 

MAIN 

SRUM 

Mass  flux  integral 

MAIN 

SRU2 

Momentum  flux  integral 

MAIN 

SRU2M 

Momentum  flux  integral 

MAIN 
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List  of  FORTRAN  Symbols  (Continued) 


t‘,:> )  R/N 

Sy mbo 1 

Meaning 

Related  Subroutines 

SS 

SPL  array 

PNLC 

SSPL 

Shock  noise  SPL  array 

SHOCK 

STC 

Azimuthal  shear  stress  x ^ 

MAIN 

STR 

Radial  shear  stress  xr 

MAIN 

STRFR 

Radial  coordinate  stretching  factor 

MAIN 

STRI'X 

Axial  coordinate  stretching  factor 

MAIN 

SIX 

Axial  shear  stress  tx 

MAIN 

SUE 

Reference  velocity 

MAIN 

SOM 

Sum 

OUTPUT 

SUMNOY 

Sum  of  noy  value 

PNLC 

SUMSPL 

Sum  of  SPL 

PNLC 

SOM  l 

Sum 

CRD 

S0M2 

Sum 

CRD 

SU8 

Integral  of  source  strength 

MAIN 

SU8M 

Integral  of  source  strength 

MAIN 

SV2 

Square  of  velocity 

MAIN 

SI 

Intermediate  variable 

LSPFIT 

T 

Temperature 

ERF,  MAIN 

TA 

Intermediate  variable 

MAIN 

TAA 

Ambient  static  temperature 

MAIN 

TAO 

Intermediate  variable 

MAIN 

TAU 

Total  shear  stress  t 

MAIN 

TAUR 

Azimuthal  average  of  r 

MAIN,  SLICE 

TC2 

Intermediate  variable 

TPNLC 

TC3 

Intermediate  variable 

TPNLC 

TE 

Exit  static  temperature 

MAIN 

TEMP 

Normalized  temperature  T/Ta 

SLICE 

TERM 

Directivity  factor 

SLICE 

TH 

Angle 

MAIN 

THOR 

Critical  angle  6cr 

SHOCK 

TERM 

Directivity  factor 

SLICE 

THE 

Angle  H 

SLICE,  CRD 

THETA 

Observer  angle  6,  radians 

SLICE,  OUTPUT 

THETI) 

Observer  angle  Oj,  degrees 

SLICE,  OUTPUT,  SHOCK 

THO 

£o 

MAIN 

THT 

Observer  angle  0j,  radians 

SHOCK 

Tl 

Intermediate  value  of  enthalpy 

flux 

MAIN 

TOP 

Intermediate  variable 

LSPFIT 

TSR 

Static  temperature 

MAIN 

TSTD 

Circumferential  asymmetry  test 

parameter 

MAIN 

TSTtl 

Circumferential  asymmetry  test 

parameter 

MAIN 

IS.  lit. 

Clrcumferenti.il  asymmetry  test 

parameter 

MAIN 

TSTL 

Circumferential  asymmetry  test 

parameter 

MAIN 
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Table  4-1 


List  of  FORTRAN  Symbols  (Concluded). 


FORTRAN 

Symbol 

Meaning 

Related  Subroutines 

TT 

Stagnation  temperature 

MAIN 

TTR 

Stagnation  temperature 

MAIN 

TURBIN 

Turbulence  intensity,  u' 

MAIN 

U 

Mean  velocity 

MAIN 

UAP 

Intermediate  variable 

MAIN 

UAVG 

Mass-average  of  U  at  x 

MAIN 

UC 

Convection  velocity  Uc 

SHOCK 

UE 

Exit  plane  velocity  Uj 

MAIN,  SHOCK 

UGLY 

Intermediate  variable 

MAIN 

UJET 

Reference  exit  velocity 

MAIN 

UMAX 

Maximum  local  value  of  U  at  x 

MAIN 

UND 

Normalized  value  of  U,  U/UREF 

MAIN 

UNITS 

Constant  for  units  conversion 

MAIN,  OUTPUT 

UR 

Azimuthal  average  of  U 

MAIN,  SLICE 

UREF 

Reference  exit  velocity 

MAIN 

U8 

Intermediate  value  of  source  strength 

MAIN 

U8I 

Integral  of  source  strength 

MAIN 

VA 

Intermediate  value  of  momentum 

MAIN 

VAO 

Intermediate  value  of  momentum 

MAIN 

VI 

Intermediate  value  of  momentum 

MAIN 

VMAX 

Maximum  of  velocities  inside  and  outside 

MAIN 

VMIN 

Minimum  of  velocities  inside  and  outside 

MAIN 

VO 

Flight  velocity  Ua 

SHOCK 

VR 

Velocity  ratio  VMIN/VMAX 

MAIN 

WITHIN 

Dummy  variable 

LSPFIT 

X 

Axial  distance  x 

MAIN,  SLICE 

XCBDY 

Centerbody  axial  coordinate 

MAIN 

XD 

Intermediate  variable 

LSPFIT 

XE 

Exit  plane  axial  coordinate 

MAIN 

XND 

Normalized  axial  coordinate  X/Deq 

MAIN 

XOR 

Variable  x/R 

SLICE 

XI 

Intermediate  variable  for  curve  fitting 

LSPFIT 

XI 3 

Intermediate  variable  for  curve  fitting 

LSPFIT 

X4 

Intermediate  variable  for  curve  fitting 

LSPFIT 

X43 

Intermediate  variable  for  curve  fitting 

LSPFIT 

Y 

Intermediate  variable  for  curve  fitting 

LSPFIT 

YC 

Intermediate  variable  for  curve  fitting 

LSPFIT 

YI 

Intermediate  variable  for  curve  fitting 

LSPFIT 

Y3 

Intermediate  variable  for  curve  fitting 

LSPFIT 
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Figure  4-3.  Computer  Program  Flow  Chart. 
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(b)  Subroutine  SLICE 


Figure  4-3.  Computer  Program  Flow  Chart  (Continued) 
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(c)  Subroutine  OUTPUT 

Figure  4-3.  Computer  Program  Flow  Chart  (Concluded) 


4.3.4  LSPFIT 


Subroutine  LSPFIT  is  a  curve-fitting  routine  which  utilizes  least- 
squares  polynomial  fits  of  second  order  to  perform  interpolation,  differentia¬ 
tion  and  integration  of  input  discrete-point  data.  The  calling  statement  is: 

CALL  LSPFIT (X,  Y,  N,  XC,  YC,  NC,  NF,  A) 

where  (X,  Y)  are  the  input  data  coordinates  (N  values  of  each),  XC  are  the 
values  of  X  where  output  is  requested,  YC  are  the  output  functions,  NC  is  the 
number  of  output  data  points,  and  NF  indicates  the  type  of  output  desired. 

The  coding  for  NF  is  as  follows: 


NF  =  0,  YC  are  interpolated  values  of  Y 

NF  =  1,  YC  are  derivatives  of  Y 

NF  =  -1,  YC  is  the  integral  of  Y  from 
XC(1)  to  XC(J),  1  <  J  <  NC. 


The  parameter  A  is  the  second  derivative  of  Y.  Subroutine  LSPFIT  is  used  in 
MAIN  to  interpolate  input  plug/centerbody  geometry  coordinates  at  various 
axial  stations  in  the  flow  field,  and  to  obtain  radial  gradients  of  density 
from  the  computed  density  profiles. 

4.3.5  SLICE 

Subroutine  SLICE  directs  the  mixing  noise  calculation  for  each  axial 
slice.  The  calling  sequence  is  as  follows: 

CALL  SLICE  (X (KA) ,  DSIG(KA),  DX,  M) 


where  X ( KA )  is  the  axial  location,  DSIG(KA)  is  the  radial  step  size,  DX  Is 
the  axial  slice  thickness,  and  M  is  the  number  of  radial  points  in  the  slice. 
The  flow  parameters  (which  are  circumferentially  mass-averaged  values)  are 
transferred  through  labeled  COMMON  statements.  Subroutine  SLICE  computes  the 
acoustic  arena  geometry  parameters  THETA,  THETD,  RADIUS  and  initializes  SPL 
( I » -J )  to  zero  during  the  first  call,  skipping  this  calculation  on  succeeding 
calls.  The  normalized  radial  profiles  of  velocity  (MACH)  and  temperature 
(TEMP)  are  evaluated,  followed  by  a  calculation  of  source  strength  amplitude 
DS  and  characteristic  frequency  FS  for  each  radial  volume  element. 
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Subroutine  SLICE  computes  the  acoustic  shielding  function  profiles  G2(J), 
the  number  of  turning  points  NTP ,  and  their  locations  RSIG.  Subroutine  CRD 
is  then  called  to  calculate  the  acoustic  shielding  exponentials  and  quadrupole 
directivity  functions.  Subroutine  SLICE  then  sums  up  the  mixing  noise  con¬ 
tributions  from  each  radial  volume  element,  factoring  in  their  individual 
source  strengths,  characteristic  frequencies,  spectrum  shapes,  directivities, 
and  shielding  factors.  The  resulting  noise  spectrum  from  each  slice  is  stored 
as  the  variable  DSPL(I,J),  where  I  denotes  the  observer  angle  index  and  J  is 
the  1/3-octave  frequency  band  index.  Upon  completing  the  calculation  for  a 
given  slice,  SLICE  returns  control  to  MAIN. 


4 . 3 .  b  CRD 

Subroutine  CRD  computes  the  shielding  function  integrals  and  quadrupole 
directivity  factors  for  a  given  axial  slice  as  a  function  of  radial  source 
location.  The  radial  distributions  of  normalized  velocity  (MACH)  and 
temperatures  (TEMP)  and  shielding  function  (G2)  are  transferred  to  CRD  through 
labeled  COMMON  statements.  The  calling  statement  is: 

CALL  CRD 


At  each  source  radius,  subroutine  CRD  interrogates  the  data  to  determine 
which  of  the  six  shielding  conditions  in  Figure  4-4  applies,  and  computes  the 
appropriate  shielding  integral  (8oi»  S02>  or  Rl2)  -ind  the  appropriate 
directivity  factors.  After  all  radial  source  volumes  have  been  evaluated, 

CRD  returns  control  to  SLICE. 


4.3.7  OUTPUT 

Subroutine  OUTPUT  performs  the  final  acoustic  calculations  and  prints 
out  tiie  far  field  SPI.  spectra,  OASPL,  PNL  and  PNLT  directivities.  The  calling 
sequence  is  as  follows: 


CALL  OUTPUT  (EMACH,  DJET,  RJET,  UJET,  UNITS) 


where  EMACH,  DJET,  RJET,  and  UJET  are  the  characteristic  (usually  reference) 
jet  Mach  number,  diameter,  density  ratio  and  velocity,  respectively.  The 
parameter  UNITS  is  a  conversion  factor  for  converting  from  lbf/ft2  to 
dynes/cm^  relative  to  0.0002  dynes/cm2.  Subroutine  OUTPUT  converts  the 
narrowband  spectra  from  SLTCE  into  1 /3-octave  levels.  Subroutine  SLTCE  then 
calls  SHOCK  to  compute  SSPL  spectra  (shock  noise)  and  adds  these  to  the 
turbulent  mixing  noise  spectra  to  obtain  the  total-noise  spectra.  The 
corresponding  power  spectrum  (PWI.)  is  then  computed,  and  subroutine  ATMOS  is 
1  hen  called  to  correct  all  SPL  spectra  for  atme.pheric  attenuation.  Sub¬ 
routines  PNI.C  and  TPN1.C  are  then  called  to  calculate  perceived  noise  level 
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PNL  and  tone-corrected  noise  level  PNLT.  Finally,  overall  sound  pressure 
level  OASPL  is  computed,  and  all  of  those  acoustic  parameters  are  then 
printed  out.  Subroutine  OUTPUT  then  returns  control  to  MAIN. 


4.3.8  SHOCK 

Subroutine  SHOCK  computes  the  broadband  shock-associated  noise  spectra 
at  each  observer  angle.  The  calling  statement  is  as  follows: 


CALL  SHOCK 


All  parameters  are  transferred  into  and  out  of  this  subroutine  through 
labeled  COMMON  statements.  Subroutine  SHOCK  computes  the  1/3-octave  SPL 
spectra  for  each  nozzle  boundary  element  which  has  a  non-zero  shock  cell 
number  input,  NCELL  >  0.  The  individual  boundary  contributions  are  summed  on 
a  mean-square  pressure  basis  and  added  to  the  mixing  noise  spectra. 


4.3.9  ATMOS 

Subroutine  ATMOS  corrects  the  input  SPL  spectra  for  atmospheric  attenu¬ 
ation  effects  using  standard-day  atmospheric  absorption  factors  for  70% 
relative  humidity  and  59°  F  ambient  conditions.  The  calling  sequence  is  as 
follows : 


CALL  ATMOS  (SPL,  RADIUS) 


where  SPL(1,J)  is  the  sound  pressure  spectrum  array,  I  denotes  the  index  for 
observer  angle,  J  denotes  the  index  on  frequency,  and  RADIUS(I)  is  the 
nozzle-to-observer  distance  array.  The  atmospheric  absorption  in  dB  per  1000 
ft,  from  Reference  4,  is  corrected  to  the  proper  distance  RADIUS(J),  and  the 
result  is  subtracted  from  SPL(I,J).  The  array  of  SPL(I,J)  returned  to  OUTPUT 
is  the  corrected  array. 


4.3.10  PNLC 

Subroutine  PNLC  computes  the  perceived  noise  level  in  PNdB  at  each 
observer  angle  from  the  input  1/3-octave  spectra.  The  calling  sequence  is  as 
f  o 1 1 ows : 


CALL  PNLC  (SS,  FAC,  PNDB,  OASPL) 
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where  SS  is  the  input  array  ol  either  1/J-oclave  or  octave  SPL  values,  PAG  is 
a  constant  equal  to  0,15  for  1/3-octave  and  0.3  for  octave  levels,  PNDB  is 
the  output  PNL,  and  OASPL  is  the  conventional  overall  level.  The  method  used 
to  calculate  PNL  is  taken  from  Reference  5.  The  OASPL  output  from  sub¬ 
routine  PNLC  is  discarded  because  it  only  computes  the  summation  for  the 
first  24  values  of  SS.  Tills  Is  sometimes  Insufficient  for  scale  model 
condition,  where  the  frequency  range  of  interest  can  cover  as  many  as  thirty- 
three  1/3-octave  frequency  bands. 


4.3.11  TPNLC 

Subroutine  TPNLC  determines  a  pure-tone  correction  factor  to  the  PNL 
value  as  a  function  of  the  1/3-octave  SPL  spectrum.  The  calling  sequence  is 
as  follows: 


CALL  TPNLC  (SPL,  PTCOR) 


where  SPL  is  the  input  1/ 3-octave  spectrum  and  PTCOR  is  the  correction  to  be 
applied  to  PNL  to  account  for  the  presence  of  tones  in  the  spectrum.  Sub¬ 
routine  TPNLC  reads  in  SPL  and  returns  PTCOR.  The  tone  correction  and 
detection  procedure  is  based  on  the  method  proposed  in  Reference  7. 


4.4  PROGRAM  USAGE  AND  LOGIC 


A  complete  description  of  the  program  input  variables  and  input  format 
is  given  in  Section  4.5.  A  list  of  notes  and  suggestions  on  running  the 
program  is  also  included.  A  description  of  program  output  format,  including 
warning  flags  and  diagnostics,  is  given  in  Section  4.6.  A  sample  case  listing 
(including  input  data  card  images)  is  given  in  Section  4.7  for  a  7-tube 
suppressor  nozzle,  one  of  the  data-theory  comparison  cases  presented  in 
Reference  2.  A  complete  FORTRAN  source  listing  of  the  program  logic  is 
given  in  Section  4.8. 

Program  users  should  be  completely  familiar  with  Appendix  A,  since  there 
are  many  pitfalls  which  can  be  avoided  by  giving  attention  to  the  recom¬ 
mendations  presented  therein.  The  program  flexibility  permits  analysis  of 
nozzle  planforms  of  any  imaginable  shape,  so  long  as  certain  input  rules  and 
guidelines  are  followed.  When  non-axisymmetric  nozzles  are  run,  a  completely 
three-dimensional,  turbulent,  compressible  flow  field  analysis  is  performed, 
and  input  mistakes  can  be  costly  in  terms  of  computer  processor  time.  The 
user  should  make  initial  checkout  runs  for  complex  nozzles,  running  just  one 
or  two  axial  slices  at  first,  to  ensure  that  all  input  is  as  desired,  before 
running  a  complete  jet  plume. 

The  program  is  designed  to  serve  as  a  diagnostic  tool,  in  addition  to 
functioning  in  the  standard  jet  noise  prediction  mode.  Individual  slice 
calculations  can  be  made  by  suitable  input  selection,  running  each  slice  (or 
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axial  station)  as  a  separate  case.  This  mode  permits  evaluation  of  the 
relative  contributions  of  each  slice  at  each  frequency  and  observer  angle. 
Various  components  of  the  acoustic  model  can  be  bypassed  to  assess,  for 
example,  the  separate  effects  of  convection,  acoustic  shielding,  etc.  The 
program  can  also  be  used  to  predict  only  the  jet  flow  field,  if  desired. 
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4.5  DESCRIPTION  OF  INPUT 


The  input  data  is  supplied  through  NAMELIST  input  format,  with  the 
exception  of  the  alphanumeric  title  data  card,  which  precedes  the  input 
NAMELIST  data.  Any  number  of  successive  cases  can  be  run  consecutively, 
limited  only  by  the  user's  execution  time  available.  Each  successive  case 
requires  a  title  card  (80  -  character  label  in  columns  1  -  80),  followed  by 
the  INPUT  NAMELIST.  The  data  from  preceding  cases  remain  in  storage,  so 
only  those  variables  which  are  to  be  changed  from  the  preceding  case  input 
value  need  be  included  in  the  INPUT  file  of  succeeding  cases. 

A  suggested  input  preparation  format  is  given  in  Table  4-2.  Those 
variables  marked  by  an  asterisk  (*)  have  preset  values  built  into  the  program, 
and  need  not  be  input  unless  the  user  desires  to  override  the  preset  values 
with  a  different  one.  The  definitions  of  each  of  the  input  variables  given  in 
Table  4-2  are  listed  in  Table  4-3.  Again,  preset  variables  are  marked  by  an 
asterisk  (*).  The  values  of  those  variables  which  are  preset  are  given  in 
Table  4-4.  The  format  of  Table  4-3  is  such  that  a  note  number  (where  approp¬ 
riate)  is  given  for  each  variable  which  corresponds  to  the  note  number  in 
Section  4.5.1  ("Notes  on  Input").  These  notes  give  further  elaboration  on  how  to 
specify  and  prepare  the  input  data. 
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Table  4-2.  Suggested  Input  Format 


Column 

2 

(80  -  CHARACTER  TITLE  CARD,  COLUMNS  1-80) 

$ INPUT 

KX*  =  _ ,  NEST  =  _ ,  LPHI  *  =  _ ,  ISYM  =  _ 

IQUIT*  =  _ ,  NN*  =  _ ,  NCASE*  ■=  _ ,  NBREF* 

NPRINT*  =  _ ,  NCBDY  =  _ , 


N0V  = 

X 

DSIG  = 

BETAIN* 

DELTIN* 

AMUIN* 

RMIN 


XE  =  0 

ALP0  =  0 

LEAV  =  1 

NUM  =  1 

KN  1 


DEQ  = 
DS 

NCELL  = 

PI¬ 

TT 
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Table  4-2.  Suggested  Input  Format  (Concluded) 


Column 

2 


DALP(1 ,2)  = 

J 

9  9 

9  9 

DALP(1 ,3)  = 

9 

9  9 

9  9 

(etc .  ,  for 

boundary  4,  5,  6,  .. 

. . .NEST) 

RA(1,  2)  = 

9 

9  9 

9  9 

RA(1,  3)  = 

* 

9  9 

9  9 

(etc . ,  for 

boundary  4,  5,  6,  .. 

. .  .Nest) 

CM*  = 

,  CH*  = 

,  CMVR*  = 

,  CMMC*  = 

9 

CAM  = 

,  CP  = 

,  PS  = 

,  ALFA*  = 

9 

DTHM  *  = 

,  RU2M*  = 

,  AK*  = 

,  BK*  = 

* 

STRFR*  = 

,  STRFX*  = 

,  AT0TAL  = 

9 

ALPHMC*  = 

,  BETAMC *  = 

V 

NUMANG  = 

,  DIST  = 

,  FMIN  *  = 

,  FMAX*  = 

ALPHT*  = 

9  9 

9 

9  9 

XCBDY  = 

9  9 

9 

9  9 

9 

RCBDY  = 

9  9 

9 

9  9 

9 

$ 

(NEXT  CASE,  IF  ANY) 
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Table  4-3.  Input  Variable  Definitions. 


Variable 

Note 

Description 

KX* 

Number  of  axial  stations  to  be  analyzed;  a 
maximum  of  24  stations  is  permitted. 

NEST 

1 

Number  of  closed  boundary  contours  defining 
the  nozzle  exit  geometry;  a  maximum  of  110 
is  permitted. 

LPHI 

7 

Number  of  symmetric  leaves  (repeating  seg¬ 
ments  in  the  nozzle  exit  planform. 

ISYM 

Nozzle  symmetry  indicator;  ISYM  =  1  for  ax- 
symmetric  nozzles  or  completely  asymmetric 
nozzles,  =  0  otherwise. 

IQUIT 

Maximum  number  of  radii  at  which  flow  field 
is  calculated  (£200) . 

NN* 

12 

Acoustic  Calculation  option  indicator. 

NCASE* 

Number  of  cases  to  be  run  consecutively. 

NBREF* 

Reference  condition  boundary  number. 

NPRINT* 

13 

Aerodynamic  station  printout  indicator. 

NCBDY 

9 

Number  of  centerbody  input  coordinate  points 
A  maximum  of  40  is  permitted. 

N0V 

Minimum  number  of  radii  at  which  flow  field 
is  to  be  calculated,  for  each  axial  station 
(KX  values  required). 

X 

11 

Axial  location  of  each  axial  station,  ft. 

(KX  values' required) . 

DSIG 

11 

Radial  step  size  to  be  taken  for  flow  field 
calculation  at  each  axial  station,  ft. 

(KX  values  required). 

BETA IN* 

15 

Axial  shear  stress  turbulence  constant 
(KX  values  required). 

DELTIN* 

Azimuthal  shear  stress  turbulence  constant 
(KX  values  required). 
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Table  4-3.  Input  Variable  Definitions.  (Continued) 


Variable 

Note 

Description 

AMUIN* 

Azimuthal  velocity  gradient  turbulence  fre¬ 
quency  constant  (KX  values  required) . 

RMIN 

9 

Minimum  radius  for  flow  field  calculation 
at  each  axial  station  (KC  values  required). 

XF. 

8 

Axial  location  of  exit  plane  of  each 
boundary,  ft.  (NEST  values  required). 

ALP0 

2 

Reference  angle  aQ  from  which  the  coordi¬ 
nates  of  each  boundary  point  are  specified, 
radians  (NEST  values  required) . 

LEAV 

1,4 

Number  of  symmetric  leaves  (repeating  seg¬ 
ments)  of  each  boundary  (NEST  values 
required) . 

NUM 

1,5 

Number  of  input  points  (coordinate  pairs) 
to  be  supplied  for  each  boun.ary  (NEST 
values  required). 

KN 

1 

The  number  of  the  boundary  which  encloses 
a  given  boundary  (NEST  values  required). 

DEQ 

16 

Equivalent  flow  area  diameter  of  each 
boundary,  ft.  (NEST  values  required). 

DS 

16 

Shock-cell  spacing  characteristic  dimension, 
usually  hydraulic  diameter,  of  each  boun¬ 
dary,  ft.  (NEST  values  required). 

NCELL 

16 

Number  of  shock  cells  for  each  boundary 
element  (NEST  values  required). 

PT 

6 

Stagnation  pressure  inside  each  boundary, 
lbf/ft  (NEST  values  required). 

TT 

6 

Stagnation  temperature  inside  each  boun¬ 
dary  °  R  (NEST  values  required) . 

DALP(I,.J) 

2,3,5 

Angular  increment  Aa  from  preceding  boun¬ 
dary  point  which  locates  the  given  boun¬ 
dary  point  I  on  boundary  J,  radians  (omit 
boundary  number  1,  ambient  field). 
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Table  4-3.  Input  Variable  Definitions  (Continued). 


Variable 

Note 

Desci Lpl ion 

RA  (I,J) 

2,3,5 

Radial  coordinates  of  boundary  point  I  on 
boundary  J,  ft.  (omit  boundary  number  1, 
ambient  field). 

CM* 

10 

Empirical  jet  momentum  diffusion  rate 
spreading  parameter  C^j . 

CH* 

10 

Ratio  of  enthalpy- to-momen turn  spreading 
parameters  Cj, / Cm . 

CMVR* 

10 

Momentum  spreading  parameter  velocity 
ratio  influence  coefficient. 

CMMC* 

10 

Momentum  spreading  parameter  Mach  number 
influence  coefficient. 

GAM 

Specific  heat  ratio  y  =  Cp/Cv. 

CP 

Specific  heat  at  constant  pressure  Cp, 
in  (f t-lbf) /(slug  -  °  R) 

PS 

Ambient  static  pressure,  lbf/ft  . 

ALFA* 

Turbulence  characteristic  frequency  con¬ 
stant. 

DTHM  * 

7 

Maximum  allowable  increment  in  angular 
coordinate,  (d<j>)max«  for  flow  field  calcula¬ 
tion. 

RU2M  * 

2 

Minimum  value  of  jet  momentum  flux,  (pU  )  .  , 
below  which  the  flow  is  not  calculated. 

AK* 

Sound  pressure  level  proportionality  con¬ 
stant  for  mixing  noise  calculation. 

BK* 

Sound  pressure  level  proportionality  con- 
for  dipole  density-gradient  noise  calcula- 
t  ion. 

STRFR  * 

11 

Radial  coordinate  stretching  factor  for  use 
of  automatic  mesh  calculation. 

STRFX  * 

11 

Axial  coordinate  stretching  factor  for  use 
of  automatic  mesh  calculation. 
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Table  4-3.  Input  Variable  Definitions  (Concluded). 


Variable 

Note 

Description 

AT0TAL 

2 

Nozzle  Total  exit  flow  area,  ft  . 

ALPHMC  * 

14 

Convection  Mach  number  weighting  factor. 

BETAMC* 

14 

Convection  Mach  number  weighting  factor. 

NUMANG 

Arena  selection  indicator;  NUMANG  =  1  indi¬ 
cates  constant  radius  arc,  NUMANG  =  2 
indicates  sideline  parallel  to  the  jet 
axis. 

DIST 

Arc  or  sideline  distance,  ft. 

FMIN  * 

Minimum  frequency  for  which  acoustic  cal¬ 
culations  are  required,  Hz  (£50);  an 
integer  variable. 

FMAX* 

Maximum  frequency  for  which  acoustic  calcu¬ 
lations  are  required,  Hz  (£100,000);  an 
integer  variable. 

ALPHT* 

Convective  amplification  factor  turbulence 
constant  at;  15  values  required,  one  for 
each  observer  angle  0j  from  9j  =  20°  to  160 
in  10°  increments. 

XCBDY 

9 

Centerbodv  input  point  axial  coordinate, 
NCBDY  values  required. 

RCBDY 

9 

Centerbody  input  point  radial  coordinate, 

NCBDT  values  required. 
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Table  4- 4.  Preset  Input  Values. 


Variable 


Value 


AK 

ALFA 

ALPHT 

ALPHMC 

AMUIN 

BETAIN 

BETAMC 

BK 

CH 

CM 

CMMC 

CMVR 

DELTIN 

DTHM 

FMAX 

FMIN 

IQUIT 

KX 

LPHI 

NBREF 

NCASE 

NN 

Nl’RINT 

RU2M 

STRFR 


0.08 
1.0 
15*  0.5 
0.5 
24*  0.2 
24*  4.0 
0.325 
0.0 
1.15 
0.075 
0.08 
0.25 
24*  4.0 
0.1 

100000 

50 

50 

15 

9999 

2 

1 

0 

1 

3.0 

0.01 

1.259921 


STRFX 


A. 5.1  Notes  on  Input 

1.  The  jet  nozzle  geometry  is  specified  by  input  of  the  number  of 
component  boundaries,  NEST,  along  with  pairs  of  coordinates,  RA  and 
DAM’,  for  each  boundary  element.  The  ambient  field  is  always 
treated  as  the  first  boundary  in  the  input  arrays  for  UE,  PT,  TT, 
LEAV,  NUM,  KN,  XE,  and  ALPO.  This  is  why  some  numbers  have  already 
been  filled  in  on  Table  4-2  in  the  first  column  for  these  arrays. 

A  nozzle  with  N  elements  has  NEST  =  N  +  1  boundaries. 

2.  The  steps  to  specifying  nozzle  geometry  input  are  as  follows,  re¬ 
ferring  to  Figure  4-1: 

a.  Obtain  sketch  or  drawing  of  nozzle  exit  cross  section  and 
select  a  coordinate  origin  which  is  optimum  from  the  stand¬ 
point  of  symmetry  and  boundary  point  specification. 

b.  Number  each  boundary,  reserving  boundary  Number  1  for  the 
ambient  field. 

c.  With  respect  to  the  coordinate  origin,  select  a  reference 
angular  location  for  each  boundary,  ALP0. 

d.  For  each  boundary,  select  points  represented  by  pairs  of 
coordinates.  The  coordinates  used  as  input  are  radius, 

RA(I,J),  and  angular  increment  from  the  preceding  point, 
DALP(I,J).  For  the  first  point,  DALP(I,J)  is  the  angular 
increment  from  the  reference  angle  ALP0.  The  index  I  is  the 
boundary  point  number,  and  the  index  J  is  the  boundary  number. 
Both  ALP0  and  DALP  are  to  be  input  in  radians,  and  RA  is  input 
in  feet. 

3.  The  last  point  on  a  given  boundary  should  be  located  at  ALP0  if 
the  boundary  has  only  one  leaf.  The  sum  of  all  DALP(I,J)  should 
equal  zero  if  the  boundary  has  only  one  leaf. 

4.  If  the  boundary  is  a  circle  about  the  origin,  only  one  point  on  the 
boundary  need  be  supplied,  and  the  value  of  LEAV  for  that  boundary 
is  set  equal  to  the  number  of  boundary  points  desired  on  the  circle. 

5.  The  program  uses  linear  interpolation  between  input  boundary  points. 
If  a  boundary  is  made  up  of  or  contains  straight  line  segments, 
only  the  end-points  of  the  straight  line  segments  need  be  input. 

6.  The  variables  PT  and  TT  refer  to  stagnation  pressure  and  temperature 
at  the  exit  plane  inside  the  boundary  of  interest.  Setting  the 
first  value  of  PT  equal  to  PS  gives  a  static  ambient  field.  The 
lirst  value  of  PT  greater  than  PS  simulates  non-zero  flight  veloc- 

i  ty. 
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7. 


The  variable  LPHI  determines  what  angular  extent  of  the  flow  field 
needs  to  be  calculated.  If  the  nozzle  geometry  is  axisymmetric , 
setting  I. PHI  equal  to  a  large  number  (such  that  2n/LPHI  is  less 
than  D  HIM)  toivos  t  lie  program  to  calculate  t  lie  I  low  lielil  at  oiilv 
one  angular  location.  the  flow  field  for  a  nozzle  containing  two 
adjacent  circular  jets,  for  example,  has  LPHI  =  4,  since  the  flow 
is  the  same  each  quadrant.  Several  examples  of  how  boundary  para¬ 
meters  are  specified  are  shown  in  Figure  4-5. 

8.  The  program  can  currently  only  handle  coplanar  nozzles;  that  is, 
every  nozzle  element  must  terminate  at  the  same  axial  location. 
Therefore  XE  must  be  the  same  for  all  input  boundaries. 

9.  The  centerbody,  if  any,  is  input  through  coordinates  pairs  XCBDY(J), 
RCBDY(J),  where  1  J  <_  NCBDY.  A  maximum  of  40  points  can  be 
input.  The  LSPFIT  subroutine  uses  this  input  to  interpolate  for 
finding  the  values  of  RMIN  at  each  axial  location  X.  The  LSPFIT 
routine  can  treat  line  segments,  both  straight  and  curved.  Typical 
examples  of  centerbody  coordinate  input  are  shown  in  Figure  4-6. 

If  there  is  no  centerbody,  the  user  can  avoid  automatic  computation 
of  the  potential  core  of  axisymmetric  nozzles  (which  has  no  impact 
on  mixing  noise)  by  specifying  RMIN  as  input,  but  with  NCBDY  =  0. 
This  option  causes  the  computation  to  begin  at  r  =  RMIN(KA),  where 
KA  is  the  axial  station  number. 

10.  The  input  value  of  CM  is  modified  for  velocity  ratio  and  Mach 
number  effects  by  the  relation 


(1  +  CMVR*VR) (1  +  CMMC*ACH) 

where  DBDX  is  the  modified  value  of  Cm,  and  VR  and  ACH  are  the 
velocity  ratio  and  Mach  number,  respectively,  of  a  given  boundary. 
The  heat  transport  spreading  parameter  is  then  calculated  from  the 
relation 


Ctl  =  CH  *  DBDX 

The  values  of  CM,  CMMC,  CMVR  and  CH  recommended  and  preset  in  the 
program  are  given  in  Table  4-4.  These  values  can  be  changed  by  the 
user  to  reflect  experimental  evidence  if  so  desired. 

11.  The  axial  locations  of  the  axial  stations  can  be  input  by  the  array 
X(KA),  where  1  £  KA  ;<  KX.  The  radial  mesh  step  size  can  also  be 
input  by  the  array  DSIG(KA).  An  automatic  grid  selection  procedure 
has  been  devised  to  obviate  the  need  for  supplying  all  values  of 
X(KA)  and  DSIG(KA).  The  only  input  required  is  the  first  axial 
station  X ( 1 ) ,  and  the  grid  stretching  factors  STRFR  and  STRFX.  The 
grid  is  then  calculated  from  the  following  relations: 
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(a)  Circular  Jet 


J  -  2 
J  =  3  - 


(b)  Coannular  Jet 


c-J  =  2 

f 

k  . 

if  I  -  l 

9 

(c)  Rectangular  Jet 


NEST 

=  2 

LPHI 

=  999 

ISYM 

1 

LEAV 

1,36 

NUM 

1.1, 

ALP0 

=  0,0, 

KN  = 

1,1. 

NEST  =  3 
LPHI  =  999 
ISYM  =  1 
LEAV  =  1,36,36, 
NUM  =  1,1,1, 
ALP0  =  0,0,0, 

KN  =  1,1,2, 


NEST  =  2 
LPHI  =  4 
ISYM  0 
LEAV  -  1,2, 
NUM  =■  1,3, 
ALP0  =0,0, 
KN  =  1,1, 


Figure  4-5.  Examples  of  How  Boundary  Parameters  are  Specified. 


Curved  Line  Segment 


NCBDY  =  7, 

XCBDY  -  X1,X2>X3>X4lX5lX6(-X5) ,X7, 
RCBDY  =  Rj ,R2.r3.r4)°>0.0. 


(a)  Example  1  -  Curved  Centerbody 


NCBDY  -  10, 

XCBI)\  -  Xj.X2.X3  (=X2)  ,  X4  ,  Xj(=X4  )  ,  Xg  ,  Xy  ,  Xg  ,  Xg(=Xg)  ,  X  jo  > 
RCBDY  -  Rj , R2 . R3 (  R2? »R4 # R5 (=R<P  > R6 » R7 ’ Rg > Rg(=R8) / R10 1 


(b)  Example  2  -  Segmented-Cone  Centerbody  with  Curved  Tip 


Figure  4-6.  Centerbody  Input  Coordinate  Examples. 
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X(KA)  =  STRFX  *  X(KA-l) 


DSIG(KA)  =  STRFR  *  X(KA) 

This  provides  a  grid  which  exhibits  larger  and  larger  step  sizes  as 
the  plume  is  developed  downstream.  Recommended  value  of  STRFR  and 
STRFX  are  preset  and  listed  in  Table  4-4. 

The  variable  NN  determines  the  type  of  acoustic  calculation  desired. 
Normal  (preset)  operation  is  with  NN  =  0,  which  give  the  complete 
acoustic  calculation.  The  user  may  desire  to  perform  diag¬ 
nostic  computations  to  assess  the  relative  importance  of  con¬ 
vection,  shielding,  etc.  By  selecting  the  appropriate  value 
of  NN,  the  various  components  of  the  acoustic  calculation  can 
be  switched  on  and  off  in  various  combinations.  Setting  NN  =  4 
gives  only  the  aerodynamic  calculation,  and  the  acoustic  calcula¬ 
tions  are  bypassed.  The  various  options  for  NN  are  listed  below: 

NN  =  0  -  complete  acoustic  calculation. 

NN  =  l  -  convective  amplification,  no  shielding. 

NN  =  2  -  no  convective  amplification,  no  shielding. 

NN  =  3  -  no  convective  amplification,  with  shielding. 

NN  =  4  -  no  acoustic  calculation,  aerodynamics  only. 

The  printout  of  aerodynamic  flow  field  data  is  controlled  by  NPRINT. 
When  NPRINT  =  0,  no  aerodyanmic  printout  is  provided.  If  NPRINT  = 

1,  aerodynamic  printout  is  provided  at  every  axial  station.  If 
NPRINT  =  2,  aerodynamic  printout  is  provided  at  every  second  axial 
station  (i.e.,  KA  =  1,  3,  5,  7,  etc).  For  PRINT  =  3,  printout  is 
provided  at  every  third  station,  etc. 

For  dual  flow  nozzles,  if  the  inner  stream  has  a  higher  velocity 
than  the  outer  stream,  use  ALPHMC  *  0.5  and  BETAMC  =  0.325  (preset 
values).  These  variables  are  weighting  factors  in  the  convection 
Mach  number  calculation,  which  is  computed  from  the  relation 

MC  =  ALPHMC  *  MACH  +  BETAMC  *  EMACH 

where  MACH  is  the  local  acoustic  Mach  number  U/Ca  and  EMACH  is  the 
exit  plane  reference  Mach  number  Uj/Ca.  If  the  outer  stream  has  a 
higher  velocity  than  the  inner  stream,  use  ALPHMC  =0.5  and  BETAMC 
=  0.325/VR,  where  VR  =  (Uouter/^inner) j •  For  multielement  suppres¬ 
sor  nozzles,  VR  =  Uj/Um,  where  Um  is  the  postmerged  potential  core 
velocity.  If  Um  is  not  known,  BETAMC  =  0.2  to  0.25  is  usually  a 
good  approximation. 
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lor  dual  flow  nozzles,  input  BETA IN  =  4.0  (preset)  for  all  values 
ol  X,  provided  the  inner  stream  velocity  is  higher  than  the  outer 
stream  velocity  at  the  exit  plane.  If  the  outer  stream  velocity  is 
higher  than  the  inner  stream  velocity  at  the  exit  plane,  input 
RETAIN  =  0  for  all  axial  stations  where  X  10*DEQ  (NBREF) ,  and 
BETA  IN  =  4.Q  thereafter,  where  NBREF  is  the  outer  stream  boundary 
number.  For  multielement  nozzles,  input  BETAIN  =  0  for  axial 
distances  less  than  10*DEQ(1),  where  DEQ(l)  is  the  equivalent 
diameter  based  on  total  flow  area  at  the  exit  plane. 

lb.  For  each  boundary  element  DEQ,  DS  and  NCELL  are  input.  The  first 

value  of  DEQ  is  the  total  flow  area  equivalent  diameter.  The  first 
value  of  NCELL  determines  whether  or  not  the  shock  cell  noise  is 
computed.  If  NCELL (1)  is  input  zero,  no  shock  noise  is  computed; 
tor  NCELL(l)  0,  the  shock  cell  noise  routine  is  called.  The  shock 
noise  of  each  boundary  element  is  computed  separately  and  added  to 
the  total  noise.  If  any  boundary  has  a  value  of  NCELL  =  0,  that 
boundary  element  is  bypassed  in  the  shock  noise  calculation.  It  is 
recommended  that  NCELL  =  8  be  used  for  each  element  unless  the  actual 
number  is  known. 

4.5.1  Examp le  Case  Input  Select  ion 

To  illustrate  how  geometric  input  parameters  are  selected  for  a  complex 
nozzle  geometry,  an  example  is  presented,  taken  from  Reference  6.  The 
example  nozzle  exit  geometry  is  shown  in  Figure  4-7.  Consideration  of  this 
figure  indicates  that  information  over  a  45°  sector  of  the  flow  field  will  be 
sufficient  to  describe  the  complete  flow  field.  This  is  one-eighth  of  a 
circle,  thus  LI’HI  =  8.  Neither  axial  total  similarity  or  dissimilarity 
exists  so  I SYM  is  0.  Counting  the  number  of  closed  contours  indicates  a 
value  of  NEST  of  12,  where  one  is  included  for  the  ambient  or  external 
field.  Values  of  FT  and  TT  must  be  provided  for  the  exit  state  existing  just 
within  each  of  these  contours.  Values  of  XE,  ALPO,  LEAV,  NUM,  KN,  DEQ,  DS , 
and  NCELL  must  be  provided  for  all  the  contours  except  the  first  which  is  the 
boundary  at  infinity.  Values  of  these  parameters  for  the  contours  shown  in 
Figure  A- 3  are  now  considered  in  the  following  discussion. 

Boundary  2:  Description  of  this  boundary  starting  at  45°  to  the  system 
baseline  is  convenient.  Thus  ALPO  =  n/4  radians.  Since  each  90°  sector  of 
tlu‘  contour  is  identical  with  the  proceeding  one,  LEAV  =  4.  Since  the  pro¬ 
gram  assumes  straight  lines  to  exist  between  successive  boundary  points, 
description  of  this  boundary  is  possible  with  only  three  points  for  each 
quadrant.  These  are  P(l,2),  P(2,2),  and  P(3,2).  Each  point  is  described  bv  (I) 
its  distance  from  the  system  origin  and  (2)  the  angle  between  (a)  the  line 
joining,  it  with  the  origin  and  (b)  the  line  joining  the  preceding  point  with 
the  origin.  Note  that  no  value  of  RA  is  given  for  the  point  P(0,2)  since  it 
will  be  identical  to  RA(3,2).  The  value  of  NUM  for  boundary  2  will  therefore 
be  1. 


Boundary  3:  This  contour  has  eight  symmetric  leaves;  thus  LEAV  =  8. 

ALPt)  of  0.0  is  as  convenient  as  any  other  value.  The  eight  points  indicated, 
I’d, 3)  through  P  ( 8 , 3 )  ,  probably  are  sufficient  to  describe  the  boundary.  Tims 
NUM  8. 
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Boundary  4:  Since  this  is  a  circle  about  the  origin,  it  can  be  divided 
into  a  convenient  number  of  leaves  and  only  one  point  need  be  given  for  each 
(NUM  =1).  If  a  hundred  boundary  points  are  desired,  set  LEAV  =  100,  DALP(1,4) 
tt/50  and  RA(1,4)  equal  to  the  circle  radius. 

Boundary  5  through  12:  Each  of  these  contours  must  be  described  indi¬ 
vidually  unless  certain  artifical  changes  are  made  in  the  arrangement.  A 
partial  representation  of  Boundary  9  is  shown  in  Figure  4-7.  Note  that  suc¬ 
cessive  points  on  the  boundary  are  obtained  by  progressing  around  the  bound¬ 
ary  in  a  counter-clockwise  fashion.  In  order  to  reduce  the  labor  of  repre¬ 
senting  each  circle  separately,  a  straight  line  can  be  drawn  connecting  each 
circle.  Two  contours  can  then  be  visualized,  one  consisting  of  the  outer 
halves  of  the  circles  and  the  lines,  the  other  consisting  of  the  inner  halves 
of  the  circles  and  the  lines.  Each  contour  has  eight  leaves  and  only  one 
need  be  represented  by  the  programmer.  Since  this  technique  requires  the 
computer  to  integrate  along  each  straight  line  twice  in  the  course  of  computa¬ 
tion,  it  will  definitely  increase  the  computational  time  over  the  method  in 
which  each  contour  is  represented  separately. 
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4.6  OUTPUT  DESCRIPTION 


The  output  format  is  generally  self-explanatory.  The  input  data  are 
first  printed  out,  using  the  same  nomenclature  previously  defined  in  Table  4-1. 
Nozzle  exit  plane  flow  conditions  (static  temperature,  velocity,  Mach  number, 
momentum  flux,  and  enthalpy  flux)  are  then  printed  out  for  each  boundary 
contour. 

At  each  axial  location  specified,  the  radial  and  tangential  distribu¬ 
tions  of  flow  field  properties  are  printed  out.  After  the  flow  field 
information,  the  noise  characteristics  of  that  particular  axial  station  are 
then  listed. 

Following  all  of  the  axial  station  flow  field  data,  a  summary  table  of 

the  noise  characteristics  (SPL  spectra,  PNL,  PWL,  OASPL)  is  given. 

Section  4.7  contains  an  input  deck  card  listing  and  output  printout  for 

a  sample  case  run.  This  particular  case  is  for  a  7-tube  nozzle  presented  in 

Reference  2.  For  brevity,  only  a  portion  of  the  total  output  is  shown;  but 
the  formats  of  the  various  output  data  are  all  included. 

Two  warning  flags  are  built  into  the  program.  The  first  is  a  case 
termination  flag,  which  occurs  whenever  an  input  total  pressure  (Pj)  is  less 
than  the  input  static  pressure  (Pg) .  The  flag  message  is  as  follows: 


****ERROR  -  MACH  NO.  SQUARE  IS  NOT  GREATER  THAN  ZERO  -  CASE  WILL 
TERMINATE**** 


The  second  flag  is  a  warning  detected  in  subroutine  SLICE,  which  occurs 
whenever  the  number  of  turning  points  (NTP)  is  found  to  be  greater  than  2. 
The  f Lag  message  is  as  follows: 


WARNING  -  NO.  OF  TURNING  POINTS  IS  CREATER  THAN  2  AT 
KA  =  X  =  ,  ITH  =  _ ,  THETA  = 


,  NTP  = 


where  KA  is  the  axial  station  number,  X  is  the  axial  location,  ITH  is  the 
observer  angle  index,  THETA  is  the  observer  angle  in  degrees  (6t),  and  NTP  is 
the  number  of  turning  points  found.  The  two  outermost  turning  points  are 
used  and  those  inboard  of  these  two  are  discarded  in  such  cases,  since  the 
acoustic  shielding  model  can  only  accommodate  up  to  2  turning  points.  The 
noise  output  at  those  values  of  Bj  where  this  warning  appears  should  be 
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treated  as  suspect,  since  the  acoustic  shielding  effects  are  not  properly 
modeled.  This  is  most  likely  to  occur  in  the  initial  mixing  regions  of 
multitube  nozzles,  where  multiple  peaks  in  the  azimuthally  averaged  velocity 
profiles  are  likely  to  occur. 
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4 . 7  SAMPLE  OUTPUT  LISTING 

An  example  case  of  a  7-tube  multielement  nozzle  is  described  here, 
selected  from  one  of  the  data/theory  comparison  cases  presented  in  Reference 
2.  The  nozzle  consists  of  a  hexagonal  array  of  0.875-inch-diameter  tubes, 
with  a  spacing/diameter  ratio  of  3.  The  acoustic  arena  is  a  9-ft-radius 
arc.  The  geometry  is  illustrated  in  the  sketch  below. 


Observer 


(I  0.875  in. 
1)  8.0  in. 


The  input  data  cards  for  this  case  are  listed  in  Table  4-5.  Note  that  all 
geometry  input  lengths  are  in  feet,  and  all  input  geometry  angles  are  in 
radians.  The  output  listing  for  this  case  follows  Table  4-5. 


136 


Table  4-5.  Input  Data  Card  Listing  Sample  Case 


5R329  01  10-06-77  16.471  ***  INPJT  DATA  CARD  LISTING  --  H*G«8  *** 

CHD  7-TUBE  Afi*2.S  Nn7ZLE  -  VJ*??00  EPS  -  TTjmfcOO  OEG-R 

$ 1NPU  T 

NEST=8,  LPHI=12,  I S  T  M  =0  /  IQU1T=100, 

RU2N=3,  0  T HM  =  0.1/  PS=2116, 

ATOT AL  =  0. 029231 ,  D E3 = 8 *0. 0 7  2 9 1 6  7,  DS*8*U. 0729167,  NCELL*8*3, 

KN  =  S  *  1  ,  X  E  =  8  *  0  , 

GAM=1.S5,  CP=6619, 

ALP0=0. 0,5.96144,0.725447,1  .772 64/2. 8198/3. 86 70/4. 9142/0.0/ 

LEA  9  =  0/  6*1/24/  iM  UM  =1,6*24/1,  KN*8*1,  XE=8*0, 

DALP ( 1/2) = 

.033596, .045590, .054084, .059862/. 063450,. 065168, 

.06 51 68/. 063450/. 059 862/. 054084/. 045590/. 033596/ 

.01 70 39,-. 0053 17, -.034 336/-. 066 98 4,-. 10 35 91, -.12656 2/ 
-.126562/-.103591 /-.068984, -.034336/-. 00531 7, .01 7039/ 

D  ALP ( 1 / 3)  = 

.0335 96/. 045590/. 0  54 084,. 059862/. 0  63 450/. 0651  68/ 

.0651 63/. 063450/. 059862/. 05 40 84, .045 590/. 033596/ 

. 0170 39,-. 00 5317/-. 034336/-. 068 98 4, -.103591,-. 126562/ 

-.126 56 2/-. 103 591 ,-.068984, -.034536, -.005 31 7/. 01 70 3 9/ 

DALPt 1/4) = 

.053596/. 045590/. 054084,. 059862/. 063450, .065168/ 

.065163/. 063450/. 059862/. 054084  /. 045590/. 033596, 

.017039, -.005317, -.034536, -.068 98 4, -.103591, -.126562/ 

-.126 562, -.103591, -.068984, -.034336/-. 005 31 7,. 01 703 9, 

DALPt 1/5 ) = 

. 0335 96 /.C45590,. 054 084, .059 862/. 063450,. 06 51 68, 

.06 51 63, .063450, .059862/. 054084, .0455 90 ,.033596, 

.01  70  39, -.03531  7, -.034356/-. 0  68 93 4, -.  133591  , -.  1  2  65  62/ 
-.1265S2/-.103591 ,-.068984, -.034336/-. 00531 7, .01 7039, 

DALPt 1,6) * 

.0335 96, .045590/. 054 084, .059 862/. 063450/. 06 5168/ 

.06 5  1  68/.  063450/.  0  59 862,. 054084  , .  045590/.  03359b, 

.017039, -.005317, -.034556, -.068 93 4, -.103591, -.126562/ 

-. 126562/-. 103591, -.  3689 84,-. 034536/-. 005317, .01  73 39, 

D  ALP  t 1 , 7  )  = 

.033596,. 045590, .054084, .059862/. 063450/. 06 51 68/ 

. 06 51 63 ,.06 54 50/. 0 59 862 ,.05 4084, .045 5 90 ,.03 35 96, 

.01  70 39, -.0053 1  7, -.0  34536/-. 068984, -.  103591  , -.  1  2 656 2/ 

-. 126562,-. 1 0 3591, -. 06898 4, -.034336,-. 005 317,. 01 703 9, 

R  A  f  1  ,  2) * 

.12392,.13145/.15759,.14212/.144S9,.14583/ 

.1448 9,. 14212, .137  59,. 13145/. 12  39 2, .11529, 

.1059 6/. 39646/. 03 748/. 07991,. 07476/. 07292, 

. 0747 6/. 37991, .08743/. 09646, .10596/. 11529, 

R  A  1 1  , 3)  = 

.12392,.13145,.13759,.14212,.14489,.14583, 

.1448 9,. 14212, .13759, .13145, .12 39 2, .11529, 

.  1059 6,. 39646, .03748, .0  79 91  , .07476, .07292, 

. 0747 6/. 37991, .03748,. 09646, .10596, .11529, 

RAt 1 ,4) = 

.12392, .13145/. 137 59, .14212,. 14 48 9, .14583, 

.14489,. 14212, .137 59, .13145, .12  39 2, .11529, 

.  10596/. 09646, .087 43, .07991  , .07476, .07292, 
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Table  4-5.  Input  Data  Card  Listing  Sample  Case  (Concluded) 


5  R  52  9  01  10-06-7?  16.471  *•»  INPUT  DATA  CARD  LISTING  --  M»G*B  *** 

.074 76,. 07991 ,.08 748,. 09646,. 10 59 6,. 1 1529, 

RA< 1 , 5) = 

.12392/.1 51 4 5,.1 5 7 59,. 142 12,. 14489,. 14 585, 

.1448  9,  .14212, .15759, .13145, .12  39 2, .11529, 

.  1059  6,  .09646, .08748, .07991  , .07476, .37292, 

.074 7 6,. 07991, .08 743, .09646, .10596, .11529, 

RAC  1 ,6)  = 

.  1239 2,. 1  3145, .137 59, .  14212, .14 48 9, .  1  4585, 

.14489,. 14212, .137 59, .13145, .12 392, .11529, 

.  1059  6,  .09646, .087 48, .07991  , .07476, .07292, 

.07476, .07991, .08748, .09646, .10596, .1  1529, 

RA  ( 1  ,  7)  = 

.  1239 2,. 1  3145, .137 59, .14212, .14 489, .14583, 

. 14489,. 1421 2,. 13 759,. 131 45,. 12392,. 11 529, 

.  1059 6,.  09646,. 08  748,. 07991  , .07476, .  07292, 

. 0747 6,. 07991, .08748,. 09643, .10 596, .11529, 

DAL  PCI, 8) =0.2618,  RA (1,81=0.0  36 458, 

ALPHMC=0.S,  BFTAMC=0.25, 

fMlN=10C,  F 1  A  *  =80  000 ,  NUMANG  =  1,  DIST=9.0, 

K  X  =  2  4  »  X=0. 0729167,  STRFR  =  0.01, 

OSIG=10*0. 03729167,  14«0,  N09 = 1 0 • 20 , 1 4 *3 , 

BF  TA  IN=1 5  *0.C, 

NPR I N  T  = 

NCA$£=1 , 

PI=21 16,7*573?,  TT«SaO,7»1feOS, 

S 
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COMPUTATION  OF  ACRG-ACOUsTU  P*-*'  PERI  XES  0*  Sijp*'*'*  SaC~.  NC»  JELLS  - _ 

case  no.  l  cm  7- r»#*-*er  &°-2,i  v*?JLf  -  vj=2?:.o  e^s  -  ttj=i*oo  ofg-p 


INPUT  0  AT  A 


K  X=  ?U 

NEST  = 

8 

LPhI=  12 

1 SY  M  -  J 

NPPINT*  6 

CM=  ,075 

CH=  1.1SC 

r,AM*  \ 

.35** 

CP=  6M9.0 

CTHMs  .KOO 

Ru2M  =  3.0COO 

PS=  2116.0 

COMPUTATION  PFSH  CONTROL  parameters............ . .  TURBULENCE  CONSTANTS _ 


SLICE  NO. 

X 

osio 

-  PM  IN 

NOV 

BETA 

-DELIA 

_ MU- - 

_ 

.-?7  3*2 

Vor.7?Q 

e.coooo 

20 

0.00 

aToo' 

.20 

2 

.  *>6  l  h  7 

.03 729 

o.oooo:- 

2u 

0.00 

4 . 00 

--  .20  . 

1 

. 11575 

.  V-729 

O.OOOOy 

2) 

0.00 

4.00 

.20 

4 

.1*543 

.3; 724 

0.00031 

?) 

0.00 

4.00 

.20 

S 

.1*374 

.  r  r  7  ?  q 

0  .  ft  0  0  C  n 

2A 

o.oo 

4.00 

.20 

ts 

.231  «■<! 

,rr  ?:>4 

(*.0  0  POO 

20 

__O.OP 

-  _f*.S0 

_ *?P__ 

7 

.  1  7 

,  =V  7P^ 

0 • 0  QGO  0 

2.3 

0.00 

4.00 

.20 

H 

.  3*740 

. 

e.oocoo 

2*3 

o.oc 

4.00 

.20 

.4*  »o  *» 

.0  7P.J 

0. *0000 

2r 

n.co 

4.00 

.20 

)  J 

.C.M3  ^3 

.0  739 

c.coooc 

20 

0.00 

4.00 

.20 

1  1 

.  7  3*.^S 

.O'”  7  3< 

O.iOOCC 

(• 

0.00 

<*.00 

.20 

12 

.0 

p. 

.O.PO 

_ 4,00 

..  _  •  20_ 

»  3 

>  .  1***7 

.  >11'/ 

-J.  ,V'f  > 

.' 

o.co 

4.00 

.20 

1  *• 

1  .<**3  41 

.  1-7, 

o. :ooo  • 

0 

0.00 

4.00 

.20 

:  8 

l  .  -  •  1  -  ? 

r  .  '  n !  ft  -• 

0 

o.oo 

4. DO 

.20 

?.  run 

.  r  2  .'  *»  1 

o.  :c,i ): 

0 

a. oo 

4.00 

.20 

17 

,?,4  Tv.jp 

•  3  24* .. 

>) .  •.  ni  0 0 

) 

A. 30 

4,00 

.20 

ta  _  . 

_ 3. 7 J  193  _ 

_  _  ..OOJOn 

_  L.OOOOJ _ 

_ 0 

_ 0.00 

_ .4,00. 

.^20. 

19 

4.  *>6667 

.04**7 

o.coooo 

0 

4.00 

4.00 

.20 

2Q 

S. 87463 

.058*0 

0.00000 

0 

.  .  4.00 

4.00 

-.20 

21 

7.a;7H7 

.0740M 

0.00000 

0 

4.00 

4.00 

.20 

2? 

<1.33333 

.09313 

0.00000 

.  G 

_  _  4.00 

_ 4.00 

_ .20. 

2  3 

11.75426 

.11769 

0.00000 

0 

4.00 

4.00 

.20 

-  24— _ 

—  .14*81574. 

.A4616. 

__Q.QQQ.00 

Q_ 

4.00 

_ 4.0.0, 

_ =20_ 

, 

— 

- XEt.2>  = _ 0.Q0  ALPOl 

21=  S. 

9614 

LEAVC  21 

=  -  1 

NUMJ. 

2L* 

24_  _ 

K.ML  2) 

=  i 

— 

.  _  D ALP  l  1.2)  = 

.3336 

RAl 

l*  2)  ~ 

.1239 

DALP  ( 

2» 

2)  5 

.0456  _ 

__RA  L_2j»  27  = 

_ *  1316 

DAI  P  1  3.  _21  = 

PAL 

2,  2)~ _ 

.  1_376_  DALPL 

<M 

21  = 

.  0599 

RAC  4,  21  = 

.1421 

> 

_ 

OALp(  S.  2) = 

.0635 

PA  l 

5.  21  = 

.1449 

DALPl 

6  • 

21  = 

.0652.  ... 

RAl  6.  2) = 

.1456 _ 

1 

—  - 

_  DALPl  7 «  2)  = 

.0652 

RAC 

7.  2)  = 

.1449 

DALPC 

a. 

21  = 

•  0635  _ 

2).«. 

.1421 _ 

Y 

OAlP(  9«  21  = 

.1599 

PA  ( 

9  «  21  = 

.1376 

OALp(10» 

21  = 

•  054  1 

PA110.  21= _ aJUS 

± 


\ 

_  -  _  DALPl  21*  2>  = 

.3456 

RAl  11, 

21  = 

.1239 

DALP 1 ] 2 1  2)  = 

.0336 

RAU2.  2)x  - 

.1153 - 

t 

_  DALP  11 Ji  21  = 

.1170 

_RAU3» 

2)  = 

.1060 

DALP  (14 »  2)  - 

-.0053 

_ RA114-L  21  =_ 

.0965 

DALP  115*..  2)  = 

_-.1343. 

RAC1Sa_21= 

.HP_7S 

['ALRUAi  Zl- 

-.0690 

RACI6,  21  = 

.0799 

J 

DALP  1 1 7 •  21  = 

-.1036 

»A11 7. 

21  = 

.0743 

DALP 1 1 8 »  21  = 

-.1266 

-.-RAILS.  2)  = 

.0729  _ 

OALp ( 1 °#  21= 

-.1266 

R  M  1  °  * 

2)  - 

.0740 

PALP C20 *  21  = 

-. 1C  36 

_  °M20,  21  = 

.0799 

( 

. . . .  DALP  (21  •  21  = 

-.at?: .. 

. .<-  '  (uli 

21  = 

_..C675 

f'ALp122i.  2)  = 

-.034  3. 

_ b4J22,_2.>.=_ 

,.0965- 

DALP (2  3  •  21  = 

-.0053 

c  A  (  2  3  * 

21  = 

.1?60 

0 ALP  ( 24  *  21  = 

.  0  1  7  Q 

RA  <24.  21  = 

.1153 

1 
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COMPUTATION  OF  AERO-ACOUSTIC  PROPERTIES  OF  SUPPPESSOR  NOZZLES 
CASE  NO,  l  CRO  7-TUHE  AR*?,3  N077LE  -  VJ=??00  FPS  -  TTJ=1600  OEP-R 


XE  ( 

3>  =  0.00  ALPO(  3)  = 

,7?S4  LEAV<  3) 

=  1 

NUM (  3>= 

?4 

*N(  3 )  = 

1 

DALP<  1*  3>  = 

.0336 

86(  1.  31  = 

.1P39 

DALP (  ?. 

3)  = 

.0456 

PA  (  ?* 

31  = 

.1315 

DALP<  3*  3 ) = 

•  0S4  l 

RA(  3*  3) = 

.  1  J76 

DALP  l  4. 

31  = 

.0599 

R A  (  4. 

31  = 

.  1  6?  1 

OAtP(  5.  3)  = 

.0635 

RA  <  5*  3) = 

.1449 

DALP (  4, 

3)  = 

.065? 

RA  t  6, 

3)  = 

.1658 

OALP<  7*  3>  = 

•  06S? 

PA)  7.  3) = 

.1449 

DALP  J  8. 

3)  = 

•  C635 

RA<  H. 

31  = 

.  1 6?  1 

OALP(  9,  31= 

.0599 

HAl  V.  3) = 

.1374 

DAt  P  (  |  0. 

3)* 

.0541 

RA  (  10* 

3)  = 

.1315 

DALPUl*  3>  = 

•  0656 

RAUlf  31* 

.1239 

DALP ( !?• 

3)  = 

.03  36 

RA ( I?. 

31  = 

.1153 

OAt  P  M  3  ♦  3)  = 

.0)  70 

RA(13.  3) = 

.1040 

DA(  p<  14, 

3 )  = 

-.005  1 

RA  <14. 

3>  = 

.0965 

DALP <15,  3>  = 

-.0343 

PA (IS*  3)3 

,08  75 

DALP  U  4* 

T>  = 

-.0490 

PA (14. 

3)  = 

.0799 

OAtPUT,  3)  = 

-.103^ 

RA (1 7 ,  3)  = 

.0748 

DAI.P  (  l  8, 

3)3 

1?46 

PA  < 1 H. 

3 )  = 

•  0  7?9 

DALp  MR.  3 )  = 

1?66 

RA  < 1 9 •  3)  = 

.0  7  4ft 

flAl.H  l?0. 

3 )  = 

-.1034 

RA  <?0« 

3)  = 

.0799 

DAlp ( Z l .  3) = 

-.0690 

RA (?  1  •  3)  = 

.0875 

0*1 P<??. 

3)  = 

-.0343 

W A ( ?? ♦ 

3)3 

,0965 

OAtP (23*  3  >  = 

-.0053 

PA(?T,  3>  = 

.1040 

DALP (?4 • 

3)  = 

•  0  l  70 

RA<<>4, 

31  = 

.1153 

XE< 

4 1  =  0-00  ALP0(  4>  =  1 

, 77?6  L£AV(  4) 

=  t 

NUM (  4 \ - 

?4 

KN  (  4)  = 

1 

OALP(  1*  4) a 

.0336 

P  A (  1*  4)s 

.1?39 

DAIP(  ?, 

4)  3 

•  0466 

RA<  ?. 

4)  r 

.1315 

0Al.P<  3.  4)s 

•  CS41 

RA(  3.  4) = 

.  11 7e> 

DALP<  4. 

4)3 

.0599 

RA  (  4. 

4)  = 

»  1  6  ?  1 

DA^M  1  5.  61  = 

,0635 

PA  (  S#  4>  = 

•  1449 

0ALP(  6. 

4)3 

.  046? 

RA  (  6, 

4)  = 

.1658 

OALP(  7 .  4>= 

.065? 

RA<  7.  4 ) = 

.1449 

0 ALP (  8, 

4  )  s 

.0635 

RA<  ft. 

41  = 

•  1 6  ?  1 

PALP (  9,  4)3 

.0599 

M A  <  9.  4)3 

.1374 

DALP  no. 

4)  = 

.0641 

RA  (10. 

4)  = 

.1315 

OALPUlt  4)3 

,  3656 

PA  (U  •  4>  = 

.  1?39 

0  At.  P  <1  <?  , 

4)3 

.0  336 

RAM?, 

4)  = 

.1153 

DALP(I3,  4)= 

.01  30 

94(13.  61  = 

•  1  «>60 

f’ ALP  (  »4f 

4)3 

-• 0o5  3 

RA ( 14, 

4)  = 

•  0^65 

PAi  PUS.  4)3 

-.0363 

RAUS.  4)  = 

.0875 

DALP (14, 

4)  = 

-.0690 

RA  (16. 

4)  * 

.0799 

PALP  M  7  *  4)3 

-.1036 

PA ( 1 7  *  4)  = 

•  0  748 

DALP (18, 

4)3 

-. 1?66 

RAUft, 

4)  = 

•  07?9 

OALPUP*  41  = 

-,l?66 

R&  <19*  4>s 

.0748 

OAl P  <  ?0  « 

4)  = 

-.1036 

RA  <  ?o ♦ 

4)  = 

.0799 

PALP < 21 •  4 ) = 

-.0690 

PA(?1.  4)= 

•  0«75 

PALP ( ?? • 

4)  = 

-.0343 

RA ( ??  * 

4  )  = 

.0965 

P AL  P ( ?  1 1  4)  = 

-.0051 

RA  <?3,  4)3 

.  1  040 

DALP (?4  , 

4)  = 

.01  70 

PA  <?4, 

4)  = 

.1163 
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COMPUTATION  or  aero-acoustic  properties  oe  suppressor  nozzles 
CASE  no,  l  CRO  7-TUHe  AR=2.3  NOZZLE  -  VJ=??00  EPS  -  TTJ=16C0  PEP-P 


XE( 

S»a 

0.00 

At  PO (  51  =  2. 

ft  1 9 A  LEAV<  51 

=  1 

NlJM  (  51* 

24 

KN(  51= 

1 

OALP(  I • 

5)  = 

•  0336 

RA  (  1  • 

51  = 

.123^ 

PALP (  2, 

51  = 

.0456 

RAC  ?, 

51  = 

.1315 

DALP<  3* 

51* 

.os^ii 

R A  (  3. 

51  = 

.  1  376 

PALPI  4. 

5l  = 

.0599 

RA(  4, 

51  = 

.1421 

OALP(  5* 

S )  = 

.0635 

RA  <  S. 

51  = 

.1449 

PALPI  6, 

51  = 

.0652 

RA  (  6, 

51  = 

.1458 

OAL  P <  7* 

51  = 

.065? 

RA  (  7, 

51  = 

.14^.9 

PALM  *• 

51  = 

.0635 

RA  (  H, 

51  = 

.1421 

DAL  P (  9  * 

51  = 

.0599 

PA  (  0. 

5)  = 

.  1376 

PALP (10, 

51  = 

,0541 

RA ( 10, 

51  = 

.1315 

DALP( 1 1 t 

S )  = 

.0456 

RAC  11  ♦ 

51  = 

.1239 

OALP ( 12, 

51  = 

.0336 

RA ( 12, 

51  = 

.1153 

PALP (13* 

51  = 

.0170 

R  A  ( 1  3  , 

51  = 

.1060 

PALP ( 14. 

51  = 

-.0051 

RA( 14, 

5)  = 

.0965 

OALP( IS. 

5)  = 

-.0343 

RA ( 1 5 ♦ 

51  = 

.0*75 

PALP( 16, 

51  = 

-.0690 

R  A  ( 1  6  , 

5)  = 

.0799 

OALP  <17* 

51  = 

-.1036 

RA  (  1  7  , 

51  = 

.074* 

nfti  P (  1 H. 

5)  = 

-. 1 266 

RA ( 1*. 

51  = 

.0729 

PALP (19, 

5  >  * 

-.1266 

RA ( 1 9. 

51  = 

.0/4* 

PM.  P  (  2 0  , 

M  = 

1  0 

RA ( 20 , 

51  = 

.0799 

n«LP<?l. 

51  = 

-.0690 

RA (2 1 , 

51  = 

.0*75 

PALP (22, 

51  = 

-.0  14  3 

RA  (22, 

51  = 

.0965 

Pai P(23. 

51  = 

-•DOS'S 

R  A  ( ?  3  , 

51  = 

.1060 

PALP (24 , 

5)  = 

.0170 

RA (24, 

51  = 

.  1  153 

XE( 

61  * 

0.00 

ALPO<  61=  3, 

8670  LC AV (  6) 

=  l 

NlJM  (  M  = 

Pu 

*N(  61= 

1 

OALP I  I « 

61  = 

.0336 

R  A  (  1, 

61  = 

.1239 

OALP<  2, 

61  = 

.0456 

R A  (  2, 

61  = 

.1315 

OALP(  3. 

61  = 

.0541 

RA(  3, 

6}  = 

.1376 

PAI.P(  4, 

6)  = 

.0599 

R  a  <  <.  • 

61  = 

•  1  R?  1 

PAL P (  5. 

6)  = 

.0635 

RA  (  5, 

61  = 

.1449 

(»A|  P<  6, 

M  = 

.0652 

RB  (  f-. 

6)  = 

.1458 

OAlp <  7. 

61  = 

.0652 

RA  (  7, 

61  = 

.1449 

0A|  P(  H, 

6)  = 

.0635 

R A  (  8* 

61  = 

.  i<i?i 

DALP<  9, 

6)  = 

.0599 

R  A  (  9, 

6)  = 

.1176 

PAl.P  (10, 

61  = 

.0541 

R  A  ( l  0 , 

61  = 

.1315 

D  A  L  P  ( 1  1  ♦ 

61  = 

■  0456 

RA  (11, 

6)  = 

.1210 

PALP  M2, 

61  = 

.0336 

RA ( 12, 

6)  = 

.1153 

PALP ( 13* 

61  = 

.0170 

RA ( 13, 

6»  = 

.  1060 

0 AL  P  (  1  4  • 

61  = 

-.0051 

RA  (14, 

61  = 

.0965 

OALP  nS, 

61  = 

-.0343 

RA (15, 

6)  = 

.0*75 

PALP (16, 

61  = 

-.0690 

RA (16, 

61  = 

.0799 

PALP ( 17, 

61  = 

-.1036 

R  A  {  1  7  , 

61  = 

.0  /4M 

PALP ( 1*, 

61  = 

-.1266 

RA(1P, 

61  = 

.0729 

f»At  P  (  1R, 

6)  = 

-.1266 

RA ( 19, 

6)  = 

.074* 

PALP (20* 

6>  = 

In  36 

RA  (20* 

61  = 

.0799 

PAt  P (PI  , 

61  = 

-,0690 

RA (21 , 

6)  = 

.0*75 

PALP ( 22 , 

6)  = 

-.0343 

RA (22* 

61  = 

.0965 

OALP  <23, 

61  = 

-.005  3 

RA (23, 

6)  = 

.  1  160 

OA(  P(24, 

6)  = 

.O|70 

RA (24, 

61  = 

.1153 
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.3 

COMPUTATION 

OF  AFRO-ACOUSTIC 

pROPERT IES 

OE  SUPPRESSOR  NOZZLES 

CASE 

NO,  1 

cph 

7-IUHE  AR=? 

.3  NOZZLE 

-  VJ=2?00  EPS  - 

TT J= 1 600  OEG-R 

AXIAL 

LOCATION 

=  .0770?  ( X/OE  0  =  l. 00000) 

M 

P 

AN  OLE 

u 

rtf  mstty 

Tt  MP. 

ll/UREF 

THRH.iNT 

.  R/DEO 

i 

.oceoi 

o .  n  o 

2199.44 

•  0*004948 

1 ?  39.S7 

1 .00000 

.00037 

.00010 

i 

.00001 

1  r  .00 

?199.44 

.0009948 

1739.47 

1 . 0  0  v  0  0 

.00038 

.90010 

i 

.1.0001 

7  .  .  00 

2199.44 

. 0  C  0  9Q4  8 

1739.47 

1  .coooo 

.00040 

.00010 

i 

.COO  01 

30.  uC 

2199. r4 

.0009948 

1739.47 

1  .coooo 

.00043 

.00010 

? 

.  R07?9 

3. or 

1199.46 

.  OOOQ^ftR 

1739.47 

1  .  0  0  0  0  0 

. OC 1 46 

.  loooo 

? 

.40729 

1  3. 00 

7  1 1>9 .44 

.0009948 

1739.47 

1.00  0  0  0 

.00169 

. icooo 

? 

.00729 

20.0  0 

71 34.4F 

.(■009948 

1739. S7 

1 .coooo 

.00186 

.10000 

? 

.  00729 

33.00 

7 1 00,46 

.3009948 

1739.47 

1.00000 

.00180 

.10000 

3 

,61448 

? .  c  c 

2199.46 

.  1009948 

1739.47 

1.00000 

.00044 

.20000 

3 

.  .;]44M 

7.4  c 

> 1 Q^.ftft 

,0009948 

1739.47 

1.00900 

.00079 

.20000 

3 

.C  1448 

17.0.3 

71PP.r4 

•  0 0  09948 

1739.47 

1  .OO^OO 

.001?? 

.20000 

3 

.1 1448 

22. SO 

2199.46 

.0000948 

1739.47 

1  .C0300 

. 0PC9P 

.20000 

3 

.  1  4  5  R 

3  3.  30 

7  1  P  P  .  4  4 

.4009948 

1739.47 

1 .00900 

.00099 

.20000 

ft 

. J?l«8 

c.co 

71P9.7C 

..(OC9PSO 

1739.31 

.999(38 

. 00?69 

.30000 

u 

.  -2188 

7. S3 

7  l  P  P  .  7  ■' 

.0009940 

1739.31 

.99988 

.00303 

.30000 

u 

.  021R6 

1  S  .  o  0 

2  199.20 

.UG09Q40 

1739.31 

, 99988 

. 00?8p 

.30000 

ft 

.  '?1PR 

77. SC 

7  1PP.  7f- 

•  0;  00 99 SO 

1739.31 

.99988 

• 00?98 

.30000 

ft 

.  12188 

3  3  .  C  •) 

2199.19 

.0009940 

1739.31 

.99988 

. 00?93 

.30000 

c. 

.02917 

O.OO 

7  1  4S  «  PR 

.■'01  (3  IPS 

1709.44 

.9746? 

. 05?47 

.40000 

s 

.  '.?P17 

s.0  0 

714S.76 

. 0  C 1 0 1 94 

17C9.S? 

.97449 

. 0S?6 1 

.40000 

ft 

.  1?017 

12.0  ) 

714S.RC 

.1.1010144 

1709,48 

.97461 

. OS?S3 

.40000 

S 

.0291 7 

18.00 

714S.RC 

.0010194 

1709.48 

.97460 

•  05?4  3 

.40000 

s 

.'?917 

7r  .  yt; 

714S.  7S 

.'301C19S 

1  709. S? 

.97446 

. 0S?6 1 

.40000 

s 

. ^?P17 

30.00 

2145.83 

.3017144 

1709.46 

.9746? 

. 04?47 

.40000 

6 

. ‘  3646 

■’.00 

1  .378.4  1 

.(011447 

1 06  3.74 

.  6?68  0 

. 1 3781 

.40000 

ft 

.  *>  Iftftft 

6.0  0 

1.371 ,48 

.’01  1494 

1  ('6  3.44 

,62355 

.13764 

.40000 

6 

•  ■*  Iftftft 

i7.cu» 

1  37.3.  R7 

.(01  164? 

144?. 79 

,  6  ?  R  ?  8 

.13784 

.40000 

ft 

.  1 3646 

ip.  CO 

1 373. 47 

.001160? 

1 ; 4?. 79 

,b?428 

.13785 

.S00O0 

^  * 

ft 

.13646 

74. OC 

1371.48 

.8011 S9S 

1063.44 

.67345 

.13784 

.40000 

I  < 

6 

.  13646 

30.03 

1378.41 

.0011497 

1  (  6  3 . ?4 

.  6  ?  6  8  Q 

.13781 

.40000 

7 

.r-4374 

0.07 

?  0  3  .  S  4 

.0014983 

776.07 

.C9744 

•04??S 

.60000 

7 

.04374 

4.00 

707.17 

.  3C 14977 

7?6.S4 

.09189 

•  04?0? 

.60000 

> 

7 

.  'A  374 

10.03 

7(37.  14 

.001497? 

7  ?4 . 44 

.09190 

.  04?U? 

.60000 

7 

.-'4  374 

1S.0C 

703.49 

.U0149H3 

7  ?4 . 0  4 

.  09?57 

•  04??4 

.60000 

1  1 

7 

,^417£ 

70.  CO 

?07. 1 7 

,c  (314977 

7?4.S3 

.09)9? 

.04?00 

.60000 

I  ^ 

7 

.  R  4  3  7  c 

7S.3C 

737.  1  7 

7?6.S3 

.0919? 

•  Q4?00 

.6O000 

* 

7 

.  04  37c 

3  •'■.OC 

7*3. 6 1 

.  r.  0  1  *9ft  1 

7  Rft  #0^ 

.  l)9?4  7 

.  04??3 

.60000 

ft 

.  OS  1  04 

■>.GO 

2.42 

.0 071714 

481 ,?6 

.00114 

. 00?S3 

.70000 

) 

R 

..’4104 

S  .  i.i  (. 

7.43 

.7  071741 

48.  .41 

.60  1  IS 

. 00?9? 

. 70000 

8 

.  0  S  1  f 1 4 

1  0 . 0  c 

7.P4 

.(,071441 

474.83 

.00134 

,00?°4 

. 7P0O0 

ft 

.  -SI  04 

is.  33 

3.47 

•  >07 ) 4?7 

47..  1  7 

.00)48 

•  0  0  ?  8  4 

.70000 

ft 

.  (  S  1 

7  3.70 

3.44 

.'07149? 

Sft  H  #  ft  ft 

.00166 

.  00 ?H  7 

.70000 

£ 

.  SIC  4 

7  S  .3 

3.86 

.  0  ?  1  7  R  8 

444.98 

.00174 

. 00?67 

. 70000 

ft 

.  'SI  04 

3'  •  C  -7 

4.17 

.071 794 

444, 79 

.00)87 

. 00?78 

. 70000 

Q 

.  -SH37 

.03 

3.74 

.  071444 

47?. 04 

.00147 

.00144 

.80000 

.  48  3  3 

4 .70 

1  .  IS 

. .07074? 

494.?! 

.0006? 

.00191 

.60000 

Q 

.  '  S«  33 

M  .  S  7 

C.,,0 

.  . ■  0 ?? 8  34 

{3ft  0^ 

8.00000 

0.00000 

.60000 

Q 

.  S  ft  ft  1 

1  7  .  H6 

o .  f  n 

.  0  0?  ?R  14 

^>4 t'  .;■>!» 

0 .  r  n  n  1)  0 

ii  .00(00 

.ho  non 

144 

_ 

I  >  •  M  fi  H  *  *  * 
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COMPUTATION  or  Af.PO-ACOtlSTIC  PROPFPTlfS  or  SUPPRFSSOP  N077LtS 

f  N0<  ,  CRD  7-TIJHe  AR=?.3  N077LF  -  VJ=??G0  FPS  -  TTJ=l600  PFC, 

AXIAL  LOCATION  =  .0779?  (X/OFO  =  1 .0000  =  > 


M 

p 

ANGLE 

II 

DFNSITt 

g 

. -SHIT 

17.18 

0.(.,0 

. 0  rt??8  3S 

g 

.  rsfl37 

71.83 

O.oy 

. 0 0??8  3S 

g 

. -SATA 

7S.71 

.007783S 

9 

.  1S8  13 

3  .■  .  uv 

O.c  « 

.9C7783S 

1  C 

. , 8  =  61 

3 . 0  v 

793.7? 

.00 1A98A 

1  0 

.  :asa? 

8.39 

1  A  8  .  **  C 

.1017879 

]  0 

."hShl 

8  .  S  7 

OS.  (8 

.  9  0  I  ft  A  7  A 

1  c 

.  .iFSf'l 

17. 8A 

77.89 

.0070377 

1  0 

.  ASAl 

17.1« 

S.80 

.1.071  A  1A 

\  f) 

.  ASA  3 

71.87 

0  .  ■  • " 

. ' 0  ??8  3S 

\  0 

.’ASaT 

?S.  71 

o. 

.,0778  78 

1  0 

.  'ASA  3 

7  3.00 

0.  '■ 

.  10778  IS 

1 1 

.1  7?43 

i.cr 

1  378. 8S 

.  901 1S98 

1 1 

.:>7?9? 

3.7S 

178?. n 

.901 1 798 

\  1 

,C7?9? 

7.  SO 

1  0 1 S . 37 

.0013  19a 

1 1 

.47793 

11. 7‘- 

(.08.89 

.00!  388 a 

1 1 

.0739? 

is.ro 

387.7A 

.3  01  AS. 18 

1 1 

. ^7 PR? 

18.7  = 

S  7  .  ?c 

..'019S8S 

1  \ 

.4739? 

73. sc 

A  .88 

..'031  A3  7 

V  \ 

.,>779? 

38. 3S 

0 . 1  0 

.  .,0??83S 

1 1 

.A 770? 

78.  (,0 

ft  .  '  3 

.1>0??«3S 

1? 

.48031 

0.0' 

318S.78 

.'■0101RA 

1? 

. )H0?1 

7.7S 

3138. S8 

. 0  0 1 0?MA 

1  ? 

.08071 

7  .  c  0 

707S.18 

. J018S7S 

)  ? 

.  484?] 

1  1  .  3S 

1717.89 

.">0  11  9S7 

1? 

.4ftp?( 

l  S  .  C  0 

1 8S9.87 

.0013773 

\? 

.08071 

18. 7S 

3=3.98 

. 001S887 

\  ? 

. 78071 

3  3  .  S  0 

83.70 

.0019AAS 

\  ? 

.  ’•  A  0  7  1 

3A.3S 

?.?1 

.0 071771 

1  ? 

.  8  0  ?  l 

30.OC 

o.-,n 

.0077839 

\  1 

. 087  so 

0.0  4 

3199.19 

.  C-  0  0  9  9  S  0 

\  3 

.  0  8  7  S  0 

1.  7  ? 

?I 98.98 

.  v  C  0  9 9  S  ? 

1  3 

.  8  7  S  0 

A.  A  7 

3  1  9a  .  At* 

.  0  0O99A7 

1  3 

. o  8  7=  0 

1  O.0C 

7178.80 

.i’OIOO^a 

\  3 

.08790 

17.37 

3077.8a 

. 00 1 OS08 

1  3 

.  0  8  7  S  0 

1  A.A7 

1S00. 1A 

.(•0  11779 

\  3 

.  8  7  S  0 

7C.00 

A 1 A  .  3S 

.90 1 7H?ft 

l  3 

.087=0 

37.77 

1  1  C  .  ?  1 

•  0  Q187A7 

1  3 

. A87S0 

7  A  .  A  7 

7.8S 

.' 071SS? 

n 

.-J87S0 

7J.V0 

0  . 0 

.  077838 

14 

.'•<78  70 

8 . 8  c 

.•>199.8  = 

.,1009988 

14 

.(9870 

7.37 

31 99.8S 

.  ,009988 

1*4 

. AQ879 

A  .  A7 

3 l °9 . 88 

.(  909988 

\  4 

. 098 79 

13. GC 

3198. 7S 

.90099S3 

14 

.  i*  Qft  79 

17.  77 

317S.18 

.9010078 

14 

.  .198  79 

1  A  .  A7 

1911.83 

.  '01078S 

14 

.  (98  79 

70.00 

10??.  8 

..01  7777 

14 

. 19879 

73.77 

70  8.7't 

.  00 1 A9?9 

TI  MP  . 

U/IJPEF 

Tl'RH.lM. 

P/Df  0 

S8  .1.0  9 

O.OOoOO 

o .  o  o  o  o  o 

.AOOOo 

bftO  ,  0>' 

O.OCftOO 

O.OuOuO 

.80040 

S8  *,  .  u  0 

9.00030 

0 . 0  0  0  0  0 

.80000 

SftO.O. 

'1.00000 

O.OoOOO 

.80000 

77S.98 

. 0938 A 

.0831  A 

.90000 

707.88 

.07A38 

. 0  37 1 C 

.90000 

AA.;.?S 

.  08  3  7S 

.03887 

.90000 

A(,  A  ,  7  7 

.  0 1 S?  3 

.01 0A7 

.90000 

984.01 

• 90?AA 

.00  33S 

.90000 

Sft  .  .  C.  4 

ft.00030 

0 . 0  0  0  0  0 

.90000 

S8-  . ,  ■; 

?  .  0  u  0  0  0 

0. 00(100 

.49000 

S8  ..(9 

O.OuvOO 

0.000O0 

.90300 

1  1A  1.  19 

.  A?A7? 

.13779 

1.0(1900 

1 389.07 

.  S  8  3  o  7 

.13990 

1  .1,0000 

498. 7S 

,8A 1  AS 

. 1 3AA8 

1 .90000 

HQ-  .  S  8 

.37978 

.10987 

1 .toooo 

787.17 

.110  3  7 

.0S481 

i  .coooo 

A?9.6l 

.0,?A00 

.01888 

1  .00000 

S  A  9 , 9 

.00793 

.00379 

1 .coooo 

98.  .00 

0.0 0000 

O.OOOCO 

1.00900 

S8-'.9f. 

o.coooo 

0 . 0  0  0  0  0 

i  .oftooo 

1709.8? 

. 97SS8 

.  (,S?88 

I .10000 

1791.18 

.9AS98 

. 0 A88 l 

1.10000 

1 171. AO 

,9?r  79 

. 1 38 AO 

1.10000 

1 1  IS. AH 

.78087 

.1SHS? 

i  .  10000 

K.08. 07 

.88188 

.  1  A  38  0 

1.10000 

798.77 

. 1AC78 

. 08S33 

1  .  1  0000 

A77.C3 

.07881 

.01880 

1.10000 

Sa  7 . A9 

.00100 

.00338 

i .  ioooo 

S8  4.0'. 

'J  .  C  0  ‘J  0  o 

3.00000 

1.10000 

1739.7'. 

.99988 

.00013 

1  .70000 

1779...7 

.99877 

. 0  0  3S0 

1 .70000 

1777.1? 

.9987? 

.Cl  18a 

1  .70000 

17?  1.8. 

. 4897S 

.08313 

1.70030 

1173.87 

.97808 

.11890 

1.70090 

1087, A! 

. A871A 

.18610 

1  .70000 

ft°l .78 

,  38  v,  3  3 

. 1 3  o?9 

1 .70000 

A71 .78 

.OS-Ill 

. 0  3SA A 

1.70000 

S77.18 

.  00  3S  7 

.00863 

1.70000 

S8I,  .9(1 

0.00030 

o.ooooc 

1 .70000 

1739. S7 

1 . 00O00 

.00140 

) .30000 

1774. S7 

1.00000 

.00)61 

1 .30000 

1779. SS 

.94499 

.00087 

1 .30000 

1778,9= 

. 994A8 

.  0  0  S  8  3 

1 .30000 

1 773.88 

.98897 

.(,8688 

1.30000 

1 187. A? 

.  8  A  4  4  S 

. ISAHft 

1.30000 

4  9  A  ,  S  7 

,  ft  Aft  A  9 

.  1  7  A?  1 

1 .30000 

738.79 

. 04ftH9 

.  969  79 

1 . 39000 

o  «  *  M  Cl  H  o«o 
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COMPUT  All  ON  OF  AFRO-ACOUSTIC  PROPERTIES  OK  SUPPRESSOR  N07ZLE  S 

CASE  NO,  1  CRD  7-TIIRE  AR=7.)  N07ZLF  -  VJ=??00  FPS  -  TTJ=|600  OEG-R 


AXIAL  LOCATION  =  .  779?  ( X/OEO  =  l.OOOUO) 


M 

R 

ANGLE 

IJ 

I3FNSI  TV 

TFMP. 

U/UREF 

TURH. INT 

.  P/DEG 

1  A 

.,-.0470 

76.67 

14. 7S 

.  07116S 

687.69 

. 00671 

.00707 

1.30000 

18 

. >0470 

1 3 . 0  3 

0.09 

•  0  0778  JS 

640 .00 

0.00030 

C  .00000 

1.30000 

IS 

.  1  0708 

3  .  C  0 

7199. 4S 

.  100.9948 

1719.67 

1.00000 

.00031 

1 .40000 

IS 

.  1  0708 

1.  C  V 

7199. 4S 

.'009948 

1719.67 

1.03000 

.00186 

1.40000 

IS 

.  10?C« 

6.00 

?199.4S 

019948 

1719.67 

l .OOOPO 

.00140 

1 .40000 

IS 

.10708 

9.00 

7199.45 

.0009948 

1719.66 

1 .00000 

.001S9 

1 .40000 

IS 

. 1 0708 

17.03 

7198. SS 

.C0099S4 

1718.7? 

.99969 

.00773 

1 .40000 

IS 

.  1  0708 

IS.  00 

7161.61 

•  -j  o  i  o  l  is 

1716.6? 

.96780 

.06164 

1 .40000 

IS 

.  1  C?0« 

18.00 

1801 .S| 

.101 09? 1 

1179.1? 

.81907 

. 1 7666 

1 .40000 

IS 

.  1  0708 

71 .00 

894.46 

.6017917 

964.61 

.17919 

.16066 

1 .40090 

IS 

.  1  O?08 

74.00 

146.71 

.3017718 

696.16 

.06670 

.04583 

1 .40900 

IS 

,10?0« 

77.00 

9.4S 

.0071459 

674.61 

.00430 

.00616 

1 .40000 

IS 

.  1  0708 

11.0  c* 

O.C  9 

.  '"3??8  IS 

54/ 

o.oocoo 

0.00000 

1 .40000 

18 

.  i  o<n« 

0.00 

?  1  99 . 4  S 

.0009948 

1779.67 

l .oocon 

.00104 

1  .50000 

18 

. 1  OPTS 

l.cu 

7199. 4S 

. Ifl 1Q948 

1719.67 

1 .00800 

.0077? 

1.60000 

18 

,10938 

6. 00 

7199. 4S 

.  1009948 

1719.67 

1.00000 

.0077? 

1 .50000 

18 

.10918 

9.0' 

7199. 4S 

.  .(00994  8 

1779.66 

1 .00000 

•OC 1 8? 

1 ,6(,000 

18 

.10918 

17.00 

7198.81 

. '009961 

1 778.96 

.99971 

.00679 

1.60080 

18 

.10918 

1 S .  0  3 

7161.18 

.  4  0  1  0  1  J  7 

1716.41 

.98760 

.06149 

1 .60000 

18 

.10918 

18.00 

17S7.9M 

.0010996 

1171  .4.1 

. 79660 

.1815? 

l .60000 

18 

. 1 0918 

71 .03 

710. 4S 

.  101  7.194 

9?l .61 

.17331 

.14631 

1 .60000 

18 

.13918 

74.0' 

98.76 

.001861S 

667.41 

.04468 

.03769 

1  .50000 

18 

.10918 

77.  CO 

4.71 

. '071816 

S64 ,70 

.0019? 

.00363 

1 .50000 

16 

,10918 

13.00 

O.C  9 

.(0778  IS 

S  4  f .00 

O.OCOOO 

o. 00000 

i  .soooo 

17 

. 1 1667 

3.00 

7199. 4S 

.0009948 

1719.67 

1 .00090 

.00706 

1  .60000 

17 

.1 1667 

7.71 

?199. 4S 

.0009948 

1739.67 

1.00900 

.00739 

1 .60000 

1  7 

.11667 

S.4S 

7199. 4S 

,C 009948 

1779,57 

1 .00090 

.00178 

1.60090 

17 

.11667 

8. 18 

7199.45 

.,(009948 

1719,66 

l.OOPOC 

.00167 

1  .6000  8 

17 

.  1  1667 

19.81 

7198.91 

.  } 0  0  99S7 

1/19,14 

.999  75 

.00557 

1  .60990 

17 

.1 1667 

13.64 

7174.93 

.0010080 

1771.34 

.96666 

.04751 

1.60000 

17 

.  1  1667 

16,16 

1 916. IS 

.'010751 

1146.71 

•  866  74 

.15670 

1  .60000 

17 

.  1  1667 

19.39 

1  389. r-1 

.001731? 

999.94 

.47741 

. 1 7474 

1 .60000 

17 

.11667 

71 .87 

73S , 1 7 

•  0 1)  1  6  6  y  8 

747,44 

. 1 0 687 

.06631 

1.60000 

17 

. 1 1 667 

74.55 

7'  •  7  1 

. 1 0  ?0  8  S6 

691 .7? 

.9094? 

,0087s 

1  .60000 

17 

. 1 1 667 

77.77 

O.'O 

.0077816 

64' .00 

C.PoOOO 

C .00000 

1  .  6000  C 

1  7 

.  1  1667 

1C  .09 

0 . •) 

•0C??81S 

64  r ,00 

0.00000 

C. ooooo 

1 .60000 

1  H 

.17196 

9. 1C 

7199. 4S 

.3009948 

1719.67 

1 . 30000 

.00164 

1.70000 

18 

. 1 7196 

7.71 

8104. 4S 

. 0009948 

1719.67 

1 .00000 

.00141 

1.70000 

18 

.17196 

S.4S 

7199.48 

. 3009948 

1779.66 

1 .90000 

.0016? 

1 . 70000 

18 

.17196 

8.18 

7199.13 

.1009960 

1779.76 

.94986 

.00436 

1.70000 

18 

.17196 

10.91 

7188.71 

.301031 1 

1 7  7 1 .46 

.9961? 

.07817 

1 .70000 

1« 

. 17196 

1  1.64 

70S  1 .70 

•  vlO  1  0466 

1 1 78.11 

.91769 

.11403 

1.70000 

18 

.17196 

16.16 

1471.87 

.301 1478 

1 174.76 

.6519? 

. 1 8468 

1  .  70000 

1  8 

. 17196 

19. C9 

478. S' 

.00146/4 

646. C  9 

.71766 

.10766 

1.70000 

18 

.17196 

71.«7 

67.84 

.0019401 

616.66 

. 37667 

.07146 

1 , 70000 

18 

. 17196 

74. 5S 

7.8-3 

.0071966 

66) .17 

.00177 

.0034? 

1 . 70000 

18 

.17196 

77.77 

n  .  »>  - 1 

.  '.0776  )6 

54  (-.00 

O.OOOoO 

0.00000 

1 .70000 

1  8 

. 17196 

1  9  .  P  ' 

r 

.  ,0778.16 

64  '  .  f  ,1 

0  .  f,  9  f.  0  0 

r  .  n  o  o  o  r. 

1  .70(100 
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COMPUTATION  OF  AERO-ACOUST I C  PROPERTIES  OF  SUPPRESSOR  N077LES 
CASE  MO.  1  CRO  7-TURE  AR=2.3  N077LE  -  VJ=2?00  FPS  -  TTJ=l6O0  DEG-P 


AXIAL  LOCATION  =  .3779?  (X/OEO  =  l.OOPOOT 


M 

P 

angle 

II 

OF  MSI TY 

TEMP. 

ll/UREF 

TURH. INT 

.  R/DEu 

IS 

.  m?s 

0 . 0  c- 

7  I  OR .  18 

.008090 

1239.3; 

.99948 

.00379 

1  .80000 

IS 

. 1 317S 

7.88 

7198. QP 

. TOl/9991 

1739.17 

.99979 

.0081? 

1 .60000 

19 

.13179 

9.00 

7197. A? 

. 300996? 

1737.79 

.99908 

.00976 

1  .80000 

19 

.13179 

7.83 

7189. AL 

.00  IGF  30 

1279.37 

.99361 

.07998 

1.80090 

19 

.13179 

l  o .  c  r 

7  1  ('  6  .  A  8 

..010371 

1 1 9A .73 

.9977? 

.08886 

1  .80000 

19 

.13179 

17.90 

1  776.93 

.UG 10963 

1  1 7A . 8  ^ 

.60771 

.15889 

1 .80000 

19 

. 1 3179 

19. CC 

1  839.70 

.'017377 

1  wO'j  .33 

.87703 

.1618? 

1 .60000 

19 

.13179 

l  7.93 

318.73 

.O0197A3 

783.79 

.  1  8  A  9  1 

.07710 

1.8000  0 

19 

.13179 

7  li .  0  C 

A  7 . 6  A 

. 0019899 

6?  J . 9 1 

.07166 

.01638 

1 .80000 

19 

.13179 

7?.9' 

7.89 

. '071988 

96... 79 

.CO  131 

.00370 

1.80000 

19 

.13179 

79.30 

o.  ;o 

..  0778  39 

98, -,.0  3 

0. OGCOO 

0. OOOOO 

1 . 8  o  0  0  0 

19 

.13179 

^7. 9  C 

9.-0 

.:'0??8  39 

98<;  ,  f;  •) 

O.OOi  ,>0 

0.00000 

1 .80010 

19 

.13179 

3  0.0  0 

0.  -  7 

.0 0778 39 

98  1  ,  <) 

1.001 00 

u. ooooo 

1.600  >0 

20 

. 1 3894 

o  .or 

71A9.88 

.'010196 

1709. 39 

.97999 

.05799 

1.90000 

70 

. 1 3894 

7.90 

717R.49 

. '01074P 

1703.69 

.96839 

.061)6 

1.90000 

70 

. 1 3«94 

9.00 

7081.96 

. >010838 

1181.3- 

.  9  3  7  A  9 

.08886 

1  .90000 

?0 

. 1 389a 

7.93 

1887.77 

.30106?? 

1  1  39.8r- 

.  8  A 9? 0 

.13739 

1.90000 

70 

.1389a 

1  0.00 

1836.17 

.•.011863 

1  0  7  3 . 86 

.69299 

.16179 

1 .90000 

?'J 

.13894 

17.9' 

786.91 

.0,01  371  ? 

933.33 

. 3A  399 

.17913 

1  .90000 

70 

. 1 399A 

19.00 

731.78 

.-•,010,676 

781 .66 

.  1  f.537 

.05730 

1  .90000 

70 

. 1 389 A 

17.93 

37.77 

.’■070169 

8)  ,  ,9V 

.11717 

.01799 

1  .90000 

7  0 

. 1 389 A 

70. C( 

7.99 

.107701  7 

88 -.3? 

.00118 

.00319 

1  .90000 

70 

. 1 389a 

77.99 

0.  •  8 

.  .  0776  39 

98.,,  r.  ) 

3.00100 

6.00000 

1 .90000 

70 

. 1 389a 

79.  18 

0.  o 

.00776  39 

98  >,  )( 

C.OOUOO 

0. 00000 

1 .90000 

70 

. 1 389 A 

77.9." 

0.  -4 

. 80??839 

98  '!  .  A 

1.00030 

0.00000 

1  .90000 

70 

. 1 389 A 

37 .50 

0.-4 

. u0??«39 

98' .10 

0  .  G  0  •)  0  9 

0. ooooo 

1  .90000 

71 

.19983 

'.I" 

1  378.1  1 

.'01 19PP 

146 3. lo 

.6769  7 

.13781 

7.00000 

71 

.  1  9983 

7.31 

1 303.87 

..'011719 

1  v‘8?  .9? 

.99781 

. 1 3796 

?. ooooo 

71 

.19983 

8.8? 

l  1 00.88 

.'>0  1  71  70 

1613.17 

,9Gu8 3 

. 1 3898 

7.0(1000 

71 

.18981 

6.9? 

773. <37 

.'•01  3  1  A  9 

9 36.09 

.39149 

.11871 

7. OOOOO 

71 

.18983 

9.73 

8  1 8  .  ft  0 

.v,0)8:k? 

871.91 

.18873 

.08099 

7. OOOOO 

71 

.18983 

11.84 

19 A.  1  8 

.'017639 

8  9  9  .  *0  8 

. C  70  09 

.03861 

7. OOOOO 

71 

. 1 a9R3 

1  3.88 

37.1? 

.1  070704 

81  >.  .  3? 

.01688 

.01701 

7.CO0O0 

71 

. 18983 

18.19 

A.  07 

.007)878 

98A.99 

. 0  0??8 

.00378 

7. OOOOO 

71 

. 18983 

18.  Ah 

0. 

.  07? 8 39 

68, ,0j 

C. 00000 

0. 0O000 

7.0  0000 

71 

. 1 AS8  3 

79.77 

0 . 8  9 

. ' 077839 

98  -  .in 

0.90000 

y.OoGOO 

7.0C0G0 

71 

.18983 

7  3  .  8 

I).'-' 

. ' 0778 39 

9  A -.I,,, 

1.090G0 

0.00000 

7. OOOOO 

71 

. 18983 

79. 38 

0  .  ••  4 

.  .  0776  19 

98  ,.  jfi 

o. ooooo 

•'.00000 

2. OO000 

71 

.18983 

7  7.1-p 

0.  ]<) 

. 1077839 

-  8-  .00 

2.0000  0 

0.00000 

7. -0000 

71 

.18983 

li'.or 

0.38 

.'077839 

980. 0 J 

0 . 0  0  u  •)  0 

0.00000 

7. OOOOO 

7? 

.19313 

5.  IT  A 

703.77 

.•>01  6986 

779.89 

.69787 

. 08770 

7.  1  0000 

7? 

.19313 

7.  1  A 

181.19 

. -017738 

719. AO 

.08736 

.03879 

7.1  0000 

7? 

.19313 

4.79 

111 .48 

.'018011 

8«A ,89 

.06001 

.03089 

7. 1 0000 

7? 

.19313 

8.  A  3 

7  3. 7Q 

.  ->0191  A3 

6 8 A. ?9 

.C  331? 

.01917 

7. 1  00O0 

7? 

.19311 

9.97 

11.17 

.  > 070 A? 7 

603,66 

.( Ip  1  7 

.00999 

7  .  1  0  0  o  0 

77 

.19313 

18.71 

<».  '3 

.  ,0?I9?a 

977.60 

.004  1  ? 

.00  396 

7.  1  0090 

7? 

.19313 

1  7  •  9  6 

1  .  1  T 

.071983 

961 ,AA 

.00480 

.00774 

7.10000 

7? 

•  lr-313 

19.3-, 

n .  -  3 

.■•07,>8  19 

9  A  :  .  1  1 

fi .  r  o  o  o  o 

o.o  in  or. 

7.16000 
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COMPUTATION  OF  AFPO-ACOUST  IC  PPOPFPT  IPS  Of  SUPPPESSOP  N077LES 


CAST 

NO.  ) 

CPO  7-TI|Rf  AR=7.3 

N077LF 

-  VJ=??00  FPS  - 

TT  J= l E00 

axial 

LOCATION 

=  ,.7?9?  (X/OEO  =  l.OAQJ' 

■  1 

M 

P 

ANC-l  K  tl  IiPNSITY 

TEMP. 

IJ/URFF  TURR.INT 

.  P/OFu 

7? 

.  1  S  3 1  7 

1  7  .  1  A 

1  .  '•  '• 

. 7077830 

04  -,.  Cu 

C. 000 00 

D.OOuGt 

7.1 0000 

?? 

.  1S31  7 

H.7Y 

0  . 

. 3077830 

040. 00 

0  »  0  C  0  0  0 

O.GOCOC 

7.10000 

7? 

. IS317 

7  1  .  A  3 

1.” 

.' 077830 

04  .  ,r 

n.ooooo 

0.00000 

7.10000 

?? 

. 1S71  7 

7  3.07 

o.\  ' 

.,:C7?‘330 

S4. .38 

o  *  o  o  o  o  o 

0.00000 

7.10000 

7? 

.10717 

70.71 

0.  '  ' 

•  ■  3??h30 

04. 

0 . 0  0  0 1)  0 

(1,00000 

7.10000 

7? 

.10717 

77.  ME 

0.  ■' 

. ; 077830 

04  .  0  0 

E.COOOO 

0.00000 

7.10000 

7? 

.10717 

1 J . ')  0 

0.'r 

.  8>??h  m 

04  y  .  C  0 

3.  OCC'OO 

0 .00000 

7. 10000 

73 

.  IMii.? 

3.00 

4  .  E  1 

.  "'071  800 

OE  3.70 

. 0C2  l  0 

.  0  0?f*o 

7.70000 

77 

.10007 

7.10 

3.77 

.  '071  <344 

0E1  .<37 

.0017? 

.007PS 

7.70000 

73 

. IOuO? 

o  ,?P 

7.10 

.  o  0  7  7  0  7  E 

008. OE 

.  COC'RP 

.00776 

7.70000 

73 

.  1  eoa? 

E  .  4  7 

0.  .E 

. rc7?7jo 

04;-.  C>) 

).CC  I'OO 

0.00000 

7.70000 

7  3 

. 1 eja? 

i.^7 

0. 

.  i' 0  778  70 

04. .0' 

« .  0  0  v  0  0 

0.00000 

7.70000 

77 

.  10007 

1  ..71 

r  *•  t 

•J  • 

,70??h3O 

c-4< 

O.UCGOO 

o.ooooo 

7.70000 

73 

. IE  00 7 

17. EE 

0  .  .  f 

.  •.'C  778  30 

04' .or 

■j.00000 

0.00000 

7.70000 

77 

. 1E0O7 

lO./'^ 

0.  ' 

. :c?783o 

04'. 00 

0 . 0  0  0  0  0 

0.00000 

7.70000 

77 

.  1  Mo? 

17.1.. 

o.. 

.  .MV7H30 

04  r, 

o.ooooo 

o.ooooo 

7.70000 

73 

.  1604? 

1  r.?8 

.  0  7  7  8  3  o 

04  .  0  3 

0.00000 

0.00000 

7.70000 

73 

.  lEOo? 

71.43 

0  • 

.  <077830 

04 7. CO 

0 . 0  0  0  J  0 

0.00000 

?  .  7  O  0  0  0 

73 

.  1  E()4? 

7  3.07 

n.  ■ 

. ‘077830 

04  ■  .  r,  J 

O.OoyOO 

o.ooooo 

7. 71.000 

73 

.  1  E^a? 

?r  .  7  1 

. e  077m  30 

.  00 

0. 00300 

O.OCGCO 

7.70000 

77 

. 1604? 

77.  on 

0. 

. < C??6 30 

047.07 

0.00000 

0.00000 

7.70000 

7  3 

. INTO? 

1  '  .  0  ^ 

0.  .1- 

.f  •)??«  30 

O4...00 

O.OOoOO 

C. 00000 

7.7000U 
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CIPCUMFEWENTI  ALLY-AVEPAGFl)  PARAMETERS 


NP 

RAOIIIS 

MACH  NO. 

TEMP. 

INTENSITY 

frequency 

1 

.000  1 

1 .96*? 

7.7955 

.  66  ')C  6F- 1 7 

0. 

7 

.  1000 

1.966? 

7.7955 

.74773F-06 

3. 

3 

.?ooo 

1 .966? 

7.7956 

.96416F-08 

c. 

4 

.  7000 

1  .9661* 

7.7961 

.71 3R4F-04 

7. 

5 

.4?00 

1 .91 a? 

7.7398 

.16103F.Q5 

7416. 

6 

•  SOO" 

1 .7776 

1  .9688 

.171 9HF ♦  Q8 

75958. 

7 

.  a  o  o  •> 

.1809 

1 .345? 

.51 J4  3E ♦ 04 

16507. 

a 

.7^00 

.  oo?a 

1 .0583 

.39HP5F-04 

4377. 

9 

.  a  o  o  o 

.0019 

1 .0777 

.41734F-06 

574. 

10 

.9000 

.1777 

1.7755 

. 1 0304F ♦ 04 

5139. 

11 

1 .OOCO 

.  0979 

1 .8757 

.  1  1  874F  *08 

17317. 

1? 

1.1000 

1  .SS73 

7.1133 

.?9016F*08 

5707. 

1.7 

I .7000 

1 . 773? 

7.7014 

.4771 3F*08 

7695. 

1  4 

1 . 70CO 

1 .7745 

7.7747 

.4O008E*08 

1676. 

15 

1.4000 

> . 79h4 

7.7343 

.•ospre^or 

1376. 

16 

1  .  S  O'  r,  0 

1 .6074 

7.7367 

.479 1 8F  *08 

1313. 

17 

1  .  a  o  o  r 

1.7967 

7.7339 

#  3  8  6  ?  6  F *08 

1513. 

ia 

1 . 7  0  CO 

1.7761 

?. 7766 

.43643F  *08 

18P8. 

IP 

1  .  Hfi  C  0 

1 .7716 

7.7014 

.300  37F  *08 

7593. 

?o 

1 .9000 

1 .6589 

7.1139 

. ?4 1 84F »0H 

4755. 

?l 

7.  >000 

.ag?Q 

1.8767 

.  14  34  3F *08 

8933. 

a? 

7. 1  00A 

.1776 

1 .7750 

•  1  6  .I46E  ♦  04 

34  36. 

73 

7 . 7  0  0  0 

.0079 

1.0397 

.19614F-04 

733. 

WARNING 

-  NO . 

OF 

TURNING 

P0INT6  Is 

gpeatf.r  than 

7  AT 

K  A=  1 

X  = 

.0779? 

ITH=  IT 

1HFTA= 

110.00 

NTP  = 

3 

WARNING 

-  NO. 

OF 

turning 

POINTS  IS 

grfatep  than 

7  AT 

XA  =  1 

X  = 

.0779? 

ITH=  11 

The  TA  = 

170.00 

NTP  = 

3 

WARNING 

-  NO. 

OF 

TURNING 

D0 1  NTS  IS 

GREATfR  T  HA  >* 

7  AT 

<  A=  1 

X  = 

.0779? 

1 T  H=  1? 

THfTA= 

130.00 

NTP  = 

3 

WARNING 

-  NO. 

OF 

TURNING 

POINTS  IS 

greater  than 

7  AT 

II 

<x 

* 

X  r 

.0779? 

F  T  u=  13 

THE  T  A  = 

140.00 

NTP  = 

3 

WARNING 

-  NO. 

OF 

TURNING 

POINTS  IS 

GREATER  THAN 

7  AT 

X  A=  1 

X  = 

.0779? 

ITH=  16 

THE  T  A  = 

150.00 

NTP  = 

3 

WARNING 

-  NO. 

OF 

TURNING 

POINTS  IS 

GREATER  THAN 

7  AT 

*  A  =  1 

x  = 

.0779? 

1 T  Hr  16 

THF  TAr 

160.00 

NTP  = 

3 
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«(  11=  .07P9  IIHM  11=  .?N899E*?I  FM(  Jl=  . ?8 1 9E • 0 1  UAVG<  11=  1887.59  UMAX  (  1)=  P199.86 

V  ARM I  NO  -  NO.  OF  TURNING  POINTS  IS  GPFATFP  THAN  ?  AT 

*A  =  ?  X  =  .09187  1  TH=  10  THFTA=  110.00  NTP*  3  .......  ... 


WARNING  -  NO.  OF  TURNING  POINTS  IS  GREATER  THAN  7  AT 

KA  =  ?  X*  .09187  !TH=  n  THE  T  A=  1?0.00  NTP*  3 


WARNING  -  NO.  OF  TURNING  POINTS  IS  GPFATFP  THAN  ?  AT 
K  A*  ?  X=  .09187  1TH=  1?  1  HIT A=  130.00  NTP* 


EARNING  NO. _0F  TUPNING .POINTS  IS  GREATER  THAN  ?  AT  _ 

x  A  *  7  *=  .09187  ITH=  13  IHETA=  180.00  NTP*  3 


EARNING  -  NO.  OF  TURNING  PO|’  TS  JS  GRFATFR  THAN  ?  AT 

x  a*  2. _  x=.  _.ayit<7_  ith=  i8.  iheta=  isu.oo  _niE£ _ i 


<<  ,>)=  .0919  nA  |  |  ,’l=  .r»3MSF«Pl  FI-  (  ?>  =  .PS?9f.M  UAVG<  7 » =■  1 7G7.S.I1  U*AX<  2)=-2199.A6— . 

WARNING..-  MO.  OF  TURNING  POINTS  IS  Gur  A I  E°  THAN  7  AT  _ 

XA  =  1  X=  . | 1 9  7S  ITH=  n  TH»1A=  11". ""i.  NTP»  3 


W&PMMG  -  NO.  OF  TURNING  POI'TS  IS  <.R  F  A  T r  9  THAI  7  AT 

xA  =  3.  ..  x  =  .  .  11575  .....  1  T h=  1L  .  .  I  ME  1 A  *  I.S.OO  __  _NI££ _ 1 


FARMING  -  MO.  OF  TURNING  POINTS  IS  GREAT! p  THAN  ?  AT 

►A=  3  X*  .  1  ISIS  ITH:  1?  TMfTA=  IU.  00  NTP=  3 


WARNING  -  NO.  OF  TURNING  POINTS  IS  GREATER  THAN  7  AT 
K  A  *  3  x=  .11S7S  ITH=  13  T  HI T  A*  180.00  NTP*  3 


WAPNtNG  -  NO.  OF  TUPNING  POINTS  IS  GREATER  THAN  p  AT 
i  x  A  =  3  X=  .11S7S  ITH*  1 8  THETA=  ISf.OC  NTP”  3 

> 


»(  31=  . 1  I S7  U8l(  3 1 =  ,?NT9?r*?l  FP<  3*=  ,?889f.*Ol  UAVGl  3>=  1878.56  UMAX (  31=  ?)99.85 


WARNING  -  NO.  OF  TURNING  POINTS  IS  GREATER  THAN  >  AT 
x  A=  8  X=  * 18SH  3  ITH=  13  THETA*  1  1C . 00  NTP*  .3 

J 


WARNING  -  NO.  OF  TURNING  POINTS  IS  GREATER  THAN  ?  AT 

XA*  8  X*  . 1 8S8  3  ITH=  11  THtTA*  IPO. 00  NTP*  1 
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. 1 1854 

/.k: 

1 APO.P7 

.  101  IN  8 

1  A  A  .  ?  1 

.  ft  A  S  A  3 

.  15099 

1  .98010 

PI 

.  1  IPSA 

l'.OO 

1  1-'  1  .PS 

.  7  0  1  ?  1,  S  P 

1  >?1. 1 1 

.ssi so 

.  ISftP? 

1  .  911000 

PO 

.  1  IPSA 

1  P .  S  C 

11 1ft. 77 

.50  1  PAS  7 

SQ1.AA 

.  S  )  ft  8  A 

. 1SS77 

1 .9001,0 

p  ■ 

.  1  IPS A 

1  s  .  0  3 

1r-i,,7A 

.  01PHOG 

SSA. 1 7 

. A  3  3SS 

. 1 4844 

1 .9f-0u0 

Pc 

.  1  IPSA 

l  7.c  ■ 

/AG.  IS 

.  '0  1  is  10 

si  1  ,7G 

.  .3  A  8  .14 

. 1 3455 

1 . 9000(1 

PO 

.  1  IPSA 

Ps. O' 

SwS.i  p 

.  I01A1PC 

u  G 1  ,  r,  7 

.Pftftf 1 

.  1  ISAS 

1 .9(.000 

PC 

.  1  IPSA 

PP.s- 

p  ps , A  ? 

.  OlSPio 

•:  S.AS 

.  1 s3a? 

.  1 9  3  5  P 

1 .90000 

PO 

.  1  IPSA 

PS.i  5 

POG.PS 

.'•01  Cl  3P 

7  ft  A  ,  IS 

.  1 14Sft 

.  0  7  0  7  S 

1  .90000 

Pc 

.  1  IPSA 

P7.s- 

POP.'  1 

.00’ A  7  —  H 

7  7A  ,  AH 

.  7-9457 

.  0 APPft 

1 . 9  0  o  (  0 

?o 

.  1  IPSA 

I-.  i: 

1 7A.PG 

.'•01  asm  7 

7  PS. A  7 

.C7SP3 

.0.10  1  9 

1.90000 

PI 

.  1  ASP  7 

1  .’JO 

1PP7.AF 

.  Ao  i  ?;./s 

1  -  P  ..ha 

.585 15 

.13417 

P. '  0000 

PI 

.  1  ASP  7 

-'.11 

IP73.P1 

.  'OIPKS 

I  18.15 

.S7AM8 

. 1 3S IP 

P.1'0000 

PI 

,  1  4  A  -  1 

P  ,  A  7 

IPPt .-G 

. on  1 P?1 g 

1 ( OS.  IS 

.SSS  IS 

. 1 3771 

P  .  0  0  0  0  0 

P  l 

.  1  AGP  1 

A.9P 

11  A  3. PC 

.001  PACO 

SQA .A A 

.5184(1 

. 1 3SSft 

P.00000 

7  ' 

.  1  ASP  1 

1.P  1 

1  11. c  < 

.  OIPg/p 

S7  1  , n A 

. AftSHS 

. 1 39 Pa 

P  .  0  0  0  0  0 

PI 

,  1  A  S  H  7 

1  1  .SA 

0(5.48 

.  101  1 0  A  A 

SAA .SS 

.  A  11  AS 

. 1 350ft 

P  .  0  0  o  0  0 

PI 

.  1  A  G  P  7 

1  l.PS 

7  A  A  .  G  7 

.01 isis 

SI  .74 

.  34  7ft A 

. 1 PAAft 

P.  0,1000 

?1 

.  1  AS-  7 

IMS 

7 1  A,  7  A 

.  5  0  1  A  1  M 

A7' ,7A 

. ?8 1 77 

.11170 

P.00000 

PI 

.  |  A  S  P  7 

1  "  •  A  A 

API  ,  p  1 

. ; 0 1 AGl 7 

A  P  g  ,  A  7 

.PI  APS 

.09777 

P.  (’00  0  0 

p  1 

.  1  AS“7 

P'-  .  7  7 

157.4] 

.  C  0  1  s  7  G  A 

7HP.P1 

.  lftP7J 

.0801  ? 

P.10000 

?  1 

.  1  ASP  1 

P  1 . 0  K 

PSA .77 

.,’,01  A  AAA 

7  7S.SA 

.  1  1S7A 

.  (i  ft  P  3ft 

P  .  COOliO 

PI 

.  1  AS  -  7 

PS.  IP 

1  7S.  | 

.001 757  1 

7  0  A  .  S  A 

.37957 

•  0  4  S  8  a 

P  .  0  0  0  0  0 

PI 

. 1 aS-7 

■■'po 

1  P  1  .  f.A 

.'■0  1  7-) 

APP.I  7 

.CSS  3  0 

.  0  !  1  s  ? 

?  .  u  0  o  0  0 

PI 

.  1  AS"  1 

i 

11  1  .  —  A 

.•01  IP  AS 

ft  7S . 71 

.  )Aft  1  A 

.  OPo  1  1 

P  .  0  0  o  0  0 

p? 

. 1S71 7 

i  .  V. 

017.77 

.01  PSAA 

SS? .PH 

. A?ft  17 

.11981 

P  .  1  (I()l  0 

pp 

.  IS  H  1 

P  .  1  A 

SPA . SS 

.10  1  PSA  7 

/*  Ci  %  u  Q 

•  A  P  0  S  1 

.  1  IS  79 

p .  l  o  o  o  o 

p? 

.  1  S  1  1  7 

-  .  PS 

A  p  7.-1 

.,'0  1  71  '  1 

-A  1.(4 

.AC  ) AS 

. 1  ISA  3 

P.  1  0010 

p? 

.171)1 

'  .41 

-PH . A M 

.  ■  01  71. ’A 

spft.SS 

.17849 

. 1 1 7SS 

P. 1 00 00 

pp 

.  I  S  1  |  7 

",<■7 

7  A  0  .  S  A 

.1  01  1598 

Sift. 80 

.341. 7S 

,  1  l  Aft  \> 

P  .  10(11() 

pp 

.  1  s  1  1  7 

1  .71 

ASS  ,  s.j 

.01)1  is  <0 

AA1  ,  A  1 

.PSh  70 

. 1  Oh  70 

P. 1 1000 

pp 

.1-11  < 

1  P."A 

SSG.  I  7 

.  1(1  1  A  A  A  A 

AS) .PO 

.PSPSft 

.10008 

P.  1  (  OI’O 

>  ? 

.  t  s  1  1  1 

l  1  .  S' 

P  s  A  .  A  A 

.  ,  015087 

A  1  7  .  S  1 

.  P  0  ft  ft  7 

.089  11, 

? .  1  f  n  0  n 
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COMPUTATION  OF  AFMO-AOOUST IC  PPOPfrtips  OF  SWFPPFSSOP  N077LFS 


CAS  F 

NO,  1 

CPn  7-TllHF  AP=7.3  N0771  F 

-  VJ=??00  FP4  - 

TTJ=I60C 

AXIAL 

M 

L OC AT  I  ON 

P 

=  .  79167  (X/OIO  =  4  .  o  A  A  o  ) 

ANGLE  1)  Of  NS  I  T  Y  TFMP. 

tl/UPFF  TU9R.INT. 

F-/DFO 

?? 

.14313 

17.14 

347.44 

.'014777 

781 .47 

. 16743 

. a  76  74 

7.1(1000 

?? 

.14313 

I9.?9 

770.44 

.001 444  3 

744.33 

. l?3dl 

.06780 

7.1 O0O0 

7? 

. 1431 3 

71.43 

147.7? 

.'01 7344 

710,84 

.08367 

.04961 

7. 1 P0O0 

7? 

.14313 

7  3 . 4  7 

137.04 

. 501814? 

673.47 

.06,-7? 

.  0  37516 

7.10000 

7? 

. 1431 3 

74 . 7  1 

03.70 

. 00 1 P844 

444.37 

, C  4  ?6  0 

.07674 

7.10000 

?? 

.14313 

7  7.44 

44.4  7 

.0013313 

438.44 

.07331 

.0  1  7»C, 

7. 1 7000 

7? 

. 1431 3 

?  .  V  i 

43.64 

.1,0  14484 

634.1? 

.  074  3  7 

.01137 

7.10000 

73 

. 1604? 

'  .03 

413.13 

.  G 1 4 1 34 

868 . 7  r 

.77676 

.0941? 

7.70000 

73 

.  1  A  04  7 

7.14 

433. ?M 

.  ''Mi  1  4,>4? 

6  6  4.8 

.77433 

.03360 

7 . 7  0  0  0  0 

?? 

.  1  6Qi.1i 

*  .  34 

4  7  S  .  '  4 

.•-014  387 

6  S  7 .06 

.76144 

.19136 

7.70000 

?i 

.  1  Ffv 

4.43 

4  1(1.  34 

.'•014433 

843.31 

.74113 

.  1,8886 

7.70000 

73 

. 1 BC4? 

4.47 

'.77.73 

.'014374 

833.7? 

.71433 

.  064  0  1 

7.70000 

73 

.  I  A  A  4  7 

1'  .71 

404.78 

.’•0144  )3 

600.78 

.  1  8334 

.0  777(1 

7.70000 

7  3 

. 1 S047 

17. R4 

337.14 

.'■01 ‘>34? 

773.46 

.14373 

.06887 

7.7OC00 

73 

.1604? 

14.3'' 

?A9.  3  7 

.1-01  4448. 

744.48 

.  17746 

.74979 

?  •  7  A  0  A  0 

73 

. 1 604? 

1  /.  14 

?( 6.7? 

.  ‘  V 1 77  33 

714.3.' 

.  03.793 

. 04904 

7.70000 

73 

.  1604? 

1  4.?4 

147.71 

.,01 7348 

486.64 

. 063? J 

.03887 

7.70000 

73 

. 16(4? 

71.4  3 

1  7.77 

. '018481 

665.08 

,  C  4  3  0  0 

.07941 

7 . 7ti  00  C 

7  3 

.  I  A  0  4  7 

73.47 

7  3.44 

.  :013  147 

4  36.7;, 

.  0  5)  7  73 

.07169 

7. 70000 

73 

. 160*7 

?4.71 

4  8 . 7  ,3 

. CO  1 3344 

61  7.67- 

.07)3.3 

.01 884 

7.70000 

73 

. t  *04? 

77.84 

31.74 

.4U7G 74 7 

604.4? 

.01447 

.00948 

7 . 7  0  0  0  0 

7  3 

.  1  604? 

3  7.03 

74.44 

.1  070484 

60 1 . 36 

.(1161 

•  O  0  6  4  C/ 

7.70000 

74 

.16771 

T  .  0  3 

347.  to 

.  O0lf  773 

781 ,48 

.16744 

.06467 

7.30000 

?4 

.  I  <4771 

7.0  3 

34 | .4| 

.1014  4  34 

773.34 

. 14386 

.06431 

7. 300 00 

?•. 

.  1  *3  7  7  ) 

4  ,  O  A 

3  34.4-7 

. '01434 7 

7  73.74 

.14748 

.06?  7  J 

7.30000 

7*4 

. 1 4771 

4.0  0 

3  1  A  .  8 

.50161 74 

763.36 

. 14038 

.04991 

7.30000 

Ji, 

.14771 

•' .  o  r, 

777.  7 

.  '017  481 

748.  1  7 

. 17437 

.0449? 

7. 34000 

74 

.1477! 

l  1 .  or, 

740. n 4 

.(  01484) 

771.33 

.  1  C  3  1  7 

.04096 

7 . 3  0  0  0  C 

74 

. 14771 

17.74 

700.47 

.'I8177J3 

71  P.77 

.09173 

.  04404 

7. 30OO0 

74 

.14771 

1  4 . 0  c 

1  4  1  .  »  7 

.30  7 8 3 0 

631.34 

.07760 

.03844 

7. 30000 

74 

.1477] 

14.0  (j 

176.11 

.4018347 

671 . 77 

. 04  7  54 

.03189 

7 . 3  0  0  C  0 

74 

.14771 

14.CA 

84 , 74 

.5018  V) 

'  c  1.41 

.74318 

.07440 

7. 300  CO 

7  4 

.14771 

7  '.Cl 

4  8  ,h7 

.'013434 

473,46 

.0717? 

.01344 

7. 30000 

74 

.14771 

77.7 

4  7.31 

.  '070',-  33 

414. 74 

.071  78 

•  (1 1  446 

7 . 3  0  0  0  0 

7-4 

.14771 

p  ,4  ^  j3  ( 

37.7  7 

.0070431 

6  0  0.6-1 

.  J 1 46  7 

.  (1  1  J  3  I 

7. 3<'00o 

7  4 

.14771 

7'  .GO 

70.0  7 

•  5  U  ?<..  34  7 

488,64 

.Of  944 

.00  709 

7 . 3(1080 

34 

.14771 

13.73 

.' 071 7 13 

^8)  .13 

.00674 

.00464 

7.30000 

;•  4 

. 14771 

3 ,  .  c  r. 

1  0.8~ 

.■,-071,-33 

>->78 .34 

.00639 

.00777 

7 . 3(i 0 10 

3  4 

.  1  7  V,  r 

-  .  30 

186.1? 

.-017  .34 

7  04.8  1 

. 0846? 

. 7  7974 

7.4 1000 

■>  c. 

.  1  7B(\ 

7.  ;ii 

1  87.4,4 

. .01 7438 

7'  3.08 

.  08  56; 4 

.0 3877 

7.4(1000 

3  4 

.  1  7  V 

**  •  C  v 

173.41 

.  017441 

638 . 3  , 

.07884 

.77761 

7.44000 

P6 

.  1  74  A 

'  •  ? 

lr-9.  ]  4 

.  ’01  78-80 

630.73 

.07736 

.(,34  3? 

7.40000 

.  1  7V  -i 

“  .  3  0 

1  4  0 . 4  3 

.•'0181  (4 

673.03 

•06394 

. 3  17  37 

7 . 4  0  0 1 0 

1  L 

.  1  7',:. 

1  o  .OC 

•  C  016448 

667.68 

.  C  4  4  6  0 

. 078 78 

7.40000 

3  c. 

.1  7^  O  A 

1’." 

33.  *  4 

.'>0  1804', 

6  44 , 14 

.  044  a  7 

. 07488 

7.4 Apoo 

•  I  t 

.  1  74  V,- 

1  4  .  A  , 

7  8.44 

.1  013374 

633.77 

.1-7466 

•  0  ?  0  7  C 

7.4.1000 

1  4 

.  1  74.  .4 

I'.J. 

43. -.1 

.  013714 

634.43 

.07774 

.01667 

7.40000 

V. 

.1  7411 

1 

4  4.1  1 

.  1  4  3 

611.86 

.  (1 7 1, 116 

.  :1  1  7  3  1 , 

? .  4  a  n  ,1 0 
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COMPUTATION  of  afro-acoustic  PROPERTIES  OF  SUPPRESSOR  nozzles 
CASE  NO.  1  CRH  7-TURF  AP=2.3  N07/LE  -  VJ=2?00  FPS  -  TTJ=1600  OEF— R 


AXIAL  I.  Of  AT  ION  =  .29167  (X/OEO  =  4.000'j)> 


M 

P 

ANGLE 

U 

DENSITY 

TEMP. 

U/UREF 

TURR. I  NT 

.  R/r;EO 

?s 

. 17SO0 

a  ■> .  L  o 

3  1  .  1  A 

a  3  0  A  0  S  A  3 

SPP.OP 

.C141H 

.00964 

2.40000 

AS 

. 1 7so: 

PP.hr. 

A  1  .  1  A 

.3020967 

SR7.5S 

.  0  0  9  6  0 

.00694 

2.40000 

?s 

a  1  7  S  0  0 

Ph  a  n  :• 

!  3.  77 

. *021 3P4 

S77.7I 

•  '  0  6  A  A 

.004P1 

2.40000 

?S 

.  1  7S  0  0 

Ph  .  0  (, 

A  ,  hh 

.0G21A3A 

SEP  .  PA 

.00394 

.00334 

2.40000 

PS 

a  1  7S00 

?  a.  or 

5. Pi 

,  0  A  1  A  A  A 

564  a  PA 

.00246 

.  0023A 

2.40000 

PS 

a  1  7G0C 

30.00 

P.18 

.0021677 

SA3.A4 

.00190 

.0019? 

2.40000 

Ph 

.  1  A  2A9 

7  a  A  p 

.TO  IPS  PA 

A4R. A7 

.r3p40 

.02061 

2.50000 

Ph 

a  1  A  ??0 

1  .  A7 

AS.  -P 

a  0  0  1  p  1  3  1 

a44 . 56 

.  1.  3H3l9 

.02GS2 

2.50000 

Ph 

.  )  RAA 9 

t.  7S 

A  0  .  PP 

.(  0  !  PA  1  7 

A4) ,A4 

.  (.  3  AH  0 

,  0  1  9  7  4 

2.50000 

?h 

. 1 AP?R 

S.hP 

7P  a  AP 

a  EC  IP  AS  7 

A37.0? 

.  0  3  3  9  A 

.01666 

2.56  000 

Ph 

a  1 

/.  SC 

AA.P  1 

.  P  0 1  PRS A 

A3.-  ,A7 

.03,19 

.0  1697 

2 .50000 

Ph 

a 1»2AQ 

9.37 

S7.  IS 

.1  C1P7AA 

AA3.  IP 

.  0  2  S  9  R 

.OIS) 0 

2.50000 

Ph 

a  I  H ?PQ 

I  l  .AS 

P  7  .  A  S 

.  <  flAOOPA 

615.45 

a  C  2  1  A  A 

.01303 

2.50000 

Ph 

a  1  6?A9 

n.  i  a 

3«.  AP 

a"C AC  337 

f'Od.33 

.01/41 

.01067 

2.50000 

Ph 

a l"??0 

15.33 

AP  .  AS 

.  0?O  AP  1 

SP7.4  i 

.01  34  H 

•  0OR79 

2.5O0O0 

Ph 

a  1  6??9 

1  A  .  M  7 

A  A.  I  7 

. : 0  AC  PP  A 

S  A  H  .  R  1 

.Oil  OR 

.00693 

2.50000 

Ph 

a  1  R2A9 

1  8 . 7  R 

IN’  3 

.  r  0  A 1  A  A  P 

S«  .1 . 44 

a  0‘  0  7  2  P 

. 0  05  74 

2.50000 

Ph 

a  1  R  ?  ?  9 

Pj.hp 

11  a  IP 

a  '■  0  2  1  S  ('•  1 

S73„Sv; 

.006O7 

.00396 

2  a  5  0  0  0  0 

Ph 

. 1  R??9 

aa.s; 

7  .  PP 

..  ii  A  1  7  P  A 

SA  7 , OS 

.  0  0  3  3  R 

.00293 

2.50000 

Ph 

a  1  A PPS 

?P  a  .37 

4.7a 

."0A1PSS 

SAl .63 

.'0216 

.00224 

2.50000 

Ph 

a  1  6AA9 

A  A  .  A  R 

A.  a  a 

a^uAAl |a 

SS7.SA 

.'0131 

.0017? 

2. 6O0O0 

Ph 

.  1  A PPS 

A  8  a  1  A 

1  afcA 

0A?1 ap 

SGS.7A 

a  0  0  0  7  6 

.00IS1 

2.6C0D0 

Ph 

. 1 PPPS 

3  0  a  0  f 

1  a  IP 

,  '  (i  A  A  t  7  P 

S  S  A  .  A  R 

a  0  0  0  S  4 

.0013? 

2,6(1000 

PI 

. 1 AOSA 

A  a  C  0 

3  6  .  C  A 

.  .  C  2  C  4  1  1 

t  -'4.1  A 

.01640 

.00971 

2.60000 

PI 

.  1  ARC W 

1  a  A  7 

3S.  A  A 

. 702343A 

Ac  3.37 

.  C  1 6  0  R 

.  0  0  9  S  7 

?  .  6  0  0  0  0 

PI 

. 1 H9RS 

3. 7R 

33. RA 

.COACS.O 

AO  I .SO 

.01S24 

.00917 

2.60000 

PI 

. 1 R9RR 

S.hP 

30.7  A 

a  O  2  0  A  (  0 

SQp.RQ 

a  C  1  3  9  7 

.O1.HSR 

2.60000 

PI 

. 1  8  95  8 

7.5  r 

A  7  .  1  A 

a  r  0?  0  7  3 A 

SPp . 64 

.01233 

.00776 

2 . 6  0  0  0  0 

PI 

.  1  APRS 

P.37 

A3.'  A 

.  •  ■  o  ?  o  p  •:  s 

S«9.H4 

.01046 

.00661 

2.60000 

PI 

a  1  ARRP 

1  1  .AS 

l  A  .  P  3 

.T0A10A7 

SP4.7S 

. 00661 

.00567 

2.60000 

PI 

. 1 HRRA 

1  3.  1  A 

IP. PA 

.3  0  A 1  A  a  S 

S7P.  3 A 

. 00680 

.004R5 

2 . 60  0  0  0 

PI 

.  1  A9RA 

1  A  a  C  0 

1  1  a  3P 

.  •'OAI  4RA 

S  7  3 . 9 

.00S1  R 

.0039 A 

2.600O0 

P  7 

.  1  APR A 

1  A. 87 

A  a  3  <; 

a  ■'OAl  AMR 

568. A? 
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4.3  PROGRAM  SOURCE  CODE  LISTING 


This  section  contains  the  FORTRAN  IV  source  code  listing  for  the 
aeroacoustic  prediction  model,  suitable  for  running  on  the  CDC  7600  computer. 
The  listing  of  subroutines  is  in  alphabetical  order,  as  follows: 


1. 

MAIN  Program  (MGB) 

2. 

ARRCCOS 

3. 

ATMOS 

4. 

CRD 

5. 

ERF 

6. 

LSPFIT 

7. 

OUTPUT 

9. 

PNLC 

9. 

SHOCK 

10. 

SI.TCF. 

11. 

TPNI.C 
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5.0  CONCLUDING  REMARKS 


Two  computer  programs  capable  of  predicting  the  jet  noise  of  high  veloc¬ 
ity  exhausts  from  nozzles  of  arbitrary  geometry  are  presented.  The  com¬ 
puterized  procedures  presented  herein  provide  reasonably  accurate  methods  of 
predicting  maximum  sideline  PNL  as  well  as  EPNL  (with  and  without  flight 
effects)  over  the  range  of  flow  conditions  and  observer  angles  of  interest. 
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